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CHAPTER  I 
Introduction 

1.  Purpose ♦ 

It  is  the  purpose  of  this  report  to  provide  information 
on  the  shipping  densities  in  1°  x 1°  squares  throughout  the  world. 
This  report  explains  the  techniques  and  procedures  by  which  these 
densities  were  generated.  The  data  generated  in  meeting  this 
objective  is  stored  in  machine  readable  format  and  is  available 
from  Planning  Systems  Incorporated. 

2.  Scope . 

This  report  describes  in  rather  substantial  detail  the 
procedures  by  which  shipping  densities  in  1°  x 1°  squares  are 
available.  The  basic  data  base  which  was  employed  was  the  1972 
Lloyd's  Shipping  Index  indicating  the  disposition  of  18,000  ships 
per  week  in  the  2,425  Lloyd’s  ports.  The  1972  data  base  was 
updated  to  1977  using  a trend  and  growth  analysis  which  is  dis- 
cussed in  some  detail  in  the  report.  Five  different  types  of 
vessels  are  considered  in  the  analysis.  These  are  merchant 
vessels,  tankers,  super-tankers,  fishing  vessels,  and  oil  rigs. 

The  distributions  of  the  individual  types  of  vessels  are  given  for 
each  month  of  the  year,  each  quarter  of  the  year,  and  an  annual 
average.  This  results  in  17  distributions  per  type  of  vessel. 

The  report  describes  the  methodology  for  distribution  of  the 
updated  port  traffic  throughout  the  world.  It  is  heavily  depen- 
dent upon  a preceding  volume  issued  by  Planning  Systems  Incorpor- 
ated "Ocean  Route  Envelopes",  (ORE)^  published  on  19  April,  1977.  The 
resulting  distributions  consist  of  well  over  12,000,000  individual 
pieces  of  data.  Due  to  this  enormous  data  base,  it  is  impossible 
to  present  all  distributions  and  possible  combinations  in  a printed 
format  within  the  body  of  this  report. 
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Therefore,  only  examples  are  presented.  However,  any  dis- 
tribution or  combination  of  distributions  may  be  obtained  in 
machine  readable  format  directly  from  Planning  Systems  Incorpor- 


ated. 


Background. 


The  primary  purpose  of  this  project  was  to  quantify  the 
distributions  of  all  ships  in  the  oceans  of  the  world  at  various 
times  of  the  year.  One  of  the  necessary  tasks  to  accomplish  this 
objective  was  to  employ  and  develop  seasonal  route  envelopes  which 
would  allow  the  movement  of  ships  between  the  world's  ports.  The 
routes  themselves  are  subject  to  a very  detailed  development  and 
have  been  presented  in  a previous  report:  Ocean  Route  Envelopes 

published  by  PS I. 

World  shipping  densities  are  seasonally  dependent  and  are 
due  to  tv>  ?•  Imary  factors: 

1.  The  seasonal  variation  of  the  routes  by 

which  ships  move  between  world  ports  . 

2.  The  seasonal  variation  in  port  traffic. 

The  first  part  of  this  problem,  addressed  in  the  ORE  report, 
dealt  exclusively  with  the  logic  and  mathematical  quantification 
of  the  development  of  the  temporal  world  shipping  routes.  The 
second  part  of  the  problem  i3  the  world  port  activity.  This 
activity  has  been  provided  through  a numerical  quantification 
of  1972  Lloyd's  Shipping  Index.  The  index  gives  the  locations 
of  approximately  18,000  ships  per  week.  In  the  preparation  of 
this  report  it  was  necessary  to  update  this  1972  data  to  the  year 
1978.  A number  of  techniques  were  employed  based  upon  industry- 
wide sources.  Therefore,  the  shipping  densitites  are  considered 
to  be  1978  information.  This  report,  Historical  Temporal  Shipping 
(HITS),  discusses  the  existing  Lloyd's  data  base,  the  trend  and 
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growth  analysis,  the  distribution  methodology  to  establish  the 
shipping  densities,  and  the  type  of  outputs.  In  addition  to  the 
data,  it  provides  information  on  the  development  of  the  computer 
programs  which  were  necessary  in  developing  such  a large  and  exten- 
sive data  base. 
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CHAPTER  II 

The  Lloyd  * s Data  Base 


1.  Logical  Structure  of  Data  Base. 

The  DOT  Marine  Traffic  Generator  Input  Tapes  provided  the 
majority  of  the  raw  data  on  ship  movements  necessary  to  generate 
the  shipping  density  for  this  study.  These  tapes  were  originally 
prepared  by  the  Department  of  Transportation  for  use  by  the  Marine 
Traffic  Generatorf  and  included  data  from  Lloyd's  Shipping  Index 
for  1972. 

These  tapes  were  used  in  the  study  by  Program  Gamma  to 
accumulate  the  total  time  spent  by  ships  passing  between  each 
pair  of  ports  m the  DOT  Marine  Traffic  Generator  port  table. 

Later  programs  then  converted  this  information  into  shipping  den- 
sity fields. 

The  Lloyd's  Shipping  Index  contains  a weekly  record  of 
movements  of  many  merchant  vessels  (approximately  18,000).  It 
records  the  egress  port,  entrance  port,  and  port  of  last  call 
for  each  ship,  along  v ith  the  date  of  egress.  An  additional 
date  field  may  fill  one  of  several  purposes,  depending  upon  the 
setting  of  various  flags.  Possible  interpretations  are  date  of 
arrival  at.  entrance  port,  '■''ate  of  arrival  at  port  of  last  call, 
date  of  departure  from  port  of  last  call,  date  of  passing  port 
of  last  call,  or  dace  of  presence  in  the  port  of  last  call. 

Each  ship  is  identified  by  name,  flag,  gross  tonnage,  year  built, 
and  type  (tanker  or  non*  tanker) , and  is  recorded  each  week  in 

9 

approximately  the  same  position  within  the  Shipping  Register. 

Each  of  the  52  tapes  contains  movement  reports  received 
during  the  corresponding  week  of  1972.  However,  the  actual  date 
of  the  ship  movements  recorded  can  be  as  much  as  several  months 
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prior  to  the  date  of  report  itself.  Frequently,  the  same  data 
will  be  repeated  in  several  successive  weekly  tapes  without  change. 
Furthermore,  as  the  length  of  a voyage  can  be  appreciably  longer 
than  the  weekly  reporting  intervals,  no  one  tape  will  necessarily 
contain  sufficient  data  to  permit  reconstruction  of  the  ship's 
route  and  travel  time.  At  the  opposite  extreme,  only  one  voyage 
with  a maximum  of  three  ports  (only  two  with  dates)  could  be 
included  per  ship  per  week,  so  that  numerous  short-haul  trips  were 
completely  ignored.  This  latter  phenomenon  was  aptly  illustrated 
by  ships  that  were  recorded  as  having  arrived  in  port  in  one 
weekly  issue  and  then  as  departing  from  a different  port  in  the 
following  week's  issue,  with  no  intermediate  voyage  recorded. 


A possible  fundamental  weakness  in  the  data  base  used  arises 
from  the  selective  nature  of  the  reporting.  Lloyd's  was  almost 
completely  dependent  for  its  data  upon  sailing  reports  from  the 
various  harbour  masters.  The  completeness  of  the  data  thus  depends 
in  most  cases  upon  the  diligence  of  a single  reporting  source  and 
the  completeness  of  his  information.  Numerous  reports  exist  in 
the  shipping  index  of  ships  sailing  off  to  nowhere  or  arriving 
from  points  unknown.  There  are  no  checks  against  consistent  under- 
reporting of  traffic  from  a specific  port  and,  in  particular, 
sailing  reports  for  Communist  owned  shipB  between  most  Communist 
ports  are  conspicuously  absent.  Date  information,  too,  is  fre- 
quently missing.  This  weakness  in  the  data  collection  is  reflected 
by  the  incompleteness  of  the  final  shipping  index.  The  number  of 
ships  reported  each  week  (approximately  13,000)  is  much  smaller 
than  the  total  number  of  elligible  ships  from  Lloyd's  and  no 
one  ship  is  necessarily  monitored  continuously  throughout  the 
course  of  the  year. 


5 


KHZp*?rm  ?*!■*,  (.'-V'TCVTC'.'*  n-Vt' 


In  order  to  extract  the  maximum  amount  of  information  from 
Lloyd's  shipping  index,  it  was  necessary  to  correlate  corresponding 
reports  from  successive  weekly  issues.  This  was  done  by  matching  the 
ship  reports  on  the  basis  of  the  identifying  information  (name,  weight 
type,  flag,  and  year  built).  The  reported  ship  movements  could  then 
be  compared  to  eliminate  overlapping  reports,  supplement  incomplete 
information,  and  simulate  some  unreported  voyages.  As  movement 
reports  (except  for  expected  destination)  were  normally  received 
after  the  movements  reported,  the  tapes  were  processed  backwards, 
beginning  with  the  shipping  index  tape  for  the  last  week  of  December 
1972,  and  proceeding  backwards  to  January  of  that  year. 

The  nature  of  the  data  is  illustrated  by  the  first  50  records 
from  the  DOT  input  tapes  for  weeks  51  and  52  shown  in  Table  2-1. 

Some  sample  match-up  procedures  from  these  lists  are  shown  in 
Table  2-2.  An  examination  of  the  missing  data  list  in  the  notes 
to  Table  2-1  reveals  122  missing  or  antiquated  items  of  data  among 
the  100  reports,  for  an  average  of  1.21  missing  or  antiquated 
data  per  report.  Quite  clearly,  it  is  only  the  exceptional  report 
which  is  complete  and  up-to-date. 
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Table  2-1  Sample  Listing  for  DOT  Input  Tapes 


DOT  INPUT  TAPE  FDR  WEEK  51,  1972 


A 

ft.  E.  S.  I 
ft.  HELD  a 
ft.  N.  KEMP  T 
ft.  P.  HOLLER  * 
ftftSE  ftNDRERSEN 
ftftP 

ftftSE  HOJ 

RBRSURLES 

RBRKRNLES 

ftBftSIN 

ftBftVft 


BCD  E F 

0 Bft  57  001579  ROTTE 

0 GE  67  000211  KOTKft 

1 LI  50  013127eNEW  V 
1 Dft  66  052673cPHILft 
0 Dft  66  000499  WESTE 
0 GE  32  000330dHftRTL 
0 Dft  55  000300  HULL 
0 RU  62  0046526G9RST 
0 RU  61  004633  ftNTWE 

0 PHzr64  003989-  BftNGK 

1 RU  71  OG3469fMftLMO 


G H I J K 

DEC  11  LONDO  ftR  DEC  12 

DEC  05  MftftST  ftR  ftNTWE  DEC  14 

DEC  06  ftRUBft 
NOV  295ESCRft 

NOV  175 TORRE  ftR  DEC  03 

DEC  09  BREME 

DEC  14  GUNNE  ftR  DEC  14 

SEP  225ftRCHft 
DEC  09  LENIN 

2KRRftC  ftR  KftRftC  JUL  25 

DEC  13  KLftlP 


ftBBOTSFORD 

1 

BR 

72 

001584c HftVRE 

DEC 

08  DftKftR 

ftR 

DEC 

16 

ftBEL  TRSMRN 

0 

DU 

71 

027614 

MELBO 

NOV 

145LDNDO 

ftR 

DEC 

16 

ftBERTHftW  FISHER 

0 

BR 

66 

002355 

ROTTE 

NOV 

295ftSHD0 

ftR 

DEC 

13 

4 

ftBIDft 

1 

DU 

53 

012226dSTftNLD 

DEC 

06  CURftC 

ftBIDIN  DftVER 

0 

RUH60 

004449CMERSI 

DEC 

03  ISTftN 

ftR 

DEC 

06 

ftBIDJftN 

0 

LI 

67 

003043 

DUftLft 

DEC 

07  BREME 

ftR 

ftBIDJ 

DEC 

10 

4 

RBIQUR 

0 

ftM 

43 

014640^CHITT 

OCT 

2556ftLVE 

ftR 

CORPU 

DEC 

08 

4 

ftBLE  REEFER 

0 

SG 

61 

Q02683CKRIST 

2TR0ND 

ftR 

DEC 

11 

ftBLON 

0 

FR 

61 

006451 

HftVRE 

DEC 

01  BUENO 

SB 

PftLMft 

DEC 

12 

4 

ftBOCOL 

0 

BR 

40 

001834eCftRTft 

NOV 

225TftMPft 

ftBOISSO 

0 

IV 

55 

004939 

ftMSTE 

2DUflLft 

SD 

ftNTWE 

DEC 

16 

4 

ftBOTft 

0 

CY 

58 

000499 

IZMIR 

NOV 

23  BREME 

ftR 

DEC 

14 

ftBRENE 

1 

RU 

64 

003390cLDNDD 

DEC 

14  BftLTI 

4 

ftBU  SIMBEL 

0 

EG 

60 

001946 

LftTTft 

2 ftLEXft 

ftR 

DEC 

09 

RBUGRWR  MftRU 

0 

Jft 

65 

001732c  i 

2 3 

ftR 

SINGft 

DEC 

03 

ftCftDIR  FOREST 

0 

NO 

69 

036362 cROTTE 

DEC 

12  NEW  0 

SD 

SHEER 

DEC 

13 

4 

ftCftMftR 

1 

LI 

59 

016308cCftRftC 

DEC 

12  3 

ftCftPULCD  MftRU 

0 

Jft 

51 

006493 

KOBE 

NOV 

055TDWN 

ftR 

REUNI 

DEC 

14 

4 

ftVRCUS B 

1 

BR 

53 

012326^SftRNI 

DEC 

07  ftRUBft 

PD 

ESCOU 

DEC 

13 

4 

RCHRIOS 

0 

LI 

50 

005455 cTftRRft 

DEC 

12  KftlSE 

SD 

CEUTft 

DEC 

15 

4 

ftCHftT I Nft 

1 

BR 

58 

012326 

HftMILO 

2 3 

SD 

LftNOR 

DEC 

14 

ftCHILLE  LftURO 

0 

1^47 

023629 

ftR  E Nft 

NOV 

295GENDft 

SD 

MftLTft 

DEC 

13 

4 

RCH1LLES 

1 

ftM 

60 

02447 leFRftNC 

DEC 

10  fttftSKfi 

4 

ftCHILLES 

0 

BR 

72 

016406cOSflK.ft 

OCT 

075  3 

4 

ftCHILLES 

0 

DU 

59 

006134 

DEMER 

DEC 

03  BREME 

SD 

MICHft 

DEC 

16 

4 

ftCHILLES 

0 

GE 

64 

000500 

ROTTE 

DEC 

14  LISBD 

4 

RCHILLET 

1 

GR 

56 

031571 

KftRLS 

DEC 

12  MftRSftB 

4 

RCRRDINR b 

1 

IT 

60 

022955 

UMtl  S 

DEC 

13  GENOft 

4 

BCHILLEUS 

0 

CY 

39 

002920 

i 

2 2 cftR 

PIRftE 

OCT 

08 

99 

ftCHILLEUS 

0 

GR 

65 

021505 

NftGOY 

DEC 

02  ftRGEN 

ftR 

MURDR 

DEC 

18 

4 

ftCILft 

1 

DU 

58 

012221 

MftNIL 

DEC 

09  CURftC 

4 

ftCMREft 

1 

DU 

59 

012222 

ROTTE 

DEC 

17  FREDE 

4 

RCDNCftGUft 

0 

CH 

65 

010869 

JOHN  N 

DEC 

13  3 

ftR 

PHILft 

DEC 

17 

4 

ftCONCRGUft 

0 

FI 

60 

005402 

HELSI 

NOV 

22  SCftSftB 

ftR 

DEC 

14 

RCONCftGUR  VftLLEY 

0 

SW 

63 

009611 

HftMBU 

DEC 

16  CftSftB 

4 

ftCQRES 

0 

PD 

58 

002963 

MOCftM 

NOV 

30  5BREME 

ftR 

LISBD 

DEC 

15 

4 

RCQUftRlUS 

7 


1 PR  58  020519°EURQP  DEC  18  BRNDR 
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Tabic  2-1  Sample  Listing  for  DOT  Input  Tapes  (contd) 


DDT  INPUT  TAPE 

FDR  1 

JEEK  52  > 1972 

A 

B 

c 

D 

B F 

G 

H 

I 

J 

K 

L 

A.E.S.  1 

0 

DR 

57 

001579  LONDO 

DEC 

20 

3 

4 

A. N. KEMP  1 

1 

LI 

50 

01812?®ARUBR 

2 

NEW  Y 

AR 

DEC 

£0 

fl.P.MQLLER  I 

1 

DA 

66 

052673 CPHILA 

NOV 

295ESCRA 

RAGE  ANDREASEN 

0 

DA 

66 

000499eTDRRE 

‘MANTY 

SD 

CDPEN 

DEC 

184 

AAR 

0 

GE 

32 

000330dGDTHE 

DEC 

21 

OLDEN 

4 

BASE  HOJ 

0 

DA 

55 

000300  GUNNE 

DEC 

15 

AALBO 

4 

ABAGURLES 

0 

RU 

62 

004652  c&  Aft  Jf 

SEP 

22  5RRCHA 

4 

ABAKANLES 

0 

RU 

61 

004638  BNTWE 

DEC 

09 

LENIN 

4 

ABAS IN 

0 

PKH&4 

003939  BANGK 

2 KARAC 

A 

RBAVA 

1 

RU 

71 

003468f GOTHE 

DEC 

£0 

KLAIP 

4 

ABBOTSFORD 

1 

BR 

72 

001534c HAVRE 

DEC 

03 

DAKAR 

AR 

DEC 

16 

ABEL  TASMAN 

0 

DU 

71 

027614  LONDO 

DEC 

£2 

MELBO 

AR 

HAMBU 

DEC 

23  4 

ABERTHAW  FISHER 

0 

BR 

66 

002355  ASHDO 

DEC 

£0 

3 

RBI  DA 

1 

DU 

.58 

0122264 ANNA 

DEC 

18 

SYDNE 

PD 

PANAM 

DEC 

22  4 

ABIDIN  LAVER 

0 

Tl/ 

*0 

004449cMERSI 

DEC 

03 

1ST  AN 

AR 

DEC 

06  5 

ABIDJAN 

0 

LI 

67 

003048  DUALA 

DEC 

07 

BREME 

SD 

RBIBJ 

DEC 

14  4 

RBI QUA 

0 

AM 

43 

014640fCORPU 

2 

EVERG 

SD 

FREEPB 

DEC 

20  4 

ABLE  REEFER 

0 

Sto 

61 

0 02683 cKR 1ST 

2 

TROND 

AR 

DEC 

11 

ABLON 

0 

FR 

61 

006451  HAVRE 

DEC 

01 

BUENO 

SD 

SALVA 

DEC 

£2  4 

ABOCOL 

0 

BR 

40 

001334eTAMPA 

2 

CARTRC 

AR 

DEC 

16  „ 

ABOISSO 

0 

IV 

55 

004939  RMSTE 

2 

DUALA 

AR 

DUNK  1 

DEC 

16  4 

ABOTA 

0 

CY 

53 

000499  BREME 

DEC 

16 

BENSH 

4 

ABRENE 

1 

RU 

64 

00336  (FcLONDO 

DEC 

14 

BALT  I 

4 

ABU  SIMBEL 

0 

EG 

60 

001946  ALEXA 

DEC 

22 

LATTA 

4 

RBUGRWA  MARU 

2 

AL 

50 

01732  c i 

2 

3 

AR 

SINSR 

DEC- 

03  4 

ACADIA  FOREST 

0 

NO 

69 

036862 cROTTE 

DEC 

12 

NEW  0 

SD 

SHEER 

DEC 

13  4 

ACAMAR 

1 

LI 

59 

0168Q8cCARAC 

DEC 

12 

3 

ACAPULCO  MRRU 

ACAVUS  ii 

ACHAIOS 

ACHAT I NR 

flCHILLE  LAURQ 

ACHILLES 

ACHILLES 

ACHILLES 

ACHILLES 

ACHILLET 

ACHILLEUS 

FlCHILLEUS 

AC  I LA 

ACMAEA 

ACONCAGUA 

ACONCAGUA 

ACONCAGUA  VALLEY 

ACORES 

ACQUARIUS 

ACRADINA 

ACRITAS 


0 JR  51  006493  KOBE 

1 BR  58  012326fSRP.NI 

0 LI  50  005455 cTARRA 

1 BR  58  012326  HAMILD 

0 IW7  023629  ARENA 

1 AN  60  02447  leRNCHO 
0 BR  72  01 64 06c OSAKA 
0 nu  59  006134  BEMER 

0 GE  64  000500  ROTTE 

1 GR  56  031571  KARLS 

0 CY  39  002920  * 

0 GR  65  021505  NAGDY 

1 DU  58  012221  MANIL 
1 DU  59  012222  ROTTE 
0 CH  65  010869  JOHN  N 
0 FI  60  005402  HELSI 

0 SW  68  009611  HAMBU 

0 PO  53  002963  MOCAM 

1 PA  58  02051 9CEUR0P 
1 IT  60  022955  UMM  S 
0 GR  71  009368  DULUT 


NOV  QS^TOWN  SD  REUNI 
DEC  07  ANNA  AR 
DEC  12  KAISE  SB  CEUTA 
2 3 pc  ESCDU 

NOV  E95 GENOA  AR 
DEC  12  LOS  A 
OCT  07  5 3 

DEC  08  BREME  AR  BNTWE 
DEC  14  LISBO  AR  LEIXO 
DEC  12  MRRSRB  SD  CORUN 
2 3 CRR  PIRftE 

DEC  02  ARGEN  AR  NUROR 
DEC  09  CURAC 
DEC  17  COPEN  AR 
DEC  13  RRICft  AR  NEW  Y 
NOV  223TRIP0  SD  LA  GD 
DEC  16  CASAB 
NOV  3 05 BREME  AR 
DEC  18  BANDA 
DEC  13  GENOA  SD  DAS  I 
DEC  08  3 PL  TRACY 


DEC  19  4 
DEC  21 
DEC  15  4 
PEC  15  4 
DEC  20 


DEC  22  4 
DEC  20  4 
DEC  22  4 
OCT  08 
DEC  18 

DEC  20 
DEC  19  4 
DEC  22  4 


LISBO  DEC  15  4 


DEC  15 
DEC  17 


4 
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Notes  to  Table  2-1 

(Samples  from  DOT  Input  Tapes  for  Weeks  51  and  52) 

Column  Headings 

A Ship 1 s name 

B Ship's  type  ( 1 » tanker,  0 = non-tanker) 

C Ship's  flag  (abbreviated) 

D Ship's  year  built  (last  two  digits) 

E Ship's  weight  (in  gross  tons) 

F Port  of  origin 

G Date  of  departure 

H Port  of  destination 

I Code  determining  meaning  of  column  (k) 

J Port  of  last  call 


Code  IV 
column (I) 

Contents  of 
column (s) 

Meaning 

SD 

used 

date  of  departure  from 
port  of  last  call 

PD 

used 

date  when  port  of  last 
call  was  passed 

AR 

used 

date  of  arrival  at  port 
of  last  call 

AR 

empty 

date  of  arrival  at 
destination  port 

IN 

used 

ship  was  in  port  of  last 
call  on  reported  date 

L Tragedy  code 


Incomplete  and  Difficult-to-Handle  Data 


.1 

2 

3 

4 

5 

6 


Mo  departure  port 
Mo  departure  date 
No  destination  port 
No  destination  arrival 

Antiquated  date 

Generic  port  name,  not 
precisely  identifiable 


(4  cases  in  samples) 
(17  cases  in  samples) 
(14  cases  in  samples) 
date 

(64  cases  in  samples) 
(22  cases  in  sample) 

(1  case) 


Correspondence  Difficulties  Between  Weeks  51  and  52 

a Ship  omitted  in  records 

from  other  week  (1  case  in  sample) 

b Ship  out-of  place  (1  case  in  sample) 

c Duplicated  data,  no  change 

(12  cases  in  sample) 

d Unrecorded  intermediate 

voyage  (2  cases  in  sample) 

e Departure  or  arrival  date 

must  be  interpolated  (4  cases  in  sample) 

f Ambiguous  route  (3  cases  in  sample) 

..ftiaspsUed-jand  JBgtaasaaa.  Data 

I.  Ship's  name  spelled  both  with  and  without  blanks 

II.  Other  divergences  between  data  for  successive  weeks 
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Table  2-2  Illustrative  Examples  of  Ship  Match-up 

from  Successive  DOT  Tapes  (for  weeks  51  & 52) 


ACHILLES 

0 

DU 

59 

006134 

DEMER 

DEC 

03 

BREME 

SD 

MICH  A 

DEC 

16 

ACHILLES 

0 

DU 

59 

006134 

BEMER 

DEC 

03 

BREME 

AR 

AMT  WE 

DEC 

£2 

ABUTA 

0 

CY 

59 

000499 

IZMIR 

MOV 

38 

BREME 

AR 

DEC 

14 

ABOTA 

0 

CY 

59 

000499 

BREME 

DEC 

16 

BENGH 

ABDCOL 

0 

BR 

40 

001834 

CARTA 

.NQV 

2£ 

TAMPA 

ABOCOL 

0 

BR 

40 

001334 

TAMPA 

CARTAC 

RR 

DEC 

16 

Same  voyage:  Departure  and  destination  ports  are  the  same 

from  both  week's  reports.  Merging  data  may  permit  missing 
dates  and  other  data  to  be  filled  in,  and  gains  routing 
information  through  successive  port-of-last-call  reports. 

Two  voyages:  Destination  port  from  weeks  51  and  departure 

port  from  week  52  match.  Merging  data  from  successive  weeks 
allows  for  disjoint  time  interval  of  voyage,  and  reporting 
interval,  and  may  permit  interpolation  of  missing  dates. 

Three  voyages:  Destination  port  from  week  51  does  not  match 

any  port  from  week  52.  An  intermediate  unrecorded  voyage 
has  occurred,  which  may  be  plausibly  reconstructed. 
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Technical  Structure  of  The  Data  Base. 


The  52  DOT  Marine  Traffic  Generator  Input  Tapes  were  trans- 
scribed  by  hand  from  the  Lloyd's  publications.  After  having  been 
initially  keypunched,  the  cards  are  edited  anc  copied  onto  52 
reels  of  7-track,  unlabeled  magnetic  tape  at  a recording  density 
of  800  BPI  and  even  parity.  Data  is  recorded  in  external  BCD 
code  in  unblocked,  84-byte  fixed-length  records,  with  data  fields 
as  specified  in  the  original  DOT  specification  (Marine  Traffic 
Generator  Documentation,  p.25)  and  as  illustrated  in  Table  2-1. 

Although  originally  prepared  for  use  vdth  an  IBM  7094  com- 
puter, these  tapes  are  fully  compatible  with  most  data  processing 
systems,  including  the  IBM  360/370  series.  Much  the  same  data  is 
duplicated  on  the  52  Marine  Traffic  Generator  output  tapes;  however, 
these  are  much  harder  to  utilize,  as  they  are  recorded  in  an  inter- 
nal 7094  format  and  use  a record  format  no  longer  recognized  by 
IBM.  Furthermore,  having  been  processed  previously  by  the  Marine 
Traffic  Generator,  these  tapes  could  not  be  reused  without  risking 
accumulation  of  errors  through  successive  layers  of  programming. 

These  input  tapes  were  edited  by  two  DOT  programs  whose 
purpose  was  to  identify  and  eliminate  major  keypunch  errors.  In 
particular,  data  punched  in  the  wrong  field,  gross  date  errors, 
and  unknown  or  misspelled  poet  names,  should  have  been  detected 
and  corrected  by  these  programs.  The  tapes  used  in  this  study 
represent  the  output  of  these  editing  programs. 

Nevertheless,  many  simple  errors  remain  in  the  tapes. 

These  include: 

1.  Variations  in  spelling  of  ship's  name,  particularly 
insertion  or  emission  of  blanks  after  abbreviations, 

2.  Right-or-left  shifted  data,  i.e. , *01782  * ae  a ship's 
weight , 

3.  Substitution  of  an  incorrect  letter  or  digit, 

4.  Duplication  of  a record  on  the  same  tape,  usually  when 
an  initially  mis-punched  record  was  corrected  and  the  ori- 
ginal not  removed. 
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Such  errors  are  particularly  serious  when  they  affect 
the  ship's  identifying  information,  as  these  fields  are  used  to 
match-up  ships  from  successive  weeks.  To  minimize  the  extent  of 
this  problem  ships  were  matched  via  a statistical  correlation 
on  the  identifying  data. 

In  principle,  each  tape  should  have  contained  data  on  the 
same  ships,  in  the  same  order.  In  practice,  however,  several 
types  of  violations  of  this  rule  could  (and  did)  occur: 

1.  Misplaced  records,  out  of  the  regular  order 

2.  Omitted  and  inserted  records,  when  no  data  was  reported 

for  a particular  ship  during  a particular  week. 

The  final  run  statistics  from  program  GAMMA  gives  a count 
of  records  rejected  due  to  excessive  keypunching  errors  and  of  the 
number  of  new  records  which  could  not  be  matched.  The  results  for 
week  47  can  be  taken  as  typical.  Here,  out  of  a total  of  13,130 
input  records,  1 was  mispunched  beyond  recovery,  20  were  duplicates, 
815  were  out  of  order,  and  4,644  did  not  match  any  entries  from 
the  previously  processed  weeks.  Hence,  only  7,650  records,  or 
slightly  over  half,  did  not  suffer  from  a serious  problem  in 
interpretation  or  matching,  even  before  problems  of  incomplete 
data  are  considered. 

The  data  for  certain  weeks  were  however,  much  poorer,  due 
to  hardware  failure.  In  particular,  the  tapes  for  weeks  37.  36, 

3S,  31,  30  and  29  contained  far  fewer  readable  records  than  ex- 
pected, less  than  10,000  each,  and  weeks  28  and  27  were  completely 
unprocessible.  This  reflects  the  deterioration  of  the  tapes  with 
age,  with  non-readable  records  being  dropped  on  input  by  the  pro- 
cessing system.  In  addition,  difficulties  are  known  to  have 
occurred  when  the  copies  of  these  tapes  used  in  this  study  were 
initially  prepared  by  DOT.  Fortunately,  the  other  tapes  were  in 
much  better  shape. 


Since  it  was  beyond  the  limitations  of  this  study  to 
replace  the  missing  data  by  re-punching  the  70 , 000  or  so  missing 
records,  little  could  be  done  to  correct  these  technical  diffi- 
culties, 

3.  Auxiliary  Tables. 

The  major  auxiliary  table  associated  with  the  DOT  input 
tapes  is  The  Port  Table,  containing  a list  of  all  ports  referred  to 
in  the  reports  contained  on  these  tapes.  Each  port  was  identified 
by  an  abbreviation,  which  associated  it  with  a specific  entry  in 
The  Port  Table.  The  Port  Table  contained  the  latitude  and  longi- 
tude of  the  port,  along  with  an  identifying  number. 

The  original  DOT  Port  Table  contained  a few  errors  or  in- 
consistencies which  were  corrected  for  use  in  this  study.  Some 
latitutde  and  longitude  information  were  Incorrect.  Many  ports 
had  several  entries  in  the  port  table  corresponding  to  differ- 
ently spelled  abbreviations  used  in  the  Input  Tapes.  These  were 
corrected  by  assigning  the  same  port  number  to  all  such  variants, 
since  they  could  not  be  deleted. 

4.  Conclusions. 

The  DOT  Marine  Traffic  Generator  Input  Tapes,  compiled  fro© 
Lloyd's  Shipping  Register,  suffers  from  several  weaknesses  for 
the  purpose  of  the  present  study.  It  is  incomplete  in  regards 
to  the  number  of  ships  monitored  and  the  movement  data  for  each 
ship.  It  is  systemtically  biased  in  this  regard  against  short 
voyages  and  some  classes  of  voyages  between  Communist  ports. 
Technically,  many  of  the  tapes  are  in  poor  condition,  two  having 
been  essentially  unusable.  Nevertheless,  it  was  the  best  data 
available  at  the  time  of  this  study,  and  every  effort  was  made  to 
overcome  its  limitations. 
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Maximum  information  was  extracted  from  the  tapes  by  comparing 
data  from  successive  weexs.  The  itinerary  of  each  ship  was 
reconstructed  as  well  as  practical,  and  the  number  of  days  spent 
by  ships  travelling  letween  each  pair  of  ports  accumulated.  This 
information  was  then  used  to  produce  the  shipping  distributions. 
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CHAPTER  III 

Trend  and  Growth  Analysis 

1.  Discussion  of  Available  Sources 

There  were  a number  of  very  informative  sources  available 
to  use  on  the  World's  Tanker  Fleet.  The  "Analysis  of  World 
Tank  Ship  Fleet" ^ for  the  years  1974,  1975  provided  data  on 
the  physical  characteristics  of  the  fleet;  size,  weight,  capacity, 
speed,  age,  country  of  registry,  etc.  as  well  as  very  useful 
information  on  the  economic  life  of  tankers  and  the  number  of 
new  construction  orders  placed  in  each  of  the  years  for  1964  to 
1975.  Some  analyses  of  special  vessels  (i.e.,  0/B/0' s, 0/0' s) were 
also  given. 

The  "Analysis  of  World  Tank  Ship  Fleet"  is  a Sun  Oil 
Company  publication.  The  Sun  Oil  Company  published  its  first 
fleet  analysis  in  1947  and  has  continued  to  issue  its  analysis 
yearly  since  1953.  Each  year  the  report  has  included  the  results 
of  a research  effort  on  tank  ship  economics.  The  research  for 
1974  was  in  the  area  of  total  carrying  capacity  of  the  fleet  and 
the  impact  of  fuel  costs  on  the  industry.  The  1975  report  contains 
a special  section  describing  a method  for  using  available  informa- 
tion to  predict  tanker  rates  in  the  spot  market. 

18 

"Shipping  in  the  Seventies"  is  a background  paper  pre- 
pared by  the  United  Nations  in  order  to  assist  consideration  of 
shipping  problems  in  general  and  the  position  of  the  developing 
nations  in  relation  to  them.  Part  One  addresses  the  perspectives 
and  problems  in  world  shipping  concentrating  on  technological 
progress,  freight  rates,  optimum  ships  and  scheduling.  Part  Two 
discusses  implementation  of  the  international  development  strategy 
in  shipping  and  ports,  with  specific  measures  for  implementation 
outlined.  This  brochure  provided  some  insights  into  the  future 
distribution  of  ship  ownership. 
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"An  Investigation:  U.S.  Import  Dependence  for  Mineral 

7 

Resources?  Super-Bulk  Carriers,  and  Deepwater  Ports  Development" 
is  a Master's  Thesis  written  by  Lt.  Commander  Corydon  Rouse 
Gifford  as  part  of  the  requirements  for  his  degree  at  the  Naval 
Postgraduate  School.  His  thesis  examines  " (1)  the  increasing 
import  dependence  of  the  United  States  on  foreign  sources  of 
major  bulk  fuel  and  non-fuel  mineral  resources  which  appear 
to  have  potential  requirements  for  deep-water  ports  or  terminals, 
(2)  principal  bulk  commodity  ports  and  commodity  movements,  (3) 
trends  in  the  construction  of  ocean  bulk  carriers,  (4)  the 
economics  of  "super"  bulk  carriers,  (5)  the  consequences  of  the 
United  States  not  providing  adequate  facilities  to  accomodate 
"super"  bulk  carriers  and  (6)  recent  events  pertaining  to  deep- 
water port  end  terminal  developments  by  the  United  States." 

We  v/ere  quite  interested  in  his  assessments  of  the  consequences 
of  various  courses  of  action  the  U.S.  might  be  expected  to  pursue; 
for  clearly  the  future  distribution  of  tankers  throughout  the 
world,  in  fact  the  entire  tanker  industry  will  be  ultimately 
related  to  the  future  course  taken  by  the  U.S. 

28 

“A  Forecast  of  1970-1985  World  Shipping"  ' by  Robert  V. 
Wiederkehr,  Sept. , 1971  is  a report  prepared  by  SACLANT  ASW 
Research  Centre  for  NATO.  This  report  is  part  of  a study  on 
the  future  use  of  shipping  in  the  North  Atlantic  by  NATO  nations. 
It  has  the  forecast  numbers,  average  speeds,  average  tonnages, 
fleet  capacities,  annual  trade,  and  fleet  productivity  of 
merchant  shipping  up  to  1985. 

"Trends  in  World  Shipping  with  Particular  Reference  to 

3 

Large  Crude-Oil  Carriers"  by  Peter  Burton,  July  1974  is  a 
second  SCALANT  ASW  Centre  manuscript.  This  document  reviews  and 
updates  the  previous  SACLANT  document.  The  forecast  figures 
by  Robert  V.  Wiederkehr  were  found  to  be  low  for  tankers  and  bulk 
carriers  and  high  for  container  ships.  This  paper  studies  the 
1974  Ptatus  and  short  term  trends  of  World  Shipping  as  a first 
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step  in  estimating  total  numbers  and  types  of  ships,  and  their 
concentration  routes  and  in  areas  of  importance  to  NATO  in  the 
eighties . 

19 

"Shipping  Statistics  and  Economics"  published  monthly 
by  H.P.  Drewry  (Shipping  Consultants)  Limited  is  a journal 
designed  to  provide  shipping  data  and  marketing  trends  to  shipping 
companies.  The  November  1975  and  March  1977  issues  were  used  to 
provide  us  with  insights  into  the  changes  in  shipping  trends  from 
1975  to  1977  and  an  overview  of  the  market  trends  in  those  years. 
These  journals  included  monthly  tables  on  numbers  of  vessels  by 
ship  type.  We  used  these  numbers  as  part  of  our  comparison 
analysis  to  arrive  at  realistic  total  ship  numbers.  The  growth 
factors  reflected  there  were  instrumental  in  arriving  at  our  own 
growth  factors. 

13 

"Recent  and  Projected  Trends  in  the  Shipping  Industry" 
a technical  paper  produced  July  28,  1976  by  Tetra  Tech,  Inc. 
analyzes  the  effects  upon  shipping  brought  about  by  the  various 
crises  of  the  mid-70's.  This  paper  updates  a 1975,  Tetra  Tech 
report. 
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"A  Statistical  Analysis  of  the  World's  Merchant  Fleets" 
by  the  U.S.  Department  of  Commerce  is  a world  inventory  of  ocean- 
going vessels  of  1,000  gross  tons  or  more.  The  world  totals 
for  the  various  vessel  types  for  the  years  1972,  1974,  and  1975 
were  taken  from  this  source  as  it  provided  us  with  a consistant 
ship  profile  and  data  collection  process. 

The  "Pacific  Shipper""^  is  a weekly  journal  devoted  to 
shipping  schedules  from  and  to  the  Pacific  ports  of  the  United 
States.  It  is  a reliable  gauge  of  the  volume  of  shipping  to 
major  ports  around  the  world  from  the  U.S.  Pacific  ports. 

In  order  to  determine  the  volume  of  shipping  to  major 
ports  around  the  world  from  the  U.S.  Pacific  ports,  an  analysis 
of  the  "pink  pages"  of  the  Pacific  Shippers  was  made.  The 
following  procedure  was  used: 
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The  Pacific  Shippers  (dated  weekly  from  July 
of  1973  through  June  of  1977)  were  used.  Only 
the  copy  for  the  last  week  of  each  month  was 
examined. 


2.  Traffic  to  the  following  ports  from  the  U.S. 
Pacific  Coast  was  examined:  Kobe,  Manila, 
London,  Rio,  Capetown,  Naples,  Piraievs, 
Bombay,  Honolulu,  Sydney,  Balboa,  Singapore, 
LeHarve , Rotte  rdam . 

3.  Since  the  traffic  leaving  the  West  Coast  ports 
(San  Francisco,  Los  Angeles,  Portland,  Seattle 
and  Vancouver)  is  essentially  identical,  only 
traffic  leaving  the  San  Francisco  - Oakland 
area  was  studied. 


4.  When  the  copy  for  the  last  week  of  a given  month 
was  analyzed,  for  example,  the  issue  dated 
June  27,  1977,  only  ships  leaving  through  the 
following  month  are  recorded.  For  June  27, 

1977,  for  traffic  to  London  (page  130) , there 
are  37  ships  leaving  San  Francisco  for  London 
from  July  1 through  July  29.  The  7 ships  leaving 
between  June  27-30  are  not  recorded  for  this 
issue  (they  were  recorded  for  the  May  30,  1977  issue). 
The  2 ships  mentioned  for  August  2 and  4 are  not 
recorded  for  this  issue  (they  will  be  recorded  for 
the  July  25,  1977  issue) . 


Figures  3-1  to  3-14  are  a collation  in  histogram  form  of 
the  data  extracted  from  the  Pacific  Shipper.  A perusal  of  the 
histograms  indicates: 

a.  Steady  volume  of  traffic  from  '73-' 77  to  Kobe, 

Rio,  Capetown,  Honolulu,  Sydney. 
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b.  A steady  increase  in  traffic  from  ' 73-,77  tf> 

Manila  and  Rotterdam. 

c.  Erratic  traffic  volume  to  Bombay. 

d.  London/  Naples  and  Piraievo  all  experienced 
an  increase  in  the  volume  of  traffic  arriving 

at  their  ports  in  August  of  1975.  This  is  probably 
due  to  an  increase  in  the  reporting  net  of  the 
Pacific  Shipper.  Balboa  shows  a similar  increase 
at  the  same  time  although  it  was  less  dramatic. 

e.  In  October  1975,  Singapore  experienced  the 
beginning  of  an  increase  in  incoming  traffic  that 
has  continued  steadily.  LeHarve  showed  the 

same  phenomenon  in  September  1975. 
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figure  3.2.  Histogram  of  Shipping  Traffic  from  U.S.  Pacific  Coast  to  Manila 
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Figure  3.3.  Histogram  of  Shipping  Traffic  Erom  U.S.  Pacific  Coast  to  London 
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Figure  3.9,  Histogram  of  Shipping  Traffic  from  U.S.  Pacific  Coast  to  Honolulu 


Figure  3.12.  Histogram  of  Shipping  Traffic  from  u.S.  Pacific  Coast  to  Singapore 


SHIPS  TO  LE  HARVE 


Figure  3.13.  Histogram  o±  Shipping  Traffic  from  U.S.  Pacific  Coast  to  Le  Harve 


Figure  3.14.  Histogram  of  Shipping  Traffic  from  O.S.  Pacific  Coast  to  Rotterdam 


2.  Methods  and  Sources  Finally  Used 

In  order  to  project  our  1972  data  from  the  Lloyd's  tapes 
to  the  year  1978  we  need  to  estimate  growth  factors  for  the 
various  ship  categories  in  the  major  ocean  areas  of  the  world. 

The  objective  is  to  compare  the  total  number  o?  ships  of  a given 
class  in  a particular  ocean  area  in  the  year  1972  with  the  same 
measurement  in  the  year  1978.  The  ratio  of  these  two  numbers 
then  being  the  growth  factor. 

The  difficulty  with  this  procedure  is  that  no  single 
definitive  source  exists  for  these  measurements.  Some  sources 
give  measurements  for  years  other  than  3.972  and  1978.  Since 
this  report  is  being  generated  before  1978  statistics  are 
available  on  any  measurements  for  that  year  are  projected 
measurements.  Each  source  has  its  own  procedure  for  collecting 
data  and  its  own  breakdown  of  vessel  type,  making  it  far  too 
inaccurate  to  use  data  from  one  source  for  1972  and  estimates 
for  1978  from  a second  source. 

We  are  most  interested  in  merchant  vessels  and  three  classes 
of  tankers.  Tanker  I’s  are  small  tankers  with  1,000  to  70,000  dwt. 
Tanker  II*  s are  large  tankers  and  combined  carriers  in  oil  and 
are  between  70,000  and  300,000  dwt,  Tanker  Ill’s  are  super 
tankers,  that  is  tankers  of  over  300,000  dwt. 

Since  it  is  sufficient  to  consider  only  a world  wide  growth 
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figure  for  merchant  vessels  (see  PSI  report,  "Ship  Density  ^ 

Estimates  for  Traffic  Along  Selected  Routes  in  the  Year  1990"  ) 

we  estimated  the  grow*  h factor  by  using  the  number  from  "Statistical 

24 

Analysis  of  World  Merchant  Fleet"  years  .1972,  1974  and  1975  and 
extended  them  linearly  to  1978.  The  growth  ratio  being  the 
simple  ratio  of  the  1972  figures  and  the  projected  1978  numbers. 

The  tanker  case  is  more  complicated.  The  growth  factors 
for  these  vessels  are  not  uniform  over  the  world.  But  whi.lo 

f\  C 

several  of  the  sources  ’give  world  wide  figures  no  source  gives 
comprehensive  data  for  specific  ocean  areas.  Certain  areas 
reflect  the  same  growth  pattern  as  the  world  in  general,  such 
as  the  North  Atlantic  (see  PSI  report  referred  to  previously) 
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Tankers 


World  Totals 


1.000  - 175,000  dwt 

175.000  - 300,000  plus  o/o 


1972 

2,931 

558 


Projected 

1978 

3,000 

894 


Growth 

Factor 

1.02 

1.60 


Specific  World  Areas 

Honolulu 

South  Pacific 

U.S.  - Japan 

Panama  - West  Coast 

North  Atlantic  & Remainder 
of  World 


Growth  Factor 
Tanker  I Tanker  II 


0.94 

1.21 

1.08 

1.15 

1.02 


Merchants 


1972 


1974 


Projected 

1978 


12,889 


13,253  13,600 


36 


ttm 


’■  ■’  I 


0.94 

1.60 

1.08 

1.30 

1.60 


Growth 
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1.06 


but  other  areas  do  not  follow  the  overall  pattern  of  the  world. 
To  overcome  this  limitation  two  types  of  growth  factors  were 
computed.  The  first  is  growth  factors  for  world  wide  shipping 
in  specific  tanker  categories.  The  second  is  growth  factors  for 
total  shipping  in  specific  ocean  areas  for  which  we  could 
obtain  reliable  data.  Specifically  we  used  the  Pacific  Shipper 
to  calculate  growth  factors  for  traffic  to  U.S.  western  ports 
from  Honolulu,  South  Pacific,  Japan  and  Panama.  We  collated 
the  available  data  for  the  years  mid  1973  to  mid  1976  and  used 
curve  fitting  techniques  to  eliminate  seasonal  fluctuations 
and  obtained  the  yearly  trends.  Next  we  extrapolate  backward 
to  1972  as  our  growth  factor  for  these  ocean  areas. 

However  since  super  tankers  were  a rarity  in  1972  but  have 
undergone  a dramatic,  but  erratic  growth  pattern  in  the  ensuing 
years  we  do  not  attempt  to  compute  a growth  factor  for  these 
vessels. 

In  computing  the  world  wide  growth  factors  for  the  two 

smaller  categories  of  tankers  we  were  fortunate  to  find  that 

29 

the  "World  Shipping  Statistics  1976"  carries  statistics  for 
the  year  1972  and  projections  for  1978.  The  projections  for 
1978  are  based  upon  their  figures  for  1976  plus  orders  for 
ships  scheduled  to  be  delivered  by  1978.  The  growth  factor 
is  a simple  ratio  between  the  two  years. 
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CHAPTER  IV 

Distribution  and  Methodology 


1.  ORE  Outputs  and  Cubic  Splices. 

This  chapter  deals  with  the  problem  of  distributing  the 
densities  across  the  routes  by  one  degree  squares.  Various  input 
tables  were  needed  to  accomplish  this  task.  The  use  of  the  Port 
Table,  Route  Trees,  Routing  Table,  Sextiie  Points  and  various 
mathematical  methods  for  distributing  the  densities  were  obtained 
from  the  ORE  report.  Because  of  the  multitude  of  ports  and  the 
number  of  routes  connecting  them,  the  use  of  the  above  tables 
provided  us  with  a systematic  way  of  handling  the  data. 

The  Port  Table  is  a listing  of  2,425  ports  that  were  taken 
from  the  Lloyd's  data  for  1972  by  the  Department  of  Transportation. 
Due  to  the  number  of  ports  in  existance  a simplification  was 
warranted  in  order  to  provide  a more  manageable  data  base,  ihere- 
fore,  a small  class  of  major  ports  were  determined  from  the  2,425 
ports.  These  69  ports  are  known  as  canonical  ports.  Also,  each 
non-canonical  port  will  be  associated  to  one  and  only  one  can- 
onical port.  This  is  determined  by  searching  for  the  closest 
canonical  port  in  the  area. 

The  information  on  the  Port  Table  is  of  the  form: 


684  FREEPB  FREEPORT,  BAHAMAS  26N  78W  17  22 


684  is  the  port  number;  FREEPB  is  the  abbreviation  of  the  port 
name  and  FREEPORT,  BAHAMAS  is  the  name  of  the  port.  On  the  right 
hand  side  we  have  the  latitude  and  longitude  of  the  ports  loca- 
tion. The  next  integer,  17,  corresponds  to  the  canonical  port 
to  which  Freeport,  Bahamas  is  associated  with.  22  represents 
the  distance  between  the  port  and  its  canonical  port.  The  units 
are  that  of  a convenient  metric. 
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Just  as  the  number  of  ports  were  substantially  reduced  by 
considering  a small  set  of  canonical  ports#  so  may  the  number  of 
routes  be  reduced  by  considering  a set  of  85  canonical  routes. 

The  Routing  Table  is  a 69  x 69  matrix  whose  (a,  8)  entry  symbolically 
represents  the  route  from  port  a to  port  8.  The  entries  contain  two 
pieces  of  information.  The  first  is  an  alphanumeric  followed  by  an 
integer#  say  ra,  which  lies  between  1 and  85.  The  possible  entries 
in  the  table  are  as  follows s 

i)  N=0  represents  no  route  between  the  ports.  This  situa- 
tion arises  when  a = 6/  that  is  port  a equals  port  6. 

ii)  ?=m  represents  the  case  where  port  a and  port  8 are  not 
connected  by  a canonical  route,  that  is,  there  exists  at  least  one 
intermediate  port  between  a and  8.  The  port  m is  the  first  inter- 
mediate port  encountered.  As  one  would  expect#  the  route  between 
a and  ra  is  a canonical  route.  Thus,  the  first  link  of  the  overall 
route  is  from  port  a to  port  m.  Now,  one  must  examine  the  route 
from  port  m to  port  8.  If  there  is  a canonical  route  then  one  is 
done,  if  not,  continue  recursively  until  the  destination  is  reached. 

iii)  D=m  designates  that  multiple  routing  exists.  This  case 
occurs  when  the  route  taken  is  dependent  upon  the  ship  type.  For 
example,  the  shortest  route  between  the  Persian  Gulf  and  Northern 
Europe  is  via  the  Suez  Canal,  yet  large  tankers  and  supertankers 
must  go  around  the  Cape  of  Good  Hope.  In  this  case  the  integer  m 
designates  the  mth  line  in  the  Multiple  Routing  Table.  This  line 
will  contain  two  entries,  the  first  entry  is  the  first  intermediate 
port  for  the  non-supertankers,  and  the  second  is  for  supertankers. 
Then  one  would  proceed  as  in  case  ii. 

iv)  If  a blank  exists  for  the  alphanumeric  field,  then  the 
following  number  indicates  which  direct  route  is  to  be  taken. 
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The  set  of  85  routes  can  be  found  by  looking  at  the  Route 
Trees  Table.  In  many  routes , seasonal  variations  affect  the  dis- 
tribution of  traffic  across  the  route#  on  the  other  hand  entire 
route  envelopes  are  shifted  by  the  changing  of  seasons.  Therefore# 
we  have  four  sets  of  Route  Trees#  a separate  one  for  each  season. 
One  must  note  that  the  structure  of  the  route  trees  does  not  change 
with  seasons;  however,  the  parameters  contained  within  the  tree 
structure,  such  as  the  latitude  and  longitudes  do. 

In  some  cases,  the  traffic  travelling  in  one  direction  uses 
a route  quite  different  from  the  traffic  moving  in  the  opposite 
direction.  One  example  of  this  is  Route  23  (Punta  Arenas  - 
Torres  Strait) . The  entrance  route  to  Torres  Strait  is  much 
farther  north  than  the  egress  route.  When  this  is  the  case#  the 
tree  is  structured  so  that  there  are  two  routes  between  the  ports. 
Since  a tree  is  a graph  without  loops,  this  implies  that  there  has 
to  be  two  leaves  for  one  of  the  ports.  One  leaf  is  the  vertex  used 
for  egress  routing,  the  other  leaf  is  used  for  entrance  routing. 

The  distinction  between  entrance  and  egress  routes  is  preserved 
in  the  tree  structure  by  the  ordering  of  the  vertex  indices.  A 
binary  branch  is  made  with  a terminal  zero  being  added  to  the  end 
of  the  egress  branch  point  identifier#  and  a terminal  one  to  the 
entrance  branch  point  identifier.  For  example#  Route  23  divides 
at  Ident.  10  (Strait  of  Magellan).  Ident.  100  is  the  egress  branch 
point  while  Ident.  101  is  the  entrance  branch  point.  All  identi- 
fiers are  written  in  base  2.  Figure  (4.1)  shows  the  tree  structure 
with  identifiers  for  Route  23, 

The  identifiers  describe  the  location  of  the  turn  point  in 
the  route  tree  structure.  These  identifiers  are  so  linked  that 
it  is  possible  to  progress  either  forward  or  backwards  along  the 
route  in  the  tree  structure.  In  order  to  move  forward  there  axe 
three  possible  circumstances:  1)  You  are  at  a turn  point#  in 

which  case  a zero  is  added  after  the  last  digit  of  your  current 
vertex  identifier  in  order  to  obtain  the  identifier  for  the  next 
vertex  along  the  route,  2)  You  are  at  a branch  point#  in  which 
case  a zero  is  added  to  continue  down  one  branch  of  the  route#  or 
a one  is  added  after  the  last  digit  of  your  current  vertex  iden- 
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Tree  Structures 
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tifier  to  continue  down  another  branch,  3)  You  are  at  a route  end 
point  in  which  case  you  must  reverse  direction  in  order  to  stove 
along  the  route. 

To  progress  backwards  along  the  route  is  much  simpler  than 
a forward  progression.  To  retrace  the  route  simply  drop  off 
the  last  digit  of  the  current  vertex  identifier  in  order  to 
obtain  the  identifier  of  the  previous  tree  vertex. 

The  second  and  third  columns  of  the  Route  Trees  contain  the 
latitude  and  longitude  of  the  route  point  represented  by  the  tree 
vertex. 

The  fourth  column  of  the  Route  Trees  refers  to  the  vertex 
width.  Its  units  are  in  latitudinal  degrees,  and  it  is  a measure- 
ment of  the  routes  width  at  the  vertex. 

The  length  of  a route  segment  is  contained  in  column  five 
of  the  Route  Trees.  This  length  refers  to  the  previous  route 
segment  when  progressing  forward  through  the  route.  The  units  are 
in  latitudinal  degrees. 

The  sixth  column  of  the  Route  Trees  contains  the  edge  width 
which  is  again  a measurement  of  the  widest  point  on  the  previous 
segment.  The  units  are  also  in  latitudinal  degrees. 

The  seventh  and  last  column  refer  to  the  back  pointers. 

When  progressing  backwards  along  a route  the  back  pointers  in- 
dicate the  line  in  the  Route  Trees  to  which  the  present  vertex  is 
connected.  Whenever  one  begins  a new  route  the  back  pointer 
is  set  at  zero. 

An  important  factor  in  determining  the  route  density  is 
the  route  envelope.  It  is  an  area  of  the  ocean  between  two 
ports  where  95  percent  of  the  ships  will  travel.  The  width  of 
the  route  envelope  may  vary  from  route  to  route.  In  going  from 
one  port  to  another  not  all  ships  follow  the  best  path.  Differ- 
ent ship  tracks  will  approximate  the  optima:-  path  to  a greater 
or  lesser  extent.  The  longer  the  distance  aecween  ports. 
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the  greater  the  probable  deviation  from  the  optimum.  This 
is  especially  true  for  routes  that  cross  large  open  oceans  with 
no  land  navigation  pcint  available.  Areas  with  unstable  weather 
conditions  also  cause  ships  to  deviate  from  the  best  path.  Be- 
cause of  the  effects  of  these  variables,  a route  width  is  required 
in  order  to  describe  a route's  envelope. 

The  distribution  is  found  along  a vertical  cut.  A vertical 
cut  refers  to  a slice  of  the  route  taken  along  a fixed  longitude. 
The  reason  for  the  cut  is  to  take  a look  at  the  way  the  distri- 
bution function  varies  across  the  route . The  densities  along  the 
section  of  the  route  under  consideration  are  projected  parallel 
to  the  route  median  and  integrated  to  form  the  density  distri- 
bution. The  sextile  points  are  points  which  divide  the  distri- 
bution into  sixths.  For  a cut  along  a fixed  longitude  there 
are  latitudes  associated  with  these  sextile  points.  The  sex- 
tile  points  of  the  distribution  and  their  corresponding  latitudes 
are  found  and  printed.  These  sextile  point  latitudes  are  used 
by  a separate  program  to  reconstruct  the  distribution  using  a 
Type  I*  cubic  spline  to  connect,  the  sextile  points. 

The  Sextile  Array,  is  a table  containing  the  sextile  points 
by  route  for  a given  season. 

These  sextile  points  are  normalized  to  the  interval 
[-1,1],  and  the  five  intermediate  sextile  points  will  describe 
the  distribution  (up  to  a linear  transformation) . On  a route 
between  2 canonical  ports  the  envelope  is  described  by  either  a 
specif  .c  Type  i 1 cubic  spline  for  the  given  route  or  a general 
default  spline. 

In  order  to  produce  a histogram  which  represents  the  dis- 
I tribution  of  traffic  across  a route,  it  is  first  necessary  to 

| determine  how  the  density  varies  across  the  route.  The  following 

| method  was  used  to  accomplish  this  task. 
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Given  a fairly  well  defined  (e.g.,  Southern  Hemisphere) 
route/  and  assuming  the  same  distribution  function  is  used  along 
the  entire  route,  choose  a group  of  vertical  columns  of  cells 
(Fig  4-2)  which  cross  the  route.  These  columns  should  be  chosen 
in  an  area  where  there  is  a minimum  of  interference  from  other 
routes.  There  should  be  no  cross  routes  in  this  region,  nor 
should  this  area  contain  any  other  routes  in  close  proximity. 

If  only  one  vertical  column  is  used,  then  this  would  automati- 
cally be  a histogram,  for  the  density  histogram  which  we  are 
seeking  would  merely  be  a projection  of  this  column  (assuming 
that  the  width  of  the  route  doesn't  change).  (Note:  The  his- 

tograms are  somewhat  more  complex  than  shown,  for  to  be  exact 
the  projection  of  a square  of  area  onto  the  perpendicular  of  the 
route  is  not  the  standard  step,  but  rather  a trapezoid  (Fig  4-3) . 
This  is  probably  sufficiently  minor  to  be  ignored,  but  is  not 
computationally  difficult.).  As  long  as  the  route's  distribution 
may  be  assumed  locally  constant,  then  the  vertical  column  of  Figure 
4-4  is  a cut  of  a “cylinder"  type  of  distribution  which  is  a 
function  of  only  one  variable,  i.e.,  it  is  a slightly  stretched 
version  of  the  distribution,  the  elongation  factor  being  sec  0. 

In  the  case  of  more  than  one  column  (which  will  have  better 
statistical  confidence  for  distribution  prediction)  the  columns 
can  be  combined  appropriately  using  tan  0 to  the  histogram 


(Fig  4-5) . With  the  elongation  removed  we  now  have  the  empirical 


representation  of  y(  = ^ G(w)  sin  ttX  say)  so  6 and  G(w)  are 
"determined",  through  a smoothing  operation  discussed  below. 


The  histograms  are  discrete,  empirical  approximations  to  the 
true  distribution.  This  underlying  distribution  function  does 
not  have  the  step  function  for  the  histogram.  From  the  physics 
of  the  situation,  one  may  safely  assume  chat  there  exists  a 
density  function  for  the  distribution,  and  that  it  is  continuous. 
To  obtain  the  underlying  distribution  function,  it  is  necessary 
to  “fit"  a smooth  curve  (the  density  function)  h(x)  to  the  step 
function  (i.e.,  the  histogram).  The  density  function  will  be  the 
smallest  closed  set  such  that  the  density  function  is  always 
zero  outside  that  set. 
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Figure  4-5%  Histogram  of  Two  Columns  of  Data 

— and  — are  individual  columns, 

— is  the  sum  of  both. 


0 
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The  use  of  interpolating  functions  are  used  to  fit  a 
smoothed  distribution  function  H(x)  to  an  empirical  distribution 
H(x).  In  order  to  use  an  interpolation  function  one  only  needs 
a net  of  points  along  the  interval.  Since  the  histogram  itself 
is  really  a discrete  data  structure  this  approach  is  especially 
suitable. 


,N 


Let  M = {(x.,  H(x.)}”  be  a mesh  of  points  from  the  histo- 

1 1 U N 

gram's  distribution,  and  assume  that  the  net  {x^}q  covers 


t ”1 » 1] i i.e, 


x,  ,<x. 
l-l  l 


x0  = _1 


*N  = 1 


The  smoothed  density  function  should  be  continuous  (by  the  nature 


of  open  ocean  routes) , and  it  is  not  unreasonable  to  assume  it 
,1 


is  C (i.e.,  has  a continuous  first  derivative)  on  the  interval 
[-1,1].  Thus  the  distribution  function  will  be  C [-1,1],  (i.e., 
it  will  have  a continuous  second  derivative  on  the  interval 
[-1,1]). 

those  functions  which  satisfy  four  criteria: 


Let  the  class  of  possible  distributions  D3  consist  of 


i. 

It  is  C2[  -1,1] 

ii. 

It  is  a distribution  function  on 

[-1,1] 

iii. 

It's  extension, to  a distribution 
real  line  is  CA 

function  on  the 

iv. 

It  interpolates  the  mesh  M 
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For  an  open  ocean  route  envelope  one  would  not  expect  a rapidly 
fluctuating  density  across  a route  cut.  Based  on  this  empirical 
observation  and  other  inherent  properties  of  the  route,  the  cri- 
teria for  minimization  is  chosen  to  be  the  integrated  square 
derivative  of  the  density  function. 


/k(h' (x) }2dx 


This  analytical  expression  is  a measure  of  the  fluctuation  in 
the  density  function.  Thus  the  chosen  interpolating  function 


H(x)  is  that  function  from  the  set  which  minimizes  the  func- 


tional 


,b rd  Hi2, 

V dT}  dx 


Thus  out  of  all  functions  satisfying  conditions  i through  iv, 
that  one  which  minimizes  this  functional  will  be  the  'smoothed' 
distribution  used  to  model  the  distribution  of  traffic  across 
a route  cut.  Using  the  fundamental  identity  for  cubic  splines 
a generalization  of  Holladay's  Theorem  can  be  developed.  This 
generalization  implies  that  if  the  optimization  is  done  over  all 
functions  satisfying  i,  iii  and  iv,  then  the  minimizing  function 
is  a Type  I'  cubic  spline.  Furthermore,  the  minimizing  function 
is  unique.  {A  cubic  spline  is  of  Type  I'  if  the  first  derivative 
is  zero  at  the  endpoints  a and  b.}  If,  in  addition,  the  Type  I' 
cubic  spline  is  monotonic,  then  it  is  also  a distribution  func- 
tion, which  is  condition  ii.  In  other  words,  it  will  satisfy  the 


four  conditions  for  inclusion  in  D3  and  it  will  be  the  (unique) 


distribution  which  minimizes  the  above  functional.  The  monotonicity 
was  assumed,  however,  in  actual  data  reduction  the  cubic  spline 
turned  out  to  be  monotonic  in  every  case. 
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The  mesh  is  obtained  by  using  seven  sextile  points  uniformly 
located  along  the  distributed  mass  - i.e.,  the  seven  sextile 
points,  which  are  found  by  solving  the  following  equation  for 


!»  ■ 


H(xi)  = i/6 


i **  0,1,2, ...  ,6 


Now  the  equations  of  the  cubic  spline  may  be  developed  in  an 
algorithmically  implementable  form.  Since  the  interval  [a,b] 
was  defined  to  be  the  smallest  closed  interval  containing  the 
support  of  the  distribution,  it  is  immediately  seen  that 


xQ  = -1 
x6  = 1 


Let  A denote  the  forward  difference  operator 


Ax.  s x. - x. 
3 3+1  3 


» 

f: 


Define  the  fifteen  quantities 

Ax . 


^ j Ax^  + Ax.. 


3-1 


U.  * 


l - x. 

3 


Ax. 


■t  . OA  , - 

2 {A^  + + Axj-1} 

for  j - 1,2, 3,4,5 


ill  I 
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The  cubic  spline  is  given  in  terms  of  the  seven  coefficients 
(m^i  = 0,1,... ,6).  Since  we  are  dealing  with  a Type  I*  cubic 
spline  two  of  these  coefficients  are  fixed,  viz 


m6  = 0 

The  remaining  coefficients  are  given  by  the  linear  system 

Cll 


This  matrix  is  nonsingular,  and  the  system  may  easily  be  solved 
by  tridiagonal  matrix  algorithms.  The  density  function  fi(x) 
of  the  distribution  H is  then  given  by 

A (x.-x) (2x  ' . +x.-3x) 

h(x)  = m.  . — 3 - — , .2 

(Ax.^)2 

(x-x,  ,)(2x.+x.  ,-3x) 

- m -gl- 

j <A*j-l) 

(x.-x) (x-x.  .) 

+ -J — pL- 

(AXj_^) 

for  xe  (Xj^rXjl 
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This  equation  will  be  refered  to  as  equation  4.1. 


Thus  the  development  of  the  route  envelope's  distribution  is 
completed,  with  the  underlying  distribution  being  modelle4  by 
a Type  I*  cubic  spline. 

Now  we  will  address  ourselves  with  the  problem  of  mathe- 
matically modelling  the  route  envelope  outline  and  the  distri- 
bution of  traffic  along  a ’cut*  across  the  route.  Given  two 
ports  A and  B the  envelope  between  the  two  ports  will  be  deter- 
mined by  two  mappings.  One  is  from  the  distance  space,  A,  to  a 
•flat'  earth  (i.e.,  the  plane  R ) and  one  from  the  'flat*  to 
the  'round'  earth  (i.e.,  the  sphere  S) . Thus  the  total  mapping 
is  the  composition  of  the  two  mappings. 

The  distance  space,  A,  is 

A ® A x Q * {(X,u>);  0<A<1,  0<u£l} 


The  parameter  w corresponds  to  the  stochastic  aspect  of  the  en- 
velope, i.e.,  which  route  is  taken.  The  value  u>  = 0 corresponds 
to  the  route  along  one  edge  of  the  envelope,  and  u = 1 corresponds 
to  the  route  along  the  other  edge.  X denotes  the  normalised 
distance  along  the  route  (i.e.,  distance  divided  by  the  total 
length) . Thus  X = 0 corresponds  to  port  A and  X = 1 corresponds 
to  port  B. 

2 

The  mapping  from  A to  R will  be  considered  in  several 
cases,  beginning  with  the  most  elementary  and  then  generalizing. 


Case  a:  Port-to-Port  Segment. 

In  this  case  both  ends  of  the  cresent  are  points  and 
there  are  no  protrusions  of  land  to  cause  land  avoidance.  On 
the  flat  earth  the  route  envelope  is  bounded  by  sine  waves.  A 
rectangular  coordinate  system  is  established  such  that  port  A 
is  at  the  origin  and  port  B lies  on  the  positive  x axis,  bet 
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Figure  4 . 6 Envelope  for  Case 
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26  be  the  maximum  width  of  the  envelope.  Let  g(yjxQ)  be  the 
density  function  of  the  distribution  of  traffic  across  the  route 


at  a fixed  distance  xQ  along  the  route  (0<Xg<L) 


Note  that 


g(.jxQ)  is  a density  function  concentrated  on  the  interval 
[-v,v]  where 

^0 

v = 6 sin  L 

Let  H<z;Ag)  be  the  distribution  function,  scaled  between 
[-1,1]  so  that 


+6z  sin  xA, 


H(  z;  Aq) 


= / g(y; 

J t*  _ f 


LXQ)dy 


-1<z<1 


-6  sin  xA. 


and  let  G{a);AQ)  be  the  inverse  function  of  H(2;Ag).  Thus » 
G(.;Aq)  has  domain  [0,1]  and  range  [-1,1).  Therefore,  the 
mapping  for  this  case  from  (A,w)  to  (x,y)  is 


X = LA 


y = 6G(w;A)sin  xA 


In 


Case  b:  Sloped  End  Condition. 

In  this  case,  the  lefthand  point  is  still  port  A but 
the  route  is  merely  a segment  of  the  total  route  (Fig  4-8) # 
thus  it  ends  at  a turnpoint  (e.g.,  strait,  cape)  in  the  ocean. 
As  before  26  is  the  maximum  envelope  width,  and  now  the  end 
of  the  segment  has  a width  of  2c  rather  than  zero.  Let  L be 
the  distance  to  the  midpoint  of  this  end  segment.  Note  that  if 
the  angle  of  the  end  segment  to  the  axis  (9)  is  not  x/2,  then 
the  'cuts'  across  the  route  envelope  corresponding  to  a fixed 
distance  A are  no  longer  perpendicular  to  the  axis.  For  each 


fixed  X e (0,1)  let  G(.jXq}  be  the  inverse  of  the  scaled  distri- 
bution across  the  routes  at  the  cut  X = XQ.  As  before,  G(.,*X0) 
maps  t0, 11  onto  [-1,11.  The  angle  0 is  measured  counterclockwise 
from  the  exterior  axis  and  is  between  0 and  s.  Thus  the  mapping 
for  case  b from  (X,w)  to  lx,y)  is 

x = (L+eG (ej X)  cos  8)X 
y » eG(w,X}X  sin  9 + $G(w;X)  sin  ttX 

where  the  constant  is  determined  from  the  following  equation 


where  y = (6-Q/2) 


and  c = e sin  9 


Case  c:  General  End  Conditions. 

If  both  ends  have  nonzero  width,  the  situation  corresponds 
to  Figure  4-9.  Both  angles  0 and  9 are  measured  counterclockwise 
from  the  exterior  axis  and  are  between  0 and  n.  The  distance 
L is  measured  from  the  midpoint  of  one  end  to  the  midpoint  of  the 
other.  The  midpoint  of  one  is  taken  as  the  origin,  and  the  mid- 
point of  the  other  as  the  positive  x axis. 

Let  c « e sin  9 

l » e sin  9 


Then  with  previous  definitions,  the  mapping  is 


e cG(o;  X)cos  0+(L+G{i*;  XHe  cos  9-e 
«=  GCu}X)U  ♦ (C-O-l  * $ sin 


COS  9 11 x 


56 


■iMfiaafi 


y 


where  <{>  is  a constant.  Proceeding  as  in  Case  b,  we  find  that 
<j>  may  be  computed  as 


<f>  ~ Y 


where 


Y = 6 - U+U/2 

2 

Thus  the  general  mapping  from  A to  R has  been  developed.  Next 

2 

is  the  mapping  from  the  image  of  A (in  R ) into  S. 

2 

The  plane  R had  units  of  length  attached  to  the  parameters 
L,  6,  e and  e,  but  all  equations  were  developed  without  recourse 
to  specific  dimensional  units.  References  to  angles  will  be  made 
in  terms  of  radian  measure.  The  spherical  earth  has  a well  es- 
tablished coordinate  system  using  latitude  and  longitude,  which 
we  will  be  using.  Adopting  the  conventions  that  Southern  la- 
titudes and  Eastern  longitudes  are  negative,  and  that  latitudes 
and  longitudes  is  expressed  in  decimal  degrees,  we  have  the  follow 
ing  coordinate  representation  of  the  earth  (minus  the  two  poles) ; 

S = { (a,0) :-90<S<90,-180<a<180} 

Let  (a^,Bj)  and  (c t^ie  two  endpoints  (longitude -latitude 
pairs)  of  the  route  segment,  then  the  route  in  R2  will  have  to 
be  mapped  such  that  its  endpoints  (or  midpoints  of  end  segments) 
map  into  these  two  points.  The  procedure  used  to  do  this  is  as 
follows. 
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Assume  that  the  points  are  not  diametrically  opposed  on 
the  sphere,  so  that  there  is  a unique  geodesic  between  them. 

If  the  route  is  not  vertical,  (which  would  send  the  great  circle 
through  the  poles)  then  the  geodesic  route  is  a segment  of  a 
great  circle,  which  has  an  equation  of  the  form 


tan  0 = y sin  a + 6 cos  a 


" V+s2  sin  (CMV 


where  the  positive  square  route  is  taken.  Substituting  the 
two  endpoints  in  'this  equation  gives 


Y = {cos  a2  tan  B1  - cos  tan  02>/T 
6 = {sin  tan  B?  - sin  a 2 tan  3-^J/T 

where  T = sin  (o^  - a2)  ^ 0 

To  parametrically  map  the  segment  (0,L)  of  the  x axis  in  the 
2 

R space  into  the  geodesic  between  the  two  points  (ct^,  3^)  and 
(a 2»  B2)  consider  three  cases. 

Case  a.  0^,  &2  not  both  zero: 

This  is  the  general  case,  where  at  least  one  endpoint  does 
not  lie  on  the  equator.  The  point  where  the  great  circle  up- 
crosses  the  equator  (i.e.,  > 0 at  the  crossing,  see  Figure 
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4.10)  is  (cXq,  0),  and  from  the  great  circle  equations  we  have 

tan  aQ  = -6/y 

The  correct  inverse  function  is  given  by  the  argument  function 

ct0  = Arg(y  - i6) 

Lot  the  arc  length  s along  the  great  circle  be  measured  from  the 
uperossing  point  (aQ,  0),  increasing  the  westward  direction.  Let 
it  be  measured  in  units  of  latitudinal  degrees  (1  degree  = 60  nm), 
so  that  the  arc  length  of  the  entire  circle  is  360°.  We  wish  to 
parameterize  the  route  segment  of  the  great  circle  using  oe[0,  1] 
such  that  o=0  corresponds  to  the  point  (a2>  B2>  and 


Lot  si  be  the  arc  length  (from  the  upcrossing)  to  (cu,  0^), 
then  by  the  law  of  cosines 

cos  s^  = cos  cos  ( a ^ - aQ) 

cos  s2  = cos  02  cos  (a2  “ aQ) 

To  uniquely  specify  the  value  of  si  e(-180°,  180°),  use  the 
relationship 

s^  = Arccos  ^cos  0^^  cos  (ot^  - aQ)j 
s1  = sgn  (31)s1 

Specification  of  the  correct  value  of  s2  is  somewhat  more  com- 
plicated since  the  route  must  be  traversed  the  'right  way' 
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Figure  4-10  Definition  of  Angle  a0  and  Upcrossing  Point 


around  the  earth,  and  there  cannot  be  any  discontinuities  of 
2ir  (or  rather  360°)  in  the  arc  length  function: 

s2  = Arccos  ^cos  8g  cos  (ctg  - otQ)^ 

s2  **  sgn  (Bg)  Sg  + k. 360° 

where  ke{  - 1,  0,  1}  is  chosen  such  that 

|s2  - sx|  < 180° 

Given  a oe[0,  1],  let  the  arc  length  parameterization  of  a third 
point  be 


s3  = S1  + (S2  - V 0 

Then  the  corresponding  point  (ct„,  g„)  along  the  segment  of  the 

o d 

great  circle  satisfies 

tan  $3  sin  (otg  - a^)  + tan  8-^  sin  (a3  - ag)  + 

« 

tan  8g  s:i-n  - ctg)  = 0 

If  8^  is  not  zero,  then  sin  s^  is  nonzero.  Since  not  both 

8-^  and  8g  are  zero,  the  following  equation  can  be  used  to  find  8 


sin  Bo  = sin  s„  sin  8,  / sin  s, 
3 3 1 1 


= sin  s3  sin  Sg  / sin  Sg 
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-y.  r,/~-  ' n*'- 


The  inverse  sine  function  is  unique  since 


|03l  < 90° 

The  above  relation  follows  from  repeated  application  of  the 
law  of  sines.  Now,  from  the  law  of  cosines 

cos  (ag  - aQ)  = cos  sg  / cos 

thus  a3  is  found  using  the  relations 

u » Arccos  (cos  s3  / cos  gg) 

a3  = a0  + sgn  ^sin  S3^U 

(Note  that  unlike  the  previous  applications  of  the  law  of 
cosines,  to  resolve  the  ambiguities  it  is  necessary  to  use  the 
sign  of  the  sine}. 

Case  b.  Both  g's  zero. 

In  this  case  both  points  lie  on  the  equator.  Define 
ct2  = a2  + k.  360° 

where  the  integer  k is  chosen  such  that 

|ci2  - aj  < 180° 

Then  the  parameterization  of  the  route,  in  terms  of  oe{0,  1]  is 

a3  ° al  + <52  **  al)0 


0 


Case  c.  Equal  a*s. 

In  this  case  both  points  lie  on  the  same  meridian#  and 
the  parameterization  is 


e3  = + (b2  " 

Hence  in  each  case  the  parameterization  of  the  segment  has  been 
specified. 

Our  attention  is  now  turned  to  the  perpendicular  to  the 
mean  route.  By  this  we  mean  the  family  of  great  circles  that 
cut  the  given  route  at  right  angles.  Before  determining  the 
equations  for  the  general  case,  it  is  first  necessary  to  give  some 
definitions. 

In  the  planar  projection  of  the  earth,  the  slope  of  a 

great  circle  route  is  and  it  is  (naturally)  not  constant. 

This  slope  represents  the  ratio  of  an  infinitesimal  change  in 

vertical  (ordinate)  degrees  to  an  infinitesimal  change  in 

horizontal  (abscissa)  degrees.  On  the  spherical  earth,  however, 

degrees  are  not  a uniform  measure  of  distance.  The  slope  of 

dS 

the  great  circle  (on  S)  is  sec  0^-,  which  represents  the  ratio  of 
infinitesimal  change  in  vertical  (North-South)  nautical  miles  to 
an  infinitesimal  change  in  horizontal  (West-East)  nautical  miles. 
The  slope  of  the  perpendicular  is  the  negative  reciprocal  of  the 
slope  of  the  great  circle  route  (i.e.,  - cos  0/^-) , thus  changing 
back  to  a degree  coordinate  system  (i.e.,  the  planar  projection), 
the  slope  of  the  perpendicular,  m^  (as  a ratio  of  infinitesimal 
degrees)  is 
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»• *AA**>m*h' 


mi  ■ - oos2(i ' if 

= . cos2e  lt«l  / i-§as_£ 

2 2 

= - cos  3 sec  3 / (y  cos  a - <S  sin  a) 

= 1 / <i> 

where 

}p  - 6 sin  a - y cos  a 

Let  the  angle  the  perpendicular  makes  with  the  horizontal4 
(in  the  planar  projection)  at  the  intersection  with  the  great 
circle  route  be  x as  shown  in  Figure  4.11.  Then 


cos  X = 4>  / 

3in  X » 1 

where  positive  square  roots  are  taken.  The  great  circle  per- 
pendicular to  the  route  at  the  point  (ay&3)  is  given  by  the 
equation 


tan  = y sin  •>  <$  cos  a^ 
where  y and  6 satisfy  the  equations 

•V 

tan  &3  * y sin  a3  + 5 cos  a3 

2 

sec  tan  x 88  Y cos  a3  - 6 sin 
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X • <3  ■ $ 

~ cos  ou 

Y =>  sin  a,  tan  $,  + 

3 3 $3  cos  ^e3 

..  sin  a, 

6 = cos  a.  tan  0,  - * — 

*3  cos  e3 


where  ij)  is  \p  evaluated  at  the  point  (a^,  $3) . 


Since  the  width  of  the  route  envelope  will  be  small 
(compared  to  the  radius  of  the  earth)  it  may  not  be  necessary 
to  use  great  circles  for  the  perpendicular  displacements,  but 
rather  use  rhumb  lines.  A rhumb  line  approximation  to  this 
perpendicular  is  given  by  the  parametric  equations 

01  = 63  + t sin  x3 
al  = a3  + T COS  *3 

where  x is  the  parameter  and  x3  is  X evaluated  at  ( 

This  equation,  of  course,  is  much  simpler  than  the  great  circle 
perpendicular. 

Now  we  have  the  necessary  tools  to  begin  distributing  the 
densities.  By  use  of  the  Cubic  Spline  Table  each  route  will  have 
a designated  Cubic  Spline  associated  with  it,  where  the  number 
of  ship  days  will  be  distributed  across  a normal  (i.e.,  a 'cut') 
to  the  median  of  the  route  envelope.  Ultimately  the  density  will 
be  distributed  to  1°  square;..  <nen  doing  this,  the  cut  of  width 
h (where  h = edgelength  x 5)  will  intersect  the  one  degree  square 
boundaries  in  a manner  indicated  in  Figure  4-12. 


V 335* 
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route  envelope 


■^typical  cut 


Figure  4-12.  Components  of  a Typical 
Route  Segment 


It  can  be  observed  that  various  proportions  of  one  degree  squares 
are  located  within  the  cut.  Now  we  will  consider  the  various 


cases. 


Li 


Case  1:  Suppose  that  one  edge  of  the  1°  square  lies  in  the  cut# 

as  in  Figure  4-13-  Given  the  coordinate  of  the  vertices 
of  the  one  degree  square,  the 
equation  of  the  side  which  in- 
tersects the  cut  is. 


(P>  * v-.fe) 


y « ox  +'  0 


Typical 
( cut 


(0,B) 


><Vyr> 

'M®  square 


Figure  4-13.  Case  Is  One  edge 
of  the  1*  square  lies  in  the  cut. 


where 


ArV\ 

\vV 


Let  h(y)  denote  the  density  function  indicated  in  eq.  4,1 

for  the  interval  ty^^y^].  Upon  integrating  the  density  function 

in  noth  directions  one  obtains 


ctx+B 


/„ 


h(y)dy  = H(ax+6)  - H(0I 


= *0  + V + V2  * V3 


where  aQ,  a^,  a^,  are  constants  dependent  on  the  cubic  *•”.  Line, 

1°  sq.  vertices,  etc. 


-6/a 


t 2 3 

Volume  = I (a^+a^x+a^  +&3x  )dx  = 


6 al3 
-a.-  + — 

2a2 


a a a 64 

l , 3 

3a  4a"* 


k. 

l 


Case  2:  Suppose  that  the  cut  lies  inside 

one  degree  square  as  in  Figure  4-14. 
Similar  to  Case  1,  y « ax  ♦ 6 
represents  the  side  of  the  square 
that  intersects  the  interior  of 
the  cut.  Let  h(y>  represent  the 
density  function.  Upon  integra- 
tion in  both  directions  one  obtains. 


a portion  of  the 


Figure  4-14.  Case  2:  the  cut  lies 

inside  a portion  of  the  1*  square. 
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l 7 ’ 
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f'6 

I h(y)dy  = d(y6>  - H(ax+$) 


ax+0 


” a0  + V + V2  + a33 


£ , . 

| | • 

I 1 


a^,  a^,  a2>  a^  are  determined  constants 

f **  2 3 

Volume  -I  (ay+a^x+a  2X  +a3x  )dx  = k^ 

Case  3t  Suppose  a rectangular  region  lies  in  the  cut.  Then 
the  volume  is 


p h(y)dy  » a. 


•’I 


,(h,y  > 

txr,y  ) ,|,P  lx3,yp) 


a.dx  = a.x  -artx  = k. 
0 0 s 0 r i 


(x  ,y  ) (h,y  ) 
^ r | r 


(h7oi — * * 

Figure  4-15.  A rectangular 
region  lies  in  the  cut. 

Case  4 ; Suppose  that  a triangular  region  lies  in  the  cut.  The 
equations  of  the  two  lines  were  determined  as  in  Case  1.  Then 
the  volume  is,  after  integrating 
in  both  directions 


cix+6 


2 3 

h(y)dy  ® ag+a^+a^  4-a^x 


Volume  ® 


f'V*: 


2 3 

x+a2x  +a3x  )dx  ® k- 


<0,6) 


jy^ax+B 


(*ri.yr) 


/VJ 

(0,6)  7^x*6 


(x _,y  )L 1 

P P 

Figure  4-16.  A triangular  region  lies 
in  the  cut. 
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Case  5;  Suppose  the  region  appearing  in  Figure  4-17  lies 


completely  in  the  cut  on 
Then  upon  integrating  the  density 
function  h(y)  in  two  part,  one 
obtains 


Figure  4.17.  A polygonal  region  lies 
in  the  cut. 


/ax+B  f 

n(y)dy  + / 

0 J 


ax+8 

yx+6 


h(y)dy 


~2~3a  ~ -2 

a.  + a.x  + a~x.  + a.,x  + a„  + a.x  + a.x 


2 3 

= a^  + axx  + a^x  + a3x 


v A 


where  a.^  = (a^+a^)  i = 0,lf2,3 


Volume  = 


= j 4 ^a0+“ix+'  “-a3x3)dx  +J 


x 


r 2 3 

(aQ+a1x+a2x  -t-a^x  )dx=k^ 

~6/ 

Y 


Case  6:  Suppose  that  the  region  of  the  1°  square  lieu  in  the 

upper  portion  of  the  cut  as  in  the  following  figure 


Figure  4-18.  Case  6:  Region  of  the  1°  square  lies  in  upper 

portion  of  the  cut. 
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Given  the  density  function  h(y),  the  corresponding  volume  is 


i 


yx+6 


h(y)dy  + 


f6' 

Jax+B 


h (y)  dy 


2 3~  - ~2~3 

= aQ  + a^x  + a2x  + a^x  + a^  + a^x  + a^x  + a^x 


/: 


r 2 3 f 

Volume  = / (an+a.x+a„x  +a^x  )dx  + j 

u l 2 J jx 


2 3 

(aQ+a1x+a2x  +a3x  )dx 


= k. 

i 


Case  7:  Suppose  that  three  edges  of  the  square  lie  in  the  cut, 

as  given  in  the  following  figure 


Figure  4-19,  Case  7:  Three  edges  of  the  square  lie 

in  the  cut. 


Let  h(y)  denote  the  corresponding  density  fucntion 
Yt 


I 


ox+<$> 


6 h(vidy  + 6 h{y)dy  + h(y)dy 


'yx+6 


*'yx+6 
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Then  the  volume  is 


ve 


V 


= J p (aQ+a1x+a2x2+a3x3)dx  + f a (aQ+a1x+a2x2+a3x3)dx  + 

*V>  <*X._ 


z 


/\  A 


(aQ+a1x-ia2x2+a3x3)  dx 


a 


Case  8;  Suppose. that  the  whole  1°  square  lies  in  the  cut.  Then 
the  corresponding  volume  will  be  determined  in  the  following  manner 


y 


Figure  4-20.  Case  8 
The  whole  1°  square 
lies  in  the  cut. 


x 


© 


© 


© 


/ox+8  fCix+S  ^oix+B 

h(y)dy  + I h(y)dy  + I h(y)dy 

px+<}>  Jpx+<j> 


fx+6 


V 


= J 8 (a0+a3x+a2x2+a3x3)dx  +•  (ao+aix+a2x2+a3x3)dx  + 


f r (aQ+a1x+a2x2+a3x3)dx  « kj 
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Case  9:  Suppose  that  tne  following  region  lies  in  tne  cut. 

Then  the  volume  is 


y 


Figure  4-21.  Case  9:  A polygonal  region  lies  in  the  cut. 


1 


ox+6  *ox+6  *»ax+B 

h (y)  dy  + l h(y)dy  + / h(y)dy 

px+({>  J YX+6  Jyx+6 


Then  integrating  in  the  x direction. 


-L 


' 2 3 A 2 3 

Volume  = | p (aQ+a^x+a„x  +a3xJ)dx  + / (aQ+a1x+a2x  +a^x'3)dx  + 

0 J x 


f x — — _ o — 3 

I r (ag+a^x+a^x  +a3X  )dx 


= k. 

l 
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Case  10:  Given  the  region  below  the  volume  is* 


y 


Figure  4-22.  Case  10:  The  cut  lies  completely  in  the  1°  square. 


The  volume  of  the  rectangular  region  is 


[ yi  - v1 

Jo  h(yldy  - a0 


then  Volume 


i: 


a„dx  a a h 


The  volume  of  the  remaining  region  car.  be  delt  with  as  in  Case  6. 
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Case  XI:  Suppose  the  following  region  occurs  then  the  volume  is 


■»*; - V "> 


y=ax+e 


(Vyt>' 


y=px+<}> 


1 y=ax+B 

l\(x  #y  ) 
/ r r 


vv 

Figure  4-23.  Case  11:  A polygonal  region  lies  in  the  cut. 


Jf  ox+0  f ax+B  f ctx+B 

h(y)dy  + I h{y)dy  +/  h{y)dy 

px+<j)  J px+<f>  Jyx+& 


Then 


fyx+6 


/X  j*  X 

0 S(a0+a1x+a2x2+a3x3)dx  + J p (a0+a]x+a2x2+a3x3)dx  + 


— — — 2 — 3 

<a0+alX+a2X  +a3X  Jdx 
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Figure  4-24.  Case  12:  Three  regions  of  the  1°  square  lie  outside 

the  cut. 


f 


ax+6 

h(y)dy 

px+<t> 


/y  yax+e 

h(y)dy  + / h(y)dy  + / h(y)dy 
px+$  J yx+6  J yx+6 


Volunin 


/ 


y6-0 

° ^a0+alX+a2X  +a3X  + 


J-  x 

r p . . 
v°(a°+ai 


2 a ^ 

x+a2x  +a3x  )dx 


♦ 


a (an*a.xra  x2+a,x3)dx 


♦ 


1 


t aA+a, x+a,x2+a,x3 ) dx 


1 1 Now  all  the  densities  have  been  distributed  across  the 

route.  Then  to  find  the  total  volume  in  the  one  degree  square# 
I ' you  would  sum  up  all  the  components  comprising  it. 


Oil  and  Fishing  Rigs 

Next  the  question  of  distributing  the  Oil  Rigs  arises. 

The  main  data  sources  are  the  Register  of  Offshore  Units  Sub- 
mersibles  and  Diving  Systems  1977--1978*4  and  The  Ocean  Industry5. 

In  The  Ocean  Industry,  the  oil  rigs  are  categorized  according 
to  types:  submersibles , semi-submerisibles , jack-ups,  drillships 
and  barges.  They  also  provide  us  with  the  geographical  location 
in  which  they  are  presently  working  and  the  depth  capability  of 
the  rig.  This  information  is  also  given  in  the  Register  of 
Offshore  Units.  Therefore  by  correlating  the  two  sources  for  the 
geographical  location,  one  was  better  able  to  pin-point  the  location 
of  the  oil  rig.  According  to  Ocean  Industry  there  are  439  units 
in  existance  right  now.  This  figure  also  takes  into  account  those 
oil  rigs  under  construction,  and  those  available  for  contract. 
Therefore  the  number  of  oil  rigs  in  use  is  approximately  360. 

Most  of  the  oil  rigs  tend  to  lie  in  distinct  geographical  loca- 
tions, namely  the  Gulf  of  Mexico,  the  North  Sea,  the  Caribbean, 

Suez,  Red  Sea,  and  the  Persian  Gulf.12 

The  sources  of  data  for  the  fishing  vessels  are  as  follows: 

. 2 3 

Soviet  Fisheries  Investigations  in  the  Northeastern  Pacific, 

Fisheries  Statistics  of  the  United  States  1970/  The  Fish  Resources 
of  the  Ocean®  Yearbook  of  Fisnery  Statistics®0  1971  Atlas  of 
the  Living  Resources  of  the  Seas?  Trends  in  Merchant  Shipping15 
(1969-1980).  In  order  to  determine  the  actual  distribution  the 
total  number  of  fishing  vessels  and  the  thousand  of  tons  of  catch 
is  required.  In  order  to  illustrate  the  magnitude  of  fishery 
resources  and  their  preserit  state  of  exploitation,  the  available 
data  have  been  grouped  by  regions  in  accordance  with  the  break- 
down adopted  by  FAO  in  the  Yearbook  of  Fishery  Statistics,  1971  Maps 
show  the  catch  for  each  statistical  area  according  to  the  four 
main  groups:  demersal  or  bottom-living  fish,  pelagic  fish, 

crustaceans  and  cephalopoda,  and  whales. 
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The  count  for  the  total  number  of  fishing  vessels  is 
obtained  from  Lloyd's  Register  of  Ships.  There  are  presently 
21,261  fishing  vessels,  of  which  76%  of  them  ere  at  sea  yielding 
16,159  vessels.  The  total  world  catch  is  70,481,700  tons*  therefore, 
the  thousand  tons  of  fish  per  catch  per  vessel  is  2.43. 

The  Yearbook  of  Fishery  Statistics  designates  the  total 
number  of  thousand  tons  of  fish  for  each  geographical  area  in  the 
world.  By  using  this  information  we  are  able  to  determine  the  number 
of  vessels  in  each  ocean  area  necessary  to  procure  this  catch. 

Knowing  the  various  times  of  the  year  in  which  the  4 classes  of 
fishes  are  caught,  enabled  us  to  determine  the  seasonality  factor. 

The  fishing  vessel  distributions  are  done  for  the  winter  and  summer 
for  each  one  degree  square. 

Additional  data  was  obtained  by  NOSIC,  Naval  Ocean  Surveillance 
Information  Center.  The  Red  Fla^  Shipping  data  was  used  to  supplement 
the  Lloyd's  Register  tapes,  in  various  geographical  regions.  The 
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Outputs 

Availability  - Administrative: 
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1. 

Any  types  of  data  vhich  can  be  culled  from  the  description 
of  the  data  base  is  ava;,!  ?.ble  in  machine  readable  format.  The 
data  can  be  provided  upon  request  in  any  type  of  machine  readable 
format:  punched  cards,  magnetic  tape,  or  disc  pack.  A request 

for  data  should  be  sent  by  letter  to; 

Dr.  Louis  P.  Solomon 
Suite  600 

7900  Westpark  Drive 
McLean,  Virginia  22101 
Tel:  (703)790-5950 

2.  Availability  - Technical: 

The  technical  data  is  available  on  7 track,  300  bpi 
magnetic  tape.  All  the  numbers  will  be  available  in  E format. 

The  actual  distributions  and  their  descriptors  can  be  understood 
by  considering  not  only  the  printout  of  the  data,  but  also  in 
conjunction  with  the  body  of  this  report. 
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CHAPTER  VI 

Discussion  of  Validity  of  Results  and  Sample  Outputs 

The  question  concerning  the  validity  of  the  1°  x 1°  square 
distribution  is  the  next  problem  that  must  be  addressed.  Pre- 
sently, there  is  a paucity  of  data  available  exhibiting  shipping 
distributions,  which  restricts  the  number  of  cases  that  can  be 
checked.  However,  there  has  been  extensive  studies  of  shipping 

traffic  in  the  North  Pacific.  They  are:  Church  Anchor,  Aircraft 

21  2 7 

Surveillance  of  Shipping,  Church  Opal:  Surveillance  of. Shipping, 

and  the  RMS  study,  Navy  Interim  Shipping  Distribution^®'^® 

The  Church  Anchor  Exercise  presents  the  results  of  extensive 
aircraft  surveillance  of  merchant  vessels  at  sea  in  the  North  East 
Pacific  during  the  month  of  September  1973.  Sufficient  observa- 
tions were  made  to  enable  a statistical  analysis  of  the  data  to 
be  performed.  The  resolution  of  the  Church  Anchor  Exercise  is 
by  5°  squares.  The  data  for  HITS  by  5°  squares  in  the  same  area, 
as  that  of  Church  Anchor,  were  analyzed  and  compared.  The  region 
considered  is  between  30°  N to  55°N  and  125°W  to  180°W  as  shown 
in  Figure  6-1.  The  top  numbers  in  Figure  6-1,  of  each  5°  x 5° 
square  represents  the  number  of  ships  observed  in  the  Church  Anchor 
exercise.  The  number  directly  beneath  it  represents  the  standard 
deviation  of  the  observed  results,  and  the  bottom  represents  the 
HITS  results.  We  have  found  that  the  HITS  distribution  to  be 
statistically  consistent  with  that  of  Church  Anchor.  The  following 
observations  were  found: 

a.  30%  of  the  data  lies  with  0,25  standard  deviation  by 
5°  square. 

b.  58%  of  the  data  lies  wt-hin  0.75  standard  deviation  by 
5°  square. 

c.  79%  of  the  data  lies  within  1.0  standard  deviation  by 
5°  square. 

d.  87%  of  the  data  lies  within  1.25  standard  deviation  by 
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Standard  Deviations 


Figure  6-3 


Next  the  ratio  of  the  HITS  data  to  the  Church  Anchor  data 
was  considered.  We  can  observe  from  the  histogram,  in  Figure  6-2, 
that  a majority  of  the  data  ratios  lie  within  .75  and  1.25. 

This  indicates  that  the  HITS  results  are  remarkably  close  to 
those  of  Church  Anchor's. 

Further  analysis  was  done  comparing  the  Church  Anchor  to 
the  HITS  data.  In  Figure  6-3  we  have  a histogram  representing 
the  normalized  HITS  density  relative  to  Church  Anchor's  densities. 
This  was  obtained  by  the  following  method: 

HITS  Data  - Church  Anchor  data 
CA  standard  deviation  for  the  square 

This  demonstrated  that  the  data  behaves  like  the  normal  standard 
deviation  curve.  This  again  validates  the  accuracy  of  the  HITS 
results. 

The  Church  Opal  Exercise  performed  similar  missions  as 
that  of  Church  Anchor.  The  surveillance  of  merchant  vessels 
took  place  in  September  of  1573  in  the  North  East  Pacific.  The 
number  of  merchant  vessels  were  tabulated  and  analyzed.  Figure  6-4 
illustrates  the  findings  of  various  studies  for  5°  x 5°  square 
designated  in  the  North  East  Pacific.  Reading  from  top  to  bottom 
in  a square,  one  can  observe  the  results  of  the  HITS,  Church  Opal, 
Church  Anchor  and  RMS  findings.  The  HITS  results  compared  very 
well  to  that  of  Church  Opal  data. 
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The  correlation  coefficients  were  next  performed  on  the 
Church  Anchor,  Church  Opal,  RMS  and  the  HITS  data.  The  fol“ 
lowing  equation  was  used  to  determine  the  zero  based  correlation 
coefficient. 


2 _ (Sxiyi)‘ 
P°- 


Studies 


Church  Anchor-Church  Opal 
versus  RMS 


Church  Anchor-Church  Opal 
versus  HITS 


RMS  versus  HITS 


Church  Anchor  versus  HITS 


RMS  versus  Church  Anchor 


CHURCH  OPAL 
REGION 


.626 


CHURCH  ANCHOR 


Mid  Region  | Full  Region 


.724 

.543 

.839 

.799 

.898 

.785 

Table  6-1  Correlation  Coefficients  (p^)  among  various  studies. 


In  the  above  table  the  Church  Anchor-Church  Opal  entry  sig- 
nifies that  the  most  reliable  data  between  the  two  of  them  for  each 
square  was  chosen  for  the  Church  Opal  region.  One  can  observe  that 
the  HITS  data  correlates  better  in  the  Church  Opal  region  than  RMS 
when  compared  to  Church  Anchor  and  Church  Opal  data.  There  is  also 
a very  high  correlation  in  the  HITS  data  versus  the  Church  Anchor 
data  in  both  the  mid-region  and  full  region. 
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In  conclusion,  based  on  the  data  sources  available  and  the 
analysis  of  the  data,  the  HITS  results  have  successfully  modeled 
the  actual  shipping  densities. 

This  last  part  deals  with  some  of  the  sample  outputs  of  this 
study.  Figure  6-5  designates  the  routes  from  the  West  Coast,  U.S. 
to  Vladivostok,  Yokohama,  and  Shanghai.  The  number,  1539,  between 
the  route  Shanghai  and  Seattle  represents  the  ship  traffic  in  total 
number  of  ship  days  for  the  month  of  October.  Program  Gamma  pro- 
duced the  number  of  ship  days  of  traffic  from  the  Lloyd's  tape  be- 
tween each  of  the  69  canonical  ports  for  three  class  types.  This 
output  is  placed  in  a table  called  the  Proto-Routing  Table.  Next 
Program  Epsilion  further  reduces  this  69  x 69  array  to  ship  days 
between  canonical  ports  that  are  connected  by  a direct  route. 
Therefore,  Figure  6-5  represents  the  total  number  of  ship  days  in 
the  North  East  Pacific  between  the  illustrated  canonical  ports. 

The  next  step  is  to  distribute  the  number  of  ship  days  across  the 
route  by  1°  square.  This  is  done  in  Program  Theta.  This  is  shown 
in  Figure  6-6.  This  is  the  1°  x 1°  square  distribution  of  merchant 
ships  for  October  in  the  Aleutian  Islands.  The  distribution  is 
done  by  processing  one  route  at  a time  and  storing  the  densities 
in  a (180  x 360)  array  called  the  Merchant  Counts  Array.  This  is 
ayain  by  class  and  by  month. 

Figures  6-7  and  6-8  show  the  number  of  oil  rigs  and  fishing 
vessels  in  the  Gulf  of  Mexico.  The  fishing  vessel’s  distributions 
are  done  for  the  summer  and  winter,  while  the  oil  rigs  are  done 
annually.  The  oil  rigs  consist  of  jack-up,  submersible,  semi- 
submersible,  drill  ship  and  barges. 
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CHAPTER  VII 


Computer  Programs 
1.  Overview  and  Program  Flow. 

This  chapter  deals  with  the  various  aspects  of  the  programs 
implemented  in  this  project.  The  program  flow,  function  of  each 
program,  mathematical  logic,  and  various  checks  and  cases  are 
considered. 

The  flow  charts  on  the  following  pages  depict  data  set  and 
program  interrelationships. 
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Function  of  Each  Program 

Program  Alpha:  It  obtains  the  ORE  Routing  Table  (50x50) 

the  Mediterranean  Routing  Table  (24x24) , and  constructs 
the  total  Routing  Table  (69x69) . Prints  it  out  and  stores 
it  on  a data  set. 

Program  Beta:  Converts  <r  iginal  form  of  route  envelope 

trees  to  the  final  form  by  converting  degrees/minutes  to 
decimal  degrees,  calculating  edge  lengths,  inserting  back 
pointers,  and  collecting  leaf  and  route  identifiers  into 
two  indices.  These  indices  are  the  Index  to  Envelope  by 
Route  and  the  Index  to  Envelope  by  Siglum.  A separate 
output  tree  is  generated  for  each  season. 

Program  Beta  1:  Builds  Siglum  Table  and  Siglum  Indices 

(Siglum  Map  and  Inverse  Siglum  Table)  from  the  Routing 
Table,  except  for  moving  the  Multiple  Routing  Table  from 
the  Routing  Table  to  the  end  of  the  Siglum  Table. 

Program  Beta  2:  It  sums  the  edge  lengths  from  the  Route 

Envelope  Trees  (final  form)  to  build  the  Port  Distance 
Table.  It  determines  the  total  path  length  for  each  Siglum. 
It  must  run  once  for  each  season. 

Program  Gamma  1:  Constructs  the  Port  Hash  Table  from  the 

Port  Table.  The  Port  Marne  Abbreviation  is  converted  to 
Hash  Index  by  taking  its  S/370  numeric  value,  dividing 
by  2999,  and  using  the  remainder. 

Program  Gamma  2t  Builds  the  Estimated  Travel  Times  Table 
by  tracing  the  Siglum  Table  and  totaling  Port  Distances 
from  the  Port  Distance  Table.  Assume  a speed  of  12  knots. 
For  ships  traveling  from  one  canonical  port  to  itself,  a 
travel  time  of  0 days  is  recorded.  This  program  must  be 
run  4 times,  once  for  each  season. 
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Program  Gamma:  It  reads  the  Lloyds  Tapes  by  week,  in 

reverse  order  (i.e.,  week  52  is  first).  It  coinpafea  the 
tape  with  the  ships  from  the  previous  tape  (i.e.,  the 
following  week)  for  error  control.  Uses  the  (updated) 
Canonical  Port  Mapping  Table  to  find  the  two  associated 
canonical  ports.  If  three  ports  are  given  (i.e.,  an  in- 
termediate port) , determine  which  segment  the  ship  will  be 
on  during  the  week.  (If  it  is  on  both,  it  is  treated  as 
two  separate  voyages).  If  the  associated  canonical  port 
indices  for  egtess  and  entrance  ports  are  a and  6 re- 
spectively, and  the  type  is  Te (1, 3]  (corresponding  to 
the  various  ship  types),  then  increment  the  (a,$,T)-th 
component  of  the  Proto-Routing  Table  by  the  number  of  days 
in  that  week  which  the  ship  spent  at  sea.  After  processing 
all  weeks  for  that  month,  write  the  Proto-Routing  Table  on 
disk.  Zero  out  the  matrix  and  begin  to  accumulate  for  the 
"next"  (i.e.,  previous)  month. 

Program  Delta  1?  Provides  a visual  check  on  the  Route  Trees 
to  catch  transcription  errors.  It  plots  each  vertex  of 
the  given  route  and  designates  canonical  ports  by  a C and 
a non'-canonical  port  by  a V.  The  great  circle  path  is 
generated  and  displayed  between  each  pair  of  vertices. 

A 5 6 allowance  is  allotted  on  either  ends  of  the  route 
for  clarity.  The  program  was  used  for  each  of  the  four 
seasons  on  all  85  routes. 

Program  Delta  2:  Checks  the  canonical  port  mapping.  The 

rson-canonical  ports  are  plotted  using  symbols  to  denote  the 
various  associated  canonical  ports.  This  provides  a quick 
checking  procedure  to  determine  if  in  a given  area  there 
are  mote  than  one  canonical  port,  associated  to  it.  It  also 
spots  typographical  errors  occurring  in  the  port  table. 

An  example  of  this  is  the  occurrence  of  a port  in  the  middle 
of  a large  landmass.  The  program  is  run  in  two  parts,  once 
for  canonical  ports  1-34  and  is  repeated  for  canonical 
ports  .35-69. 


Program  Delta  3:  Checks  the  Routing  Table.  Prevents 

the  occurrence  of  infinite  loops,  and  does  various  checks 
on  non-direct  routes.  In  the  case  of  multiple  route  the 
check  was  done  manually.  When  an  error  is  spotted, 
a corresponding  message  is  written. 

Program  Theta  1:  Generates  a set  of  x values  to  be  used 

in  Theta  for  the  distribution.  It  takes  a set  of  data 
points  and  determines  the  Type  I'  Cubic  Spline  in  which 
it  interpolates.  Then  it  uses  two  iterations  of  Newton's 
Method  to  determine  intermediate  values  between  any  two 
successive  sextile  points.  Thus  producing  a set  of  13  x 
values.  This  is  done  for  each  season. 

Program  Epsilon:  This  program  reads  a monthly  Proto- 

Routing  Table  and  from  it  generates  a Ship-Days-Table. 

Program  Zeta:  Inputs  the  oil  rig  densities  in  the  one 

degree  squares.  It  produces  an  180x360  array  containing 
the  densities.  (Only  the  yearly  field  is  produced.) 

Program  Theta:  Takes  the  Ship-Days-Table,  created  by 

Epsilon,  applies  the  growth  factor  to  the  ship  populations 
and  converts  the  data  from  ship-days  per  route  to  ship- 
days  per  one  degree  square.  The  distributions  are  done 
between  canonical  ports,  for  the  various  types  of  vessels 
for  each  month.  The  array  is  known  as  the  Merchants  Count  Array. 

Program  Iota;  Reads  the  Merchant  Counts  Arrays  (180x360) 
for  all  the  various  types  of  vessels.  It  then  sums  them 
up  to  produce  the  seasonal  and  yearly  Merchant  Counts  Arrays. 
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Detailed  Table  Description 


a.  Port  Tablet  Consists  of  a listing  of  2,425  ports  that 
were  taken  from  Lloyd's  data  base  for  1972  by  the  Department 
of  Transportation'.  The  information  on  the  Port  Table  is  in 
the  form  of  Table  7-1  below. 


Table  7-1  Sample  of  Port  Table 


505 

DF.NF.R 

06O:RARA  6 HEWER ARA 

6N 

58K 

22 

440 

510 

CHV  H 

CEN  HHl DOR 

52N 

4F 

33 

44 

511 

Ca>A 

CfAMARK 

52N 

4F 

32 

44 

512 

CER1Y 

DERBY 

17S 

1236 

6 

34« 

513 

SCUM 

C1RINCE  fi  CERt.YCJO 

4C?J 

29E 

36 

1 

514 

06Sf  A 

PUfRTC  OESEACl 

47S 

6bW 

27 

217 

5 1*3 

C.r.TR0 

CETlU)fT  C OtlRCIT  R1V6° 

42N 

83W 

20 

56 

516 

cevcN 

CEVCNRCRT  (TASMANIA) 

415 

1466 

12 

4 

517 

OFVGNT 

OEVCNPCRT  (TA$RANI Al 

415 

146F 

12 

4 

5 1,1 

CMtf  I 

A EU  CNAHI 

24N 

546 

4C 

77 

| D 

I 0 


The  description  of  the  fields  are  as  follows:  the  port  number, 

the  abbreviation  of  the  port  name,  the  actual  name  of  the  port. 

On  the  right  side  we  have  the  latitude  and  longitude  of  the 
port's  location,  the  canonical  port  number  to  which  the  port 
is  associated  and  the  distance  between  the  port,  and  its  canonical 
port.  From  the  2,425  ports,  69  major  ports  were  extracted.  They 
are  listed  in  Table  7-2,  along  with  their  latitude  and  longitude. 

b.  Routing  Tables  Consists  of  a 69  x 69  matrix  entered  in  69 
groups  of  69  elements.  Each  group  spans  three  lines  with 
24  elements  each  on  lines  1 and  2,  and  21  elements  on  line  3. 

Each  entry  contains  either  a letter  code  or  a two  digit  number. 

If  a two  digit  number  appears  then  that  is  the  route  number 
that  connects  the  two  porta,  namely,  port  a and  port  b. 

If  a letter  appears  then  it  can  be  one  of  the  following; 


Table  7-2 
Canonical  Ports 


Number 

Name 

1 

Bombay 

2 

Colombo 

3 

Calcutta 

4 

Singapore 

5 

Sunda  Strait 

6 

Cape  Leeuwin 

7 

Manila 

8 

Shang-hai 

9 

Vladivostok 

10 

Yokohama 

11 

Torres  Strait 

12 

Wilson  Promontory 

13 

Honolulu 

14 

Papeete 

15 

Seattle 

16 

San  Francisco 

17 

Havana 

18 

New  Orleans 

19 

Straits  of  Florida 

20 

New  York 

21 

Montreal 

22 

Panama 

23 

Valparaiso 

24 

Punt a Arenas 

25 

Belem 

j 26 

pio  de  Janeiro 

1 

27 

Montevideo 

28 

Pent land  Firth 

29 

Fastnet 

30 

Bishop  Rock 

31 

Strait  of  Gibraltar 

‘ 32 

Marseille 

33 

Skagina  Odde 

I 34 

Napoli 

1 35 

Cape  of  Good  Hope 

I 36 

Leningrad 

Lat. 

18*  56‘N 
6°  57'N 
22°  33’N 
1*  16  N 
6*  94'S 
34r  32'S 
14p  35'N 
31°  14'N 
43°  06*  39“!3 
35°  27'N 
10°  33’S 
39®  10'S 
21°  18'  30"N 
17°  32*S 
47®  36'N 
37®  48'  30  "N 
23®  08 'N 
29®  57'N 
24°  25'N 
40°  42  *N 
45®  30’  10*N 
8®  53’N 
33®  02* S 
53®  10’S 
1®  27'S 
22®  53'  30®S 
34®  54'  30"S 
58®  42‘N 
51*  20' N 
49®  45'N 
35®  57'ii 
43*  19’  10*N 
S7*  48 'N 
40*  50'  19’N 
34*  22 'S 
59®  55  *11 


Long. 

72*  51'B 
79*  Sl'B 
88*  19  'E 
103®  5C5E 
105®  50' E 
115*  08'E 
120®  58'E 
121®  29'  30"E 
131°  53’E 
139®  39'E 
142®  08'E 
146®  26'E 
157®  52'  15  “W 
149®  34 'W 
122°  22*W 
122®  24  ’W 
82®  20'  30 "W 
90*  03'  3C"W 
83*  OC 
74°  01  '> 

73*  33*W 
79*  31 'M 
71*  38 ’W 
70*  54 'W 
48*  30'  20"W 
43*  U'W 
56*  13'  30*W 
3*  20 'W 
9*  36'W 
6*  35  »W 
5*  45'W 
5*  21*  17*E 
10*  44'E 
14*  16*  20"E 
18*  23 *E 
30*  lS'E 


L 

[ 

L 


I 


< 


Table  7-2 (CONTD) 


Name 

Lat. 

Long . 

Arkhangelsk 

64°  32 ‘N 

40°  31 'E 

Istanbul 

41®  01*  15"N 

28°  59'  09"E 

Port  Said 

31°  16 'N 

32°  19  'E 

Al  Basrah  . 

30°  31*  30“N 

47’-  50'  30  "E 

Nome 

64®  29 'N 

165°  25 'W 

Straits  of  Malacca  (NW) 

6°  50*N 

95°  E 

Tiksi 

71°  50  *N 

130°  E 

Cheluskin 

77®  50 'N 

105°  E 

Churchill 

58°  46  'N 

94°  ll'W 

Reykjavik 

64°  09 'N 

21°  56 'W 

Dakar 

14°  41 ’ N 

17®  26'W 

Guayaquil 

2°  12'S 

79°  53'W 

Petropavlovsk 

53’  N 

158°  38 ' E 

Nagayevo 

59°  33 'N 

150®  45’E 

Barcelona 

41°  23 'N 

2°  10 'E 

Genova 

44°  25 'N 

8®  57 ' E 

Livorno 

43°  31'N 

10°  18’E 

Messina 

38°  15’N 

15®  40 'E 

Oran 

35°  40  'N 

0°  39’W 

Strait  of  Bonifacio 

41°  18‘N 

9°  11*  E 

Tarabulus 

32°  52'N 

13°  12'E 

Alexandria 

31°  15'N 

29°  55' E 

Venice 

45°  28’N 

12°  23'E 

Trieste 

45°  40*N 

13°  45’E 

Piraievs 

37°  55’N 

23°  40 'E 

Thessaionica 

40°  38‘N 

22°  53’E 

Izmir 

38°  24 ’N 

27°  04 ' E 

Dardanelles 

40°  03'N 

26°  05'E 

Iskenderun 

36°  40*N 

36°  05  *E 

Beirut 

33s  58 'N 

35’  30'E 

Tel  Aviv 

32*  04*N 

34’  42 'E 

Strait  of  Otranto 

40°  17 'N 

18®  34'E 

Strait  of  Sicily 

37°  15'N 

11®  IS’E 
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• p means  that  ports  a and  port  b are  connected  by  an 
intermediate  port  namely  the  number  that  follows  P. 

• Blank  means  that  port  a and  port  1 tre  connected  by 
a direct  route  namely  the  number  following  it. 

• B means  that  a multiple  route  can  be  taken.  The 
number  following  the  B represents  the  line  number  in 
the  Multiple  Routing  Table.  It  will  consist  of  2 
entries,  the  first  number  is  the  first  intermediate 
port  to  be  taken  by  super-tankers  while  the  second 
number  is  the  first  intermediate  port  to  be  taken  by 
tankers. 

• N means  that  there  is  no  route  between  ports.  This 
arises  when  port  a equals  port  b. 

Below  is  a sample  listing  of  the  Routing  Table. 

Table  7-3 

r,  Q 69P  2P  2P  2P  2P  2P  2P  2P  2P  2P  2P  2P  2P  2P  2P39P39P39P39P39P39P35 Pi'S 
P39P19P39P39P  39P39P39P j ?P39 P39  61P39P39P39  58  65P  2P  2P39P39PJ9Pi9P3VP3‘i 
P 2P  2P 3 5P 3 9P  39 P39P 39 P39P39  P39  P39P39P39  P39P39  P39P39P39P39P39P39 
69K  0 62  59  59P  5P  4P  kP  4P  4P  4P  5P  4P  4P  4P42P39P39P39P39P39P39 P35P35 
P39P35P35P39P39  P39P39P39P39 P39  62P39P39P39  59  59P  4 59P39P39P39P 39P 39P39 
P 4P  4P3SP39*’  39P39P39P39P39P39P.J9P39P39P39P39P39k39P39P39P39P39 
P 2 £2N  0 29P  k-P  4P  4P  4P  <*P  4P  4P  4P  4P  4P  4P  4P  2F  HP  2P  HP  2P  2P35P35 

P 2P35P35P  2P  2P  2P  2?  2 P 2P  2 62P  2P  2P  2P  2P.2P42  29P  HP  2P  2P  2P  2P  2 

P 4P  4P  2P  2P  2P  2P  2P  2?  2P  2P  2P  2P  2P  2P  2P  2P  2P  2P  2P  2P  2 

P 2 59  ?9N  0 59 P 5 71  72»  7 28  27P  5 14P11P10P13P39P39P39P39P39  4PU?ll 

P39P  5P  5P39P39P39P39P35P3°P39  62P39P39P39  59  59P10  29P39P39P39P39p39P22 
Pl3P10P39P39P39P35P39P39P39P39P3np3qP39?39P39?39P39Pi9P39P39i»39 

P ? 55P  4 Tl.vj  3 25  72  72P  4P-  4 71P  6P  4PI  IP  4P  4P  4P  4P  4P  4P  4P  4PUP35 
P 2P35P35P39P 39P39P39P 39P39P39  63P39P39P39  63  59P  4P  4P39P39P39P39P39P  4 
P 4P  4P39P39P39Pi9P39P39P39P39P39P39P39P39P39P39P39P39P39P391*39 


c.  Route  Trees t There  are  four  versions  of  the  Route  Trees, 
one  for  each  season.  There  are  presently  85  routes  connecting 
the  69  pairs  of  canonical  ports.  The  table  contains  various 
information.  Field  one  contains  the  vertex  identifier  which 
is  in  binary  form,  field  two  contains  the  latitude  and  longi- 
tude of  the  vertex  identifier  in  decimal  degrees.  Southern 
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and  eastern  hemisphere  are  designated  by  a negative  sign 
and  North  and  West  are  positive.  The  third  field  represents 
the  vertex  width;  i.e.,  the  width  of  the  envelope  at  the 
vertex.  Field  four  is  the  edge  length,  which  is  the  length 
of  the  previous  route  segment  when  progressing  forward 
through  the  route.  Field  five  is  the  edge  width,  which  is 
the  widest  part  of  the  envelope  on  the  previous  route  seg- 
ment while  traveling  forward  through  the  route.  The  last 
field  is  the  back-pointers,  they  indicate  how  the  vertices 
are  linked  together.  When  progressing  backwards  along  a 
route,  the  back  pointers  indicate  the  line  in  the  Route  Trees 
to  which  the  present  vertex  is  connected.  Whenever  one  begins 
a new  route  the  backpointer  is  set  to  zero  and  the  identifier 
is  1.  Table  7-4  is  a sample  of  the  Route  Trees. 


Table  7-4  Sample  of  Route  Trees 


Line  I 


28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 


X 

10 

100 

lcdv 

101 

U 

no 

1U 

Ul-J 

lUod 

t 

10 
ICO 
1 l)v#0 
1000.) 
too . >0 
loOv  too 

10 v 00vl  v-J 

Luo.*  idoUO 
IdiwOCOOuO 
loo*  Wuacor 
1 

10 

ltd 

10-Jd 


Fields 


II  III  IV  V VI 

k*e»  JO*  l 7o  .uO  J.bxJ  GOO.OOO.wUCOOO 
!*4J.C.U*ioO.OO'*.do  JO/. 30 4. 030020 
. 20*124.  t6».  0 C 02  5 . 12  00  002  7 

♦ ‘W.v>0*122.3&b,00l'Ol.  7tO.OOX)OiO 

♦ 37.nyvl22.4OO.COjo2‘J.3Sp.Ooj0O2U 

♦ ^3.00~130.0C2,OU121<.  >/i.Ojjou2P 

*;*>.<*  j-L 0*3.6 ‘j  0.3 0 i 10.  7-J  3. 0oboo3 

♦ 31. 00- 131.00 3. 2*119.26  2.  )oJl033 

♦ j1.O0-UO.6J  J.^5303.  I .0.2 Ji5 
♦31.21-I2l.4e  0#U0  JOi)  .37jl  .Ooboic 
♦o  i . 2 5 ♦«  79  , 75  ),  00  J00. 00  3.  OOUOOO 

♦ 00.  JUjM.U  l.O3)U0.4  7l).lO|iOi6 

♦ J?.u0*0,i0. loO. 3 3002. o^l*  2*1)039  ; 
{♦Uit.9  I ♦ o J L • 4 1 J.  3 3 jO  l . 24  ,)«  > 3 J040  ' 
>+t/.3d*105.oUl,3di25.92U>0}04l  , 
♦10.91-U5.S3b.i2l?Ua^.oor.O^  j 
*l«r.t»3-U5.i>0  ).£5,.0O.33V,25iHl^j  j 
^0v.»l-125.333«*2  )u0.93p.l2  0044 

♦ 09. 03-12  J. 75  3-12  JOUsu0,  UU045  - 

♦ 0/,ou“lio.>U  J.i  j i07.2’*jo.i>J  )04i* 

♦ 3l»£o-t03. 63M.C0  j!4.  11>.1.00|)U4 1 - 

♦ 2L.3i,»l57.00b.JUl)OJ.OOU.OoUOOO  ; 
{♦5‘».i3*165.O0il.i2)a3.4Lj*.00bO4V  | 
.♦aO.OJ  + ltf  7.8t[l.ouba5. 8R2.34»05J 
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d.  Siglum  Table:  This  table  is  again  a 69  x 69  matrix 

entered  in  69  groups  of  69  entries.  Each  group  spans  4 
lines  with  18  entries  on  the  first  three  lines  and  15  entries 
on  the  last.  The  codes  for  the  entries  are  exactly  like  those 
of  the  Routing  Table  except  for  the  case  where  port  a and 
port  b are  connected  by  a direct  route?  i.e.,  no  intermediate 
stop.  In  this  table  the  route  numbers  have  been  replaced 
by  Siglia.  To  every  ordered  pair  of  canonical  ports  that 
are  linked  by  a direct  route  a different  siglum  is  assigned. 
Presently  there  are  501  siglia.  The  first  pair  of  canonical 
ports  that  are  connected  by  a direct  route  is  the  first 
siglum,  the  second  pair  is  siglia  number  2 and  so  forth. 

Table  7-5  is  a sample  of  the  Siglum  Table. 


Table  7-5  Siglum  Table,  'ample 


’i 

0 

IP 

SP 

2P 

cP 

IP 

2P 

ap 

2P 

ap 

2P 

ap 

ap 

2P 

SP 

ap 

39  P 

39 

?' 

39P 

33P 

39P 

39P 

35P 

35P 

39P 

3<*P 

39P 

39P 

39P 

39P 

39P 

39P 

39P 

39 

SP 

39 

p 

39P 

4P 

2P 

SP 

39P 

39P 

39P 

39P 

39P 

39P 

2P 

2P 

39P 

39P 

39P 

39 

39P 

3vP 

39P 

39P 

39  P 

39P 

39P 

39P 

39P 

39P 

39  P 

39P 

39P 

39P 

39 

0 

<6- 

7 

SP 

5P 

4P 

4P 

4P 

4P 

4P 

5P 

4P 

4P 

4P 

42P 

39P 

39 

c> 

39P 

39P 

39P 

35P 

35P 

39P 

35P 

35P 

39P 

3?P 

39P 

39P 

39P 

39P 

39 

9P 

•w*  . 

39P 

3<4 

1 0 

UP 

4 

12P 

39  P 

39P 

39P 

39P 

39P 

39P 

4P 

4P 

39P 

39P 

39  P 

39 

p 

3rP 

3 9P 

39P 

39P 

39P 

39P 

39P 

39P 

39P 

39P 

39P 

39P 

39P 

39 

*■ 

> 

c 

13H 

0 

14P 

4P 

4P 

4P 

4P 

4P 

4P 

4P 

4P 

4P 

4P 

4P 

4P 

SP 

2 

2P 

2P 

SP 

2P 

35P 

35P 

2P 

35P 

35P 

2P 

2P 

SP 

2P 

SP 

2P 

2 

15  p 

2 

p 

SP 

2P 

SP 

SP 

4a 

*6P 

ap 

£P 

2P 

ap 

2P 

SP 

4P 

4P 

SP 

£P 

ap 

2 

S* 

2P 

SP 

2P 

SP 

ap 

2P 

2P 

2P 

2P 

2P 

SP 

SP 

SP 

SP 

2 

f' 

u 

17 

1SN 

0 

19P 

S 

20 

2 IP 

7 

22 

23P 

5 

S4P 

UP 

t op 

13P 

39P 

39 

P 

39P 

39P 

39 

25P 

UP 

UP 

33P 

SP 

5P 

39P 

39P 

39P 

39P 

39P 

39* 

39 

£6P 

39 

P 

39P 

39 

2? 

2SP 

to 

29P 

39P 

39P 

39P 

39P 

39P 

SSP 

tOP 

I OP 

39  P 

39P 

39P 

39 

0 

3 ?P 

34P 

39P 

39P 

39P 

39P 

39P 

39P 

39P 

39P 

39p 

39  P 

39P 

39P 

39 

P 

3 

3 OP 

4 

31ft 

0 

3£ 

33 

34P 

4P 

4 

35P 

6P 

4P 

UP 

4P 

4P 

4P 

4 

c 

4p 

4P 

4P 

4P 

UP 

35P 

£P 

3SP 

35P 

39  P 

39P 

39P 

39P 

39P 

39° 

39 

36P 

39 

C' 

3^50 

39 

37 

3SP 

4P 

4P 

39P 

39P 

39P 

39P 

39P 

4P 

4P 

4P 

39P 

3SP 

39P 

39 

•J. 

39F 

39P 

3$P 

39P 

39 P 

39P 

39P 

39P 

39P 

39P 

39* 

39P 

39P 

39P 

39 

WHS  P4Q£  IS  BEST  iUAiXTT  P8lCtICJLB4J| 
KK*  wji  1 10  BOQ 
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Siglum  Indices: 

ml <— . -I...  - '■>  ■ . 


This  table  consists  of  two  sub-tables; 
Inverse  Siglum  Table  and  Siglum  Map.  The  Inverse  Siglum 
Table  is  headed  by  ’ INV#SGLM#TBL 1 followed  by  a list  of 
number  pairs.  Each  pair  corresponds  to  a successive  siglum, 
and  contains  two  (2)  digit  numbers  specifying  the  pair  of 
canonical  ports.  The  first  entry,  " 1 2",  says  that  ports 
1 and  2 are  connected  by  a direct  route  and  is  designated  as 
siglia  1.  There  are  presently  501  sigla.  The  Siglum  Map 
is  headed  by  ' SGLM#MAP*  and  it  contains  85  routes  each  re- 
quiring three  lines.  Each  group  contains  a header  entry 
plus  45  entries  (15  per  line) . The  header  entry  (the  first 
entry  of  each  group)  specifies  the  number  of  port  pairs 
linked  by  that  route  (or  equivalently,  the  number  of  sigla 
assigned  to  each  route) . Table  7-6  is  a sample  of  the  two 
tables. 


Table  7-6  Sample  Data  from  Inverse  Siglum  Table 
and  Siglum  Map 


ii. v /buL-um 


i ;■ 

lo5 

l 3-3 

l<tO 

2 1 

2 3 

2 t 

2 

5 

235 

2 JV 

240 

24c 

3 I 

3 M 

3 35 

342 

4 2 

4 3 

4 5 

4 

7 

4 0 

410 

411 

4 l 

h t-  _i 

435 

4 3V 

4 40 

4-tt 

5 2 

5 4 

5 

6 

3 7 

‘ii  ii 

511 

5-  l 

5-,0 

i>  'j 

o 7 

cl  2 

o 3 5 

bit 

7 

4 

7 3 

7 6 

7 0 

71- 

711 

»*  L j 

75  0 

C 4 

8 5 

t 7 

•J 

> 

UlO 

aii 

Ul3 

HU 

<•2  2 

050 

0 8 

5 15 

<7  **5 

i)  t 

10 

7 

10  8 

in  12 

Uli 

i 0 1 J 

1 0 *.  o 

1022 

10  41 

104  0 

1050 

11  4 

11  5 

1 t 

7 

UU 

1113 

1114 

U22 

i A 2 -> 

1 1 ■* 

11  35 

U 6 

UlU 

1211 

1213 

121 

t 

1222 

1223 

122  v 

4 2 33 

SGLrtyMAP 


12  303 

304 

305 

0 

G 

0 

0 

C 

0 

4 <34 

6 5 

117 

0 

C 

0 

0 

c 

0 

6 59 

oO 

72 

0 

0 

0 

0 

0 

0 

2 25 

iol 

0 

0 

0 

0 

0 

0 

0 

0 

0 

u 

0 

0 


0 

0 

0 

0 

3 


0 

c 

c 

c 

c 


0 

0 

0 

0 

0 


J 

0 

0 

0 

0 


o 

0 

0 

0 


0 

c 

c 

c 

c 


0 

J 

0 

0 

0 


*BI8  PAGE  tS  BEST  QUAIIIS  PSiCtt5**tfl 

jjgg  QQPY  fUMttSBSS  ^ *““*1 ^ 
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0 

0 

0 

0 


n 1 lo  tie  122  123 


0 

0 

0 
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0 


0 

0 

u 

0 

c 

c 

0 

0 

0 

0 

0 
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c 

0 

0 

0 

c 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


i53 

0 

0 

0 

0 

0 

0 

0 

13 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

J 

u 


0 

u 

0 

0 

0 

J 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

u 

0 
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f.  Port  Distance  Tablet  This  table  consists  of  6 pairs 
of  numbers  per  line.  The  first  number  in  each  pair  is 
the  number  of  the  Siglum  to  which  the  pair  corresponds. 

The  second  is  the  distance  (in  units  of  latitudinal  degrees) 
between  the  pair  of  canonical  ports  corresponding  to  that 
siglum.  There  are  4 Port  Distance  Tables,  one  for  each  season. 
Below  is  a partial  listing  of  the  Port  Distance  Table  for 
the  fall. 


Table  7-7  Sample  Data  from  Port  Distance  Table 
Fall 


1 

15.03 

2 

78.3  & 

3 

51.29 

4 

26.23 

5 

15.03 

6 

20.94 

7 

27.44 

8 

35.46 

9 

86.35 

10 

57.07 

11 

42.08 

12 

17.  03 

13 

20.94 

14 

2 7,o2 

15 

101.93 

L t> 

17.20 

17 

27,44 

Id 

2 7 ,n? 

19 

B. 02 

20 

22.03 

21 

3 7.69 

22 

4 7.79 

23 

43.21 

24 

9d . 16 

25 

175.32 

26 

94.96 

27 

34.51 

28 
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g.  Estimated  Travel  Time  Table:  This  table  is  a 69  x 69 

matrix  where  each  entry  corresponds  to  every  possible 
canonical  port  pair,  whether  linked  by  a direct  route  or 
not.  These  entries  are  arranged  in  69  groups  of  69  entries 
spanning  4 lines.  The  entries  represent  the  expected  travel 
time  in  integral  days  betweeen  the  two  ports.  If  the  entry 
is  negative,  the  absolute  value  of  this  number  is  an  index  to 
the  Multiple  Travel  Time  Table,  which  follows  immediately. 

This  table  consists  of  a pair  of  numbers,  giving  the  estimated 
travel  time  for  both  super-tankers  and  non-super-tankers* 

This  table  is  also  seasonally  dependent,  one  for  each  of  the 
four  seasons.  Table  7-P  is  a sample  of  the  Estimated  Travel 
Time  Table  for  the  spring. 


Table  7-8  Sample  of  Estimated  Travel  Time  Table 
(Spring) 
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h.  Proto- Routing  Table:  This  table  consists  of  three 

69  x 69  arrays?  the  first  for  merchants,  the  second  for 
tankers  less  tahn  70,000  gt,  and  the  third  for  tankers 
greater  than  70,000  gt.  The  first  line  indicates  which 
month  is  being  processed,  for  example  in  Table  7-9  we  are  pro- 
cessing Week  48  i.e.,  month  December.  The  second  line  indicates 

the  ship  type;  for  example;  Tau  = 1 for  merchants.  The  69  x 69 
elements  spans  5 lines.  The  entries  represent  the  number  of 
ship  days  traveled  between  any  pair  of  canonical  ports  of  a 
given  ship  type  in  a given  month.  There  exists  12  Proto 
Routing  Tables,  one  per  month.  Table  7-9  is  a sample  of 
Proto  48. 


Table  7-9  Sample  Data  from  Proto  48 
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i.  Ship  Days  Table:  The  format  is  similar  to  that  of 

The  Proto  Routing  Table;  namely,  that  it  consists  of 
three  tables  each  having  501  entries.  These  are  again 
three  different  ship  types,  TAU  = 1,  TAU  a 2,  TAU  = 3. 

The  elements  are  listed  according  to  sigla;  i.e.,  canonical 
ports  linked  by  a direct  route,  which  represent  the  total 
number  of  ship  days  for  the  sigla;  i.e.,  between  the  pairs 
of  canonical  ports.  Following  it  lies  the  Residual  Table, 
again  by  ship  type.  The  Residual  Table  stores  the  number  of 
ship  days  to  be  distributed  for  local  traffic  and  'non- 
route' traffic  about  the  canonical  ports.  Table  7-10  is  a 
sample  of  the  ship  days  table  for  month  48;  i.e.,  December 
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Tabla  7-10  Sample  Ship  Days  Table 
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CHAPTER. VIII 


Conclusions  “ 

This  report  provides  shipping  distribution  by  1°  square 
throughout  all  the  ocean  areas  in  the  world.  The  distributions  are 
temporally  divided  by  month,  season,  and  annual  average.  The  ships 
are  broken  down  into  5 classes:  merchants,  tankers,  super  tankers, 
fishers,  and  oil  rigs. 

The  temporal  variation  is  clearly  area  dependent.  Some 
areas,  such  as  the  North  Pacific  are  very  seasonal  dependent; 
other  areas,  such  as  the  Indian  Ocean  are  almost  seasonal  inde- 
pendent. 

Fishing  vessels  are  distributed  for  only  winter  and  summer. 
These  distributions  are  based  on  nutrient  distribution  and  area  catch 
tonnage.  There  is  no  raw  data  describing  at-sea  positions  which 
allowed  a refinement  for  densities. 

The  merchant  and  tanker  distributions  are  heavily  based  upon 
Lloyd's  Shipping  Index  for  1972,  with  growth  factors  applied  to 
update  the  port  traffic  to  1978.  The  1972  data  base  was  rather 
spotty  in  inherent  quality.  Considerable  care  was  taken  in  polish- 
ing and  editing  the  individual  data  entries.  A new  data  base  of 
better  quality  would  provide  better  absolute  numbers;  however,  the 
basic  shipping  patterns  would  almost  assuredly  remain  unchanged. 

Multiple  sources  were  used  for  trend  and  growth  analysis  of 
the  1972  data  base.  The  growth  is  definitely  geographically  depen- 
dent. The  merchant  vessel  population  increased  rather  slowly; 
tanker  population  increased  by  60  percent  from  1972  to  1978.  Both 
populations  changed  in  character;  i.e.,  ship  type,  size,  tonnage, 
speed,  power  plant,  and  construction  technique. 


The  predicted  data  base  was  verified  by  comparison  with 
results  of  naval  exercises.  Agreement  within  a rather  large  area 
of  the  North  Pacific  is  within  one  standard  deviation  of  the 
measured  result.  Additional  comparisons  should  be  made  as  more 
statistically  reliable  data  becomes  available.  Until  th^t  time  it 
is  believed  that  the  results  are  easily  correct  within  a factor 
of  two  (2)  and  can  be  used  with  confidence  for  a number  of  analyses. 
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APPENDIX 


Program  Documentation 


A.  PROGRAM:  ALPHA  (RTNGMRGE) 

REFERENCE  PARAMETERS  - PROGRAMS; 

DSNAME  * CN5054.PS1. ALPHA 
VOL  = SER  = PUBL03 
UNIT  = 3330 

DCB  * (BLKSIZE  = 3120,  LRECL  = 80,  RECFM  = FB) 

Sequenced  73/80.  (Wylbur  line  numbers) 

REFERENCE  PARAMETERS  - JCL  and  CONTROL  DATA  SET: 

DSNAME  ® CN5054 .PS 1. JCL (ALPHA) 

VOL  ■ SER  ® PUBL03 
UNIT  =3330 
WYLBUR  EDIT  format 

FUNCTION:  Merges  the  original  Routing  Table  and  the  Mediterranean 

Routing  Table  to  form  the  final  routing  table.  It 
determines  the  correspondence  between  the  Mediterranean 
route  and  port  numbers  and  those  used  in  the  final 
routing  table.  The  program  is  written  in  a generalized 
form  to  allow  use  for  other,  similar  applications 
with  different  numbers  of  points,  routes,  and  multiple 
routing  tables,  and  also  difference  function  points  between 
the  original  and  additional  routing  tables.  Similar 
routes  and  ports  from  additional  and  original  routing 
tables  are  over lay ed. 

INPUT  DATA  SETS : 


FILENAME 


DSNAME 

RTNGTBL.ORG 


RTNGTBL.ADD 


CONTENTS 

Routing  Table  (original  version) 
including  Muliple  Routing  Table. 
Routing  Table  (Mediterranean 
additions) , including  Multiple 
Routing  Table  (if  any) . 


'tt-i.-q  v.w'i/Ti'-wvi  i ^r^rjiPrp 


FILENAME 


DSNAME 

RTOGTBL.ADD 


SYSIN 


CONTENTS 

(Both  ORG  and  ADD  may  contain 
any  number  of  entries  per  line, 
provided  that  the  LRECL  specified 
by  JCL  or  Volume  Table-of-Contents 
is  exactly  equal  to  the  line  length 
with  no  extra  padding  blanks  at  the 
end.  Each  run  of  the  routing 
table  must  start  on  a new  line; 
in  that  case  only  the  preceding 
line  may  be  truncated.  The  format 
is  otherwise  under  control  of  the 
Control  Data  Set  parameters,  (dis- 
cussed below) . The  multiple 
Routing  Tables  should  immediately 
follow  the  main  routing  tables,  with  one 
pair  per  line  and  no  spaces 
between  them,  with  the  same  number 
of  digits  as  specified  for  theimain 
routing  table. 

Control  Data  Set  Parameters  * 

ORG#NO#PORTS  (Number  of  ports 
in  the  original  routing  table) 

ORGINO# ROUTES  (Number  of  routes 
in  the  original  routing  table) 
ORGINOlMRTS  (Number  of  multiple 
routing  tables  in  the  original 
routing  table) 

ORGINOIDIGITS  (Number  of  digits 
in  the  routing  table  used  for  recording 
ports,  routes,  and  multiple-routing- 
table  numbers) 

ADDlNOtPORTS  (Number  of  ports  in  the 
additional  routing  table) 


W 


FILENAME  DSNAME  CONTENTS 

SYSIN  - ADD #NO# ROUTES  (Number  of  routes 

in  the  additional  routing  table) 
ADD#NO#MRTS  (Number  of  multiple 
routing  tables  in  the  additional 
routing  table) 

ADD#NO#DIGITS  (Number  of  digits 
in  the  additional  routing  table 
used  for  recording  ports,  routes, 
and  multiple-routing  table  numbers) . 
NEW#NO#PORTS1  (Number  of  port  for 
which  space  should  be  allocated 
in  the  new  routing  table) . 
NEW#NO#DIGITS  (Number  of  digits 
to  be  used  in  the  new  routing  table 
for  recording  ports,  routes,  and 
multiple  routing  table  numbers) 

These  parameters  are  read  in  Data-directed  format. 

This  is  followed  by  a list  giving  the  port  in  the  original  routing 
table  which  is  closest  to  the  port  in  the  additional  routing  table. 
There  is  one  line  for  each  port  in  the  additional  routing  table, 
and  two  entries  per  line*  The  first  giving  the  number  of  the 
closest  port  in  the  original  routing  table,  the  second  a code 
(enclosed  in  single  quote)  which  is  'IDENT'  if  ports  are 
identical,  and  ' ADJCN'  if  they  are  merely  adjacent.  Otherwise, 
this  list  may  be  in  free  format. 

This  is  followed  by  a list  of  number  pairs  specifying  which 
routes  from  the  additional  routing  table  are  to  be  overlayed  on 
routes,  from  the  original  routing  table.  The  first  element  of 
each  pair  gives  the  additional  route  number,  and  the  second  the 
original  number.  The  list  is  terminated  by  two  ‘O's.  At  least 
one  space  or  a comma  must  be  inserted  between  entries,  otherwise, 
this  list  may  be  in  free-form  format. 

iMay  be  omitted  ifTequal  to  2. 


A- 3 


OUTPUT  DATA  SETS; 

FILENAME  DSNAME 

NEW  RTNGTBL.NEW 


SYSPRINT  SYLIB  (PORTNO) 


CONTENTS 

New  routing  table,  along  with  the 
multiple  routing  table. 

Error  messages,  along  with  the  table 
of  correspondence  between  addi- 
tional and  new  port  numbers,  and  the 
multiple  routing  table  numbers 
(original  numbers  are  carried 
over  unchanged) . 


JCL  USAGE;  Assumes  that  a load  module  has  been  placed  by  the 

linkage  editor  specifying  the  NCAL  option  in  LINKLIB 
(ALPHA)  and  is  now  to  be  reprocessed  by  the  loader. 


MAJOR  VARIABLES; 
PORTtINDX. NEW#PORT#NO 


PQRTIINDX. ADD#PORT#NO 


ROUTE#  ?NDX . NEW#  ROUTE  #NO 


PORT#MAP 


ORG*  RTNGITBL,  ADD#RTNG#TE  !• , 
NEWIRTNGITBL 


CONTENTS 

An  array  referenced  by  an  additional  port 
nuirber  containing  corresponding  new 
port  numbers. 

Array  referenced  by  new  port  number 
containing  corresponding  additional 
port  number,  or  *0'  if  none. 

Array  referenced  by  additional  route 
number,  giving  number  of  corresponding 
route  in  new  routing  table. 

Array  referenced  by  additional  port 
numbers  containing  port  number  of 
closest  original  port,  and  code 
(either  "IDENT*  or  'ADJO?') 

Two  dimensional  arrays  referenced  by 
original,  additional,  and  new  port 
numbers,  respectively,  containing 
routing  tables. 


A-4 


CONTENTS 

Array  used  to  hold  multiple  routing 
tables,  referenced  by  new  multiple 
routing  table  number. 

Various  temporaries,  used  while 
building  tables. 

PROGRAM  STRUCTURE 

LINES  FUNCTION 

5/35  Reads  in  control  parameter,  checks  for  consistency 

38/141  Reads  in  list  of  original  ports  closest  to  additional 

ports,  and  determines  the  correspondence  between  addi- 
tional ports,  routes,  and  multiple  routing  table 
numbers  and  those  used  in  the  new  routing  table. 

Builds  PORT# INDEX  and  ROUTE# INDEX  for  later  use. 

145/368  Reads  in  original  and  additional  routing  tables, 

builds  new  routing  table,  and  writes  it  out.  (More 
details  are  provided  by  comments  in  the  program.) 

PROGRAM  MESSAGES:  (apart  from  new/additional  routing  table 

correspondence  list) 

’SEVERAL  ADDITIONAL  PORTS  HAVE  BEEN  IDENTIFIED  WITH  THE  ORIGINAL  PORT 
NO  XXX  INDICATING  ERROR  IN  PORT  MAPPING  SECTOR.  EXECUTION 
TERMINATING. ’ 

Response:  Correct  control  data  set 

•INVALID  LINKAGE  CODE  IN  PORT  MAPPING  VECTOR,  ADD  PORT  NO  XXX 
EXECUTION  TERMINATING. ' 

Response : Correct  control  data  set 

SYSTEM  MESSAGES:  (Probable  meaning)  0RG 

•CONVERSION'  CONDITION  RAISED,  'ONFILE'  » *ADD 

Response:  Make  sure  that  the  control  data  sot  parameters  corres- 

pond to  the  actual  format  of  the  input  routing  tables. 
Make  sure  LRECL  of  input  routing  tables  allows  for 
no  blanks  at  end  of  each  line  (except  for  truncated 
lines  at  end  of  the  row) . Checks  the  numeric  fields 
for  invalid  character  in  the  input  routing  tables. 

A-S 


MAJOR  VARIABLES  t 
MRTITBL 


( RTNG#ELEMENT , 
RTNG#ELMNT#1  (#2, #3) 
MRT#ELMNT) 


COMPLETION  CODE  SZZZ  TIME  LIMIT  EXCEEDED 

Responses  Checks*  for  infinite  recurive  loop  in  additional 

routing  table  that  are  triggered  in  lines  257/351 


A- 6 


B 


PROGRAM:  BETA1 


REFERENCE  PARAMETERS  (SOURCE  PROGRAM) 

DSNAME  * CN5054 .PS1.0ETA1 
VOL  = SER  = PUBL03 

DCB  = (RECFM  = FB,  BLKSIZE  = 3120  LRECL  = 80) 

Sequenced  73/80  WYLBUR  line  numbers. 

REFERENCE  PARAMETERS  (LOAD  MODULE) : 

DSNAME  = CN5054 .PS1.LINKLIB (BETA1) 

VOL  = SER  = PUBL03 
UNIT  = 3330 

Load  module  may  be  invoked  by  //  EXEC  PGM  = BETA1 

REFERENCE  PARAMETERS  (JCL),  (For  load  module  invocation) 

DSNAME  * CN5054:PS1.JCL(BETA1)  , 

VOL  a SER  = PUBL03 
UNIT  * 3330 
Wylbur  EDIT  format 

FUNCTION:  Builds  Siglum  Table  and  Siglura  Indices  (Siglum  Map 

and  Inverse  Siglum  Table)  from  Routing  Table,  except 

for  moving  the  Multiple  Routing  Table  from  the  Routing 

Table  to  end  of  the  Siglum  Table,  which  must  be  done  by  hand. 

INPUT  FILES; 

DDNAMS  DSNAME  CONTENTS 

RTNGTBL  CN5054 . PS1.  RTNGTBL.NEW  Final  version  of  Routing 

Tables,  must  be  precisely 
in  format  specified  in 
table  descriptions 

OUTPUT  FILES: 

SGLKNOX  CN5G54 . PSi. 3GLHNDX  Inverse  Siglum  Table  and 

Siglum  Map 

SGLHTBL  CN5054.PS1.SGLKTBL  Siglum  Table 


MAJOR  VARIABLES: 

INV#SGLM#TBL 

SGajM#MAP 

SGLM#TBL,RTNG#TBL 
PROGRAM  PARAMETERS: 


Inverse  Siglum  Table 
Siglum  Map 

Routing  Table,  on  input,  converted 
to  Siglum  Table  in  place. 


The  following  parameters  are  set  implicitly  by  the  program 
structure  and  will  necessitate  program  modification  if  changed. 


PARAMETER 

Number  of  Canonical 
Ports  (69) 


REFERENCES  (LINE  NUMBERS) 
24.27  Dimension 
52,74  I/O  Format 
54,  DO  loop  limit 


Number  of  Canonical  Routes 
(85) 


Maximum  Number  of  Sigla 
(600) 

Maximum  Number  of  Sigla  per 
Route  (45) 


19  Dimension 

20.21  Repetition  factor  in  INITIAL 
attribute 

75,4  I/O  Format 

72.22  DO  Loop  Limit 
16  Dimension 

21  Dimension 

21  Repetition  factor  in  INITIAL 
attribute 

75.4,75.5  I/O  Format 
(Note:  overflow  of  this  number  may 

be  eesily  detected  by  examining  the 
Siglum  Map.  The  first  entry  per  route 
is  the  actual  number  of  the  Sigla 
for  that  route.  If  this  is  greater 
than  45,  this  parameter  must  be 
increased. ) 


A- 8 


PROGRAM  STRUCTURE: 
LINES 


FUNCTION 


50/52  Reads  in  the  Routing  Table. 

54/60.6  Converts  all  entries  flagged  as 

representing  direct  routes  to 
Sigla,  building  simultaneously  the 
Inverse  Siglum  Table  and  Siglum 
Map. 

601/end  Writes  out  Siglum  Table,  Inverse 

Siglum  Table,  Siglum  Map 


PROGRAM  MESSAGES : 'None 

SYSTEM  MESSAGES  (WITH  PROBABLE  CAUSE) : 

'CONVERSION'  CONDITION  RAISED,  ' ONFILE 1 = RTNGTBL 

Response:  Check  the  Routing  Table  for  invalid  code.  Make 

sure  that  the  format  corresponds  to  that  given 
in  cable  descriptions. 
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c. 


PROGRAMS  BETA 


REFERENCE  PARAMETERS  - SOURCE  PROGRAM; 

DSNAME  = CN5054 . PS 1. BETA 
VOL  = SER  = PUBL03 
UNIT  = 3330 

DCB  = {RECFM  = FB,BLKSIZE=3120,LRECL=80) 
sequenced  73/80  with  WYLBUR  line  numbers. 

REFERNECE  PARAMETERS  - LOAD  MODULE; 

DSNAME  = CN5Q54 .PS1.LIMKLIB  (BETAF)  Note  Suffix 
VOL  = SER  = PUBL03 
UNIT  = 3330 

Load  module  read  for  execution  via  PGM  parameter  of 
EXEC  card;  no  link-editing  necessary. 

REFERENCE  PARAMETERS  - JCL 

DSNAME  = CN5054.PS1.  JCL  (BETA) 

VOL  - SER  = PUBL03 
UNIT  = 3330 
Wylbur  EDIT  format 

FUNCTION:  Converts  the  original  form  of  the  route  envelope  trees 

to  the  final  form  by  converting  degrees/minutes  to  decimal 
degrees,  calculating  edge  length,  inserting  back 
pointers,  and  collecting  leaf  and  route  identifiers 
into  two  indices,  the  Index  to  Envelope  by  Route  and 
the  Index  to  Envelope  by  Siglum.  A separate  output 
tree  is  generated  for  each  season. 

As  a secondary  by-product,  it  effectively  compares 
all  direct-route  linkages  in  the  routing  table  through 
the  intermediary  of  the  Siglum  Map  and  Inverse  Siglum 
Table  to  the  associated  route  envelope  trees-  It  prints 
an  error  message  if  a branch  of  a route  tree  is  never 
referenced  by  the  route  table.  Inspection  of  the 
Index  to  Envelope  by  Siglum  will  also  reveal  any 
direct- route  linkage  for  which  no  associated  branch 
in  the  route  Envelope  Tree  exists,  as  the  indicated 
line  number,  will  be  *999. 99*. 
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A ’999.99’  for  latitude  or  longitude  in  the  output 
route  envelope  trees  means  that  a vertex  has  been 
omitted  for  some  season  in  input. 


INPUT  DATA  SETS: 

DDNAME 

DSNAME 

SGLMNDX 

CN5054.PS1. SGLMNDX 

INENV 

CN50  54 . PS1  .ENVELOP] 

OUTPUT  DATA  SETS: 

ENVINDX 

CN5G54.PS1. ENVINDX 

OUTENVP 

CN5Q54. PS1 .OUTENV . 

SPRING 

OUTENVS 

CN5054.PS1 .OUTENV . 

SUMMER 

OUTENVF 

CN5054.PS1. OUTENV. 

FALL 

OUTENVW 

CN5054.PS1. OUTENV. 

WINTER 

CONTENTS 

Inverse  Siglum  Table  and  Siglum  Map 
Original  Version  of  Route  Trees 


Index  to  Envelope  by  Route  and 
Index  to  Envelope  by  Siglum 
Final  version  of  Route  Envelope, 
one  per  season. 


JCL  USAGE:  Assumes  that  the  load  module  has  been  placed  in  CN5054.PS1. 

PROGRAM  PARAMETERS: 

The  values  of  the  following  parameters,  are  set 
implicitly  via  program  structure.  The  program  must 
be  revised  and  recompiled  if  they  are  changed: 

PARAMETERS  REFERENCES  (POSSIBLY  NOT  COMPLETE)  BY  LINE 

Maximum  I of  routes  (85}  13,48  Dimension 

14  Repetition  factor  in  INITIAL  attribute 
78,245  Repetition  factor  in  FORMAT 


PARAMETERS 

Maximum  # of  Siglum  (506) 


Maximum  # of  Sigla  per 
route  (45) 

Maximum  # of  Output 

Vertices  per  Route  (75) 


REFERENCES  (POSSIBLY  NOT  COMPLETE)  BY  LINE 
15,46  Dimension 

16  Repetition  factor  in  INITIAL  attribute 
(2X  number  of  Sigla) 

77,246  Repetition  factor  in  FORMAT 
50  Dimension 

78  Number  of  lines  read  by  FORMAT  from 
Siglum  Map  (3) 

16.3,17  Dimension 

16.5,16.6  Repetition  factor  in  INITIAL 
attribute 


PROGRAM  MESSAGES: 

'PORT  W ROUTE  ENVELOPE  DOES  NOT  OCCUR  IN  ROUTING  TABLE,  ROUTE 
NO  = XXX  YYYYY ' 

XXX  is  Route  number 

YYYY  is  print  out  of  leaf,  as  taken  from  original  form  of  Route 
Tree . 


Message  may  be  printed  twice  for  the  same  leaf,  if  it  is  marked 
in  the  input  Route  Trees  as  being  used  for  both  entrance  and  egress 
and  is  used  in  Routing  Table  for  neither  purpose. 


Responses  Correct  Route  Trees  by  eliminating  unused  branch  of 
the  route,  or  correct  the  Routing  Table  by  finding 
some  use  for  the  route,  or  ignore  it  if  certain 
that  the  error  is  harmless,  (Note:  Route  21  and 
and  one  branch  of  route  12  are  currently 
unused) . 

SYSTEM  MESSAGES  AND  LIKELY  MEANINGS: 

IBM  7801:  NO  'OTHERWISE'  CLAUSE  AND  NO  'WHEN'  CLAUSES  SATIS XFED. 

Response:  Check  for  Season  Code  other  than  'P',  'S',  'F', 

'W' , or  Hemisphere  code  other  than  'N',  'S',  *E', 

'W'  in  the  input  Route  Tree  Envelope. 
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IBM  5371  DATA  EXCEPTION 

(Program  uses  decimal  arithmetic  in  many  calculations.) 

Response:  Check  for  illegal  code  in  the  alphanumeric  field  in  the 

input  route  tree  envelopes.  Also  make  sure  input 
matches  picture  specifications  in  lines  27/37, 


43/45.2. 


?’ TRANSMIT'  CONDITION  RAISED, 


' ONFILE ' 


,INENV 

^SGLMNDX 


(DCB  parameters  are  specified  within  program  and  cannot  be  modified 
by  JCL  or  volume  table  of  contents) 

Response:  Check  that  INENV  or  SGLMNDX  are  associated  with 

data  sets  with  LRECL  - 72,  BLKSIZE  = 3096.  If 
not,  reformat  data  set  (it  may  have  been  accidently 
stored  in  WYLBUR  EDIT  format)  by  using  it  and 
saving  it  with  LRECL=72 . 


'CONVERSION*  CONDITION  RAISED,  'ONFILE'  = SGLMNDX 


Response:  Make  sure  Inverse  Siglum  Table  contains  no  more 

than  the  maximum  number  of  Sigla  (currently  506) 
expected  by  program  BETA.  If  it  contains  less, 
pad  with  extra  blanks  or  blank  lines  to  simulate 
506.  (Note:  If  a terminal  error  occurs  in 

processing  a route,  the  envelope  output  file 
will  contain  all  previously  processed  routes.) 


Contain  Envelope  Indices 
Contain  Siglum  Indices 

Alternative  dummy  variables  describing 
record  from  the  input  route  envelope  trees, 

overlayed  upon  input  buffer,  depending 
upon  whether  record  is  vertex  descriptor 
or  route  header. 


MAJOR  VARIABLES: 

INDX  ENVLP_BY  ROUTE  , 
INDX~ENVLP_BYJSGLM 

INV_SGLM_TBL , 

SGL_MAP 

IN  VERTEX, 
IN~ENVELOPE 
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OUT_ENVELOPE  Array  used  for  storing  the  final  forms  of  the 

route  envelope  trees.  It  holds  the  verticles 
corresponding  to  all  four  seasons,  but 
one  route  only  (it  is  written  out  after 
processing  each  route) 

OUT_ENVELOPE_INIT  Array  used  to  reinitialize  OUT  -ENVELOPE 

at  the  start  of  processing  each  route, 
setting  latitudes  and  longitudes  to  the 
impossible  '999.99' 

OUT_VERTEX_BUF  Dummy  variable  overlayed  on  OUT-  ENVELOPE 

to  facilitate  output. 


PROGRAM  STRUCTURE 

LINES 

75/78 


101/108 


FUNCTION 

Reads  in  Siglum  Indices 
Establishes  ENDFILE  on-unit  to  terminate 
processing  after  all  input  trees  have 
been  read. 

If  this  is  the  beginning  of  processing  a 

new  route,  write  out  envelope  tree  for 
preceeding  route.  Record  vertex  number 
in  INDEX  - TO  ENVELOPE  - BY  - ROUTE 
and  re-initialize 


109/144 


144/172 


Check  if  this  vertex  is  a leaf  corresponding 
to  a canonical  port,  and  if  so,  update 
INDEX  TO  ENVELOPE  BY  SIGLUM. 

Convert  degrees  - minutes  to  decimal 
degrees,  and  moves  the  fields  which  do  not 
need  to  be  changed  to  the  output  variables. 
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173/184.7 

185/206 


LJ 

L! 


206.5/229 


232/end 


Locates  predecessor  node  and  sets  ^ack  the 
pointer.  If  the  same  vertex  is  used  for 
all  four  seasons,  copy  it  to  all  four 
copies  of  OUT_ENVELOPE . Calculate 
edge  length,  and  if  latitude/longitude 
of  both  main  vertex  and  its  predecessors 
are  the  same  for  each  season  copy  it. 

If  the  vertex  differs  for  the  various  seasons, 
search  the  seasons  previously  calculated  to 
see  if  both  main  vertex  and  its  predecessors 
have  same  latitude/longitude,  in  which 
case  copy  edge  length,  otherwise, 
calculate  it  over.  Go  back  to  100  for 
next  vertex. 

Write  out  last  route. 


D.  PROGRAM:  BETA2 

REFERENCE  PARAMETERS  (SOURCE  PROGRAM) 

DSNAME  a CN5054 . PSl . BETA2 
VOL  * SER  * PUBL03 
UNIT  = 3330 

DCB  = (RECFM  = FB,  LRECL  = 80,  BLKSIZE  = 3120) 

Sequenced  73/80  WYLBUR  line  numbers 

REFERENCE  PARAMETERS  (OBJECT  MODULE) 

DSNAME  = CN5054  .PS1.LINKLIB  (BETA2) 

VOL  • SER  * PUBL03 
UNIT  = 3330 

Objective  Module  must  be  reprocessed  by  Linkage 
Editor/Loader  specifying  entry  point  PLISTART 

REFERENCE  PARAMETERS  (JCL)  (For  object  module  execution) 

DSNAME  = CN5054 .PSl . JCL (BETA2 ) 

VOL  = SER  = PUBL03 
UNIT  = 3330 
WYLBUR  EDIT  format 

FUNCTION: 

Sums  the  edge  lengths  from  the  Route  Envelope  Trees  (final  form)  to 
build  the  Port  Distance  Table,  determining  the  total  path  length  for 
each  Siglum.  It  must  be  run  once  for  each  season. 

CONTENTS 

Route  Envelope  Trees,  for 
appropriate  season 

Index  to  Envelop  by  Siglum 


INPUT  FILES 

DDNAME  DSNAME 

ENV  CN5054. PSl. OUTENV. SPRING 

SUMMER 

FALL 

WINTER 

ENVINDX  CN5054 .PS1.ENVINDX 
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OUTPUT  FILES: 

PRTDIST  CN5054.PS1.PRDIST. SPRING  Port  Distance  Table,  for  the 

|2£“  appropriate  seasons 
WINTER 

MAJOR  VARIABLES : 

ENVELOPE  Array  holding  needed  subfields  of  the 

Route  Envelope  Trees 

INDEX#ENVLP#BY#SGLM  Array  holding  the  index  to  Envelope  by 

LENGTH  Siglum  used  for  summing  up  the  edge  lengths 

PROGRAM  PARAMETERS: 

The  following  parameters  are  built  into  the  logical  structure 
of  the  program  and  if  changed  will  necessitate  program 
modifications: 


PARAMETER 

Number  of  routes  (84) 

Maximum  number  of 
Sigla  (506) 


Maximum  Number  of  Vertices 
in  Route  Envelope  Tree 
(900) 


REFERENCE  (LINE  NUMBER) 

56.2  I/O  Format  (skipped  over) 
12  Dimension 

56.3  I/O  Format 
100  DO  loop  limit 

8 Dimension 


PROGRAM  MESSAGES 


•ERROR,  SGLM  « XXXX' 


Response:  The  Index  to  Envelope  by  Siglum  containes  an 

•O'  instead  of  a line  number  in  the  Route 
Envelope  Trees- . Correct  the  Routing  Table  or  the 
Route  Envelope  Trees  (original  form)  as 
necessary  and  re-run  everything. 

PROGRAM  STRUCTURE: 


LINES 

51/56.3 

100/112 

108.4/111.6 


FUNCTION 
Read  in  Tables 
Sum  lengths 
Write  out  length 
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PROGRAM:  GAMMA1 


REFERENCE  PARAMETERS  (SOURCE  PROGRAM) 

DSNAME  - CN5054 . PSl. GAMMAl 
VOL  = SER  - PUBL03 
UNIT  = 3330 

DCB  * (RECFM  = FB,  BLKSIZE  = 3120,  LRECL  = 80) 

1; 

Sequenced  73/80  WYLBUR  line  numbers 

REFERENCE  PARAMETERS  (OBJECT  MODULE) 

DSNAME  * CN50 54.  PSl. LINKLIB ( GAMMAl) 

VOL  = SER  = PUBL03 
UNIT  = 3330 

Object  Module  - must  be  reprocessed  by  Linkage  Editor/  Loader 
specifying  entry  point  PLISTART 

REFERENCE  PARAMETERS  (JCL) 

DSNAME  = CN5054.  PSl. JCL (GAMMAl) 

VOL  * SER  = PUBL03 
UNIT  = 3330 
Wylbur  EDIT  format 

FUNCTION:  Constructs  the  Port  Hash  Table  from  the  Port  Table.  The 

Port  Name  Abbreviation  is  converted  to  a Hash  Index  by  taking  its 
S/370  numeric  value,  dividing  it  by  2999,  and  using  the  remainder. 

REVISION  STATUS:  Make  certain  JCL  specifies  most  recent  revision  of 

PORTTBL  before  running. 

INPUT  FILES: 

DDNAME  DSNAME 

PRTTBL  CN50S4 . PSl . PORTTBL 

OUTPUT  FILES: 

PRTHSH  CN5054 . PSl . PRTHA5H  Machine-readable 

Hash  table 


CONTENTS 
Port  Table 
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HASHOUT  CN5054 

MAJOR  VARIABLES: 
PORT#HASH#TBL 
PORT#  HASH  # TBL#  BUF 

PORT  # DSCRPTN 

PORT#NAME#A, 

PORT#NAME#B 

PORT#NAME#C 

NEXT#AVL#SPACE 

PROGRAM  STRUCTURE: 
Line 

40/47 

50 

51/68 

52/55 

56/66 


PS  1.  HASHOUT  Human-readable  Hash 

Table 


Port  Hash  Table 

Dummy  variable  overlayed  on  PORT#HASH#TBL 
during  output. 

Array  storing  the  needed  fields  from 
the  Port  Table 

Various  formats  of  the  temporary  varia- 
bles used  during  conversion  of  the  port 
name  abbreviation  to  the  Hash  Index 

Cursor  used  to  manage  loading  of  over- 
flow area  in  the  Hash  Table 


Format 

On-unit  to  enable  an  error  causing  Port 
to  be  printed  out. 

On-unit  to  end  the  program  after  all 
records  from  Port  Table  are  read. 

Loop  for  each  sucessive  port 

Render  Port  name,  and  convert  to  Hash  Index 

If  another  Port  already  has  that  Hash 
Index  determine  the  proper  bin  to  use  in 
the  overflow  indices. 

Store  the  port  name  and  number  in  the 
Port  Hash  Table. 


67/68 


75/83  Write  out  Port  Hash  Table 

PROGRAM  MESSAGES: 

Following  any  system  error  message  is  a print-out  of  Port  Table 
entry  being  processed  when  the  error  occurred,  the  Hash  Index  calcu- 
lated (if  any) , and  NEXT#AVL#SPACE  cursor. 
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PROGRAM:  GAMMA 2 


REFERENCE  PARAMETERS  (SOURCE  PROGRAM) 

DSNAME  = CN5054 .PS1.GAMMA2 
VOL  = SER=  PUBL03 
UNIT  = 3330 

DCB  = (RECFM  = FB,  BLKSIZE  = 3096  LRECL=7  2) 

Sequenced  73/80  W!'LBttr  line  numbers. 

REFERENCE  PARAMETERS  (OBJECT  MODULE) 

DSNAME  = CN5054 . PS1 . LINKLIB (GAMMA2) 

VOL  = SER  = PUBL03 
UNIT  = 3330 

(Object  module  must  be  reprocessed  by  the  linkage  editor/ 
loader  specifying  entry  point  PLISYART) . 

REFERENCE  PARAMETERS  (JCL  for  Object  Module  execution) 

DSNAME  = CN5054.PS1. JCL(GAMMA2) 

VOL  * SER  = PU3L03 
UNIT  = 3330 
Wylbur  EDIT  Format 

FUNCTION: 

Builds  the  Estimated  Travel  Tiroes  Table  by  tracing  the  Siglum 
Table  and  totaling  the  port  distances  from  the  Port  Distance  Table. 

It  assumes  a speed  of  12  knots.  For  ships  traveling  .from  one  canonical 
port  to  itself,  it  assumes  a travel  time  of  0 days  (which  may  later  be 
modified  by  Wylbur  EDIT  commands) . if  roust  be  run  4 tiroes,  once  for 
each  season.  Also  it  constructs  the  Multiple  Travel  Time  Tables  and 
places  the  additive  inverse  of  a pointer  to  them  in  the  proper 
place  in  the  Estimated  Travel  Tima  Table. 


FILES: 

INPUT: 

DDNAMB 

DSNAME 

CONTENTS 

SGLMTBL 

CN5Q54 .P31 .SGLKSBL 

Siglum  an  & Multiple  Routing 

T/'>lea 

PORTDST 

CH5054 . PS1* PRTDISY . SPRING 

P zt  Distance  Tables,  as 

SUMMER 

FALL 

appropriate  for  season. 

WINTER 
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OUTPUT*. 

TRAVTME  CN5054 . PS  1 . TRAVTME . SPRING  Travel  Time  Tables 

SUMMER 

' PALL 

WINTER 

SY SPRINT  Intermediate  checks, 

UPDATE  STATUS:  ' (CHECK) ' prefixes  and  line  75.5  used  during  the 
checkout  may  be  removed,  but  are  harmless  except  for  producing 
lots  of  output.  Otherwise,  the  program  is  finished. 

PROGRAM  PARAMETERS; 

Changes  in  the  following  parameters  will  require  program 
modification. 

PARAMETER 

Number  of  Cononical 
Ports  (69) 

Maximum  number  of 
Sigla  (506) 

Maximum  number  of  Multiple 
Travel  Time  Tables 
generated  (100) 

Estimated  Inverse  Speed 
(Days/Degree  Latitude) 


Siglum  Table 
Multiple  Routing  Table 
Travel  Time  Table 
Multiple  Travel  Time  Table 
Port  Distance  Table 
Origin  and  Dest.  Ports  numbers 
Intermediate  Port  numbers 
Accumulators  for  travel  time 
Multiple  Routing  Table  branch 
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MAJOR  VARIABLES: 
SGLMtTBL 
MULT#RTNO#TBL 
TRAV  # TME  # TBL 
MULT#TRAV#TME#TBL 
PORT#DIST 
PORT# 1, PORT# 2 
PORT,  PORTM 
TIME,  TIMEM 
BRANCH 


REFERENCE  (LINE  NUMBER) 
15,20  Dimension 
52,100  I/O  Format 
58  DO  Loop  limit 
21.5  Dimension 
56  I/O  Format 
21  Dimension 
100  I/O  Format 

64,68,80,85,91  Multiplier 


PROGRAM  STRUCTURE 


LINE 

50/56 

58/97.5 

59/60 

61/73 


76/78.5 

79/96 


FUNCTION 

Reads  in  the  Table 

Loops  through  the  Siglum  Table  using  all 
possible  port  combinations 
Initializes  for  each  origin  and  destina- 
tion port  pair 

Traces  the  Siglum  Table  recursively 
until  the  Multiple  Routing  Table  is 
encountered 

Initializes  the  Multiple  Routing  Table, 
trace  for  each  of  the  two  branches 
Traces  the  Siglum  Table  recursively, 
choosing  the  proper  multiple  routing 
table,  branch  for  each  time  through. 
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G . PROGRAM:  DELTA 1 

REFERENCE  PARAMETERS  - SOURCE  PROGRAM; 

DSNAME  = CN5054 . PS2 . DELTA1 
VOL  = SER  = PUBL03 
UNIT  = 3330-1 

DCB  = (RECFM  = FB,BLKSIZE  = 3120,  LRECL  = 80) 

REFERENCE  PARAMETERS  - JCL 

DSNAME  = CN5054 .PS2.DELTA1. JCL 
VOL  = SER  = PUBL03 
UNIT  = 3330-1 
WYLBUR  EDIT  format 

FUNCTION:  Provides  a visual  check  on  the  Route  Trees 

to  catch  transcription  errors.  It  plots  each 
vertex  of  the  given  route  and  designates  canonical 
ports  by  a C and  a non  canonical  port  by  a V.  The 
great  circle  path  is  generated  and  displayed  between 
each  pair  of  vertices.  A 5°  allowance  is  allotted 
on  either  ends  of  the  route  for  clarity.  The  program 
was  used  for  each  of  the  four  seasons  on  all  85  routes 


INPUT  DATA  SETS: 
DSNAME 

CN5034.PS2.LANDSEA 
CN5054 . PS1. SGLMNDX 
CN5054 . PS1 . ENVINDX 


CONTENTS 

World  in  terms  of  periods  and  blanks 
Inverse  Siglum  Table  and  Siglum  Map 
Index  to  Envelope  by  Route  and 
Index  to  Envelope  by  Siglum 
Final  version  of  Route 
Envelope,  one  per  season 


CN5054.PS1.0UTENV. SPRING 
CN50  54 . PS1 . OUTENV .Fall 
CN5054.PS1. OUTENV .WINTER 
CN5054.PS1. OUTENV . SUMMER 
MAJOR  VARIABLES: 

Contains  Landsea  table 
Contains  Siglum  indices 
Index  Envelope  by  Siglum 


LSEA4 . WORLD 1 

SIGMP2 

INSIG 


A-24 


PROGRAM  STRUCTURE 
Lines 
1/17 

18/37 

38/70 

71/115 

115.1/281 


Function 

Dimensions  various  arrays,  i.e., 

Route  Trees,  Landsea,  Siglum  Map, 

Index  Envelope  by  Siglum 

Reads  in  all  the  necessary  files. 

A factor  of  90.5  is  added  to  the 
latitudes  and  181  to  the  longitudes 
Checks  if  we  are  beginning  a new  route, 
and  if  so,  writes  out  the  previous  route. 
Determines  the  minimum  and  maximum  longi- 
tudes 

Determines  the  great  circle  path  between 
pairs  of  vertices  according  to  the  vertex 
widths . 

Stores  the  corresponding  code  for 
canonical  ports,  non  canonical  ports,  and 
great  circle  paths. 


H.  PROGRAM:  DELTA2 

REFERENCE  PARAMETERS  - (SOURCE  PROGRAM) 

DSNAME  = CN 5 0 5 4 . PS  2 . DELTA2 
VOL  = SER  = PUBL03 
UNIT  = 3330 

DCB  = • (RECFM  * FB  ,LRECL  * 80,BLKSIZE  = 3120) 

REFERENCE  PARAMETERS  - JCL 

DSNAME  = CN5054.PS2.DELTA2. JCL 
VOL  = SER  = PUBL03 
UNIT  = 3330 
WYLBUR  EDIT  FORMAT 

FUNCTION:  Checks  the  canonical  pert  mapping.  The  noncanonical 

ports  are  plotted  using  symbols  to  denote  the  various 
associated  canonical  ports.  This  provides  a quick 
checking  procedure  to  determine  if  in  a given  area 
there  are  more  than  one  canonical  port  associated 
to  it.  It  also  spots  typegraphical  error  occurring 
in  the  Port  Table.  An  example  of  this  is  the 
occurrence  of  a port  in  the  middle  of  a large  landmass. 
The  program  is  run  in  two  parts,  once  for  canonical 
ports  1-34  and  is  repeated  for  canonical  ports  35-69. 


INPUT  FILES: 

DSNAME 

CN5034 .PS2.LANDSEA 


CN5054.PS4. PORTTBL 


MAJOR  VALUABLES: 
SUMTAB 


WORLD 


CONTENTS 

The  world  stored  in  terms  of  periods 
and  blanks 

Table  of  all  the  ports  and  the  corres- 
ponding latitude,  longitude  and 
associated  canonical  port. 


Vector  containing  the  code  for  the 
canonical  ports. 

Array  containing  the  landsea  table* 
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PROGRAM  STRUCTURE 
LINES 
1/14 
15/20 


21/32 


23/53 


FUNCTION 

Dimensioning  various  arrays 
Read  in  the  port  table  with  the  correspon- 
ding latitude,  longitude  and  associated 
canonical  port. 

Determines  if  the  assoicated  canonical 
port  lies  in  the  range  of  this  run.  Then 
it  deals  with  the  North,  South,  East, 

West  by  adding  or  subtracting  the  appro- 
priate value. 

Checks  if  two  different  canonical  port 
codes  lie  in  the  same  1°  x 1°  square,  and 
if  so,  designates  it  on  the  output  by 
an  asterick.  It  also  displays  the  full 
world . 


I . PROGRAM;  DELTA3  . 

REFERENCE  PARAMETERS:  SOURCE  PROGRAM; 

DSNAME  = CN5034 . PS2 . DELTA4 
VOL  = SER  = VOL  » PUBL03 
UNIT  = 3330 

DCB  = (RECFM  = FB#  BLKSIZE  = 3120,  LRECL  = 80) 

REFERENCE  PARAMETERS  - JCL 

DSNAME  = CN5Q34 .PS2.DETLA4 . JCL 
VOL  = SER  = VOL  = PUBLG3 
UNIT  = 3330 

DCB  = (RECFM  - FB,  BLKSIZE  = 3120,  LRECL  * 80) 

FUNCTION:  Checks  the  routing  table.  Prevents  the  occurrence 

of  infinite  loops,  and  does  various  checks  on  non- 
direct  routes.  In  the  case  of  a multiple  routes  the 
check  was  done  manually.  When  an  error  is  spotted  a 
corresponding  message  is  written. 


INPUT  DATA  SETS: 

DSNAME  CONTENTS 

CN5054 .PS1 .RTNGTB L. FINAL  69x69  array  called  Routing  Table. 
MAJOR  VARIABLES 

ROUTEM  Alphanumeric  array  containing  the  code  ' ' 

for  direct  route,  'P'  for  intermediate 
port  and  'B'  for  multiple  route 

ROUTEM  Array  containing  the  various  routes 

between  2 ports 


PROGRAM  STRUCTURE 
LINE 
1/10 
11/19 

20/end 


FUNCTION 

Dimensions  various  arrays 

Reads  in  the  Routing  Table  and  declares 

some  logical  variables 

Reads  a pair  of  canonical  ports.  It  traces 
through  the  path  and  undergoes  various 
checks  along  the  way.  If  there  are 
inconsistencies  an  error  message  is 
written. 


J.  PROGRAM:  GAMMA 

•'  ■ ••  t . . 

REFERENCE  PARAMETERS: 

(OBJECT  MODULE) 

(Source  Programs) 

DSNAME=CN5054 . PSl . GAMMA 

DSNAME=CN5054 . PSl . LINKLIB(GAMMA)  C 

VOL=SER=PUBL03 

VOL=SER=PUBL03 

UNIT=3330-1 

UNIT«3330-1 

DCB= ( RECEM=FB , LRECL=8  0 , 

Must  be  reprocessed  by  loader 

BLKSIZE=3120) 

: ' 

Sequenced  1/834  WYLBUR  line 

numbers 

MACROS  INVOKED: 

DEFINE#MONTH#TBL 

DEFINE# PORT#HASH#TBL 

TEST 

SET#INVALID#FIELD#LIMIT 

COMPARE 

SEARCH#PORT#HASH#TBL 
SEARCH# MONTH# TBL 
GET# DATE 
SHPDYS 

MAJOR  VARIABLE  (other  than 

those  limited  to  use  in  macros) 

NEW# INPT#  SHP 

Overlayed  on  each  record  from 
SHPINPT  in  buffer 

INPT#SHP 

Contains  data  moved  from  NEW# INPT 
#SHP  (both  variables  are  needed  to 
maintain  synchronization  since  2 
sequential  input  records  must  be 
compared  to  check  for  duplicated 
records) . 

VALlFLGS , NEW#VAL#FLGS 

Validity  flags  to  ID  fields  in 
NEW#INPT#SHP,  INPT# SHP  respectively. 

TRNS#SHP 

Contains  data  obtained  by  translating 
INPT# SHP  by  converting  port  name 
abbreviations  to  port  numbers  and 
month  abbreviations/day  of  month  to 
day  of  year. 
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SHP(J) 


BETA 

W 

P 

GAMMA  (J) 

PHI (J) 

F(J) 

C(J) 

RHO 

SUM 

PROGRAM  STRUCTURE: 

LINES : 

238/261 

262/293 

294/330 


J=1  to  20 

Search  table  contains  15  to  20 
ships  from  OLDSDT. 

Index  of  first  ship  in  search  table 
Index  of  last  ship  in  search  table 
Points  at  last  ship  in  search  table 
matched . 

Points  to  ship  preceding  SHP(J)  in 
search  table.  GAMMA (BETA)  = 0 . 

Points  to  ship  succeeding  SHP(J)  in 
search  table.  PHI (W) =99, 

0 indicates  position  J empty  in  search 
table . 

1 not  empty 

0 indicates  SHP(J)  unmatched 

1 SHP  (O')  matched 

Counts  inserts  to  the  search  table. 

Counts  comparisons. 

FUNCTION 

Initialization  for  the  first  processed 
week  or  the  13th  week  of  the  season* 
Initialization  for  each  week  of  process- 
ing . 

Move  previous  NEW#INPT#SHP,NEW#VAL#FLGS, 
and  NEW#ERR#INT,  respectively.  Read  new 
NEW#INPT#SHP,  check  ID  fields  for  validity 
via  TEST  macro  and  simultaneously  set 
NEW#VAL#FLGS.  Invoke  COMPARE  macro  to 
test  NEW#INPT#SHP  and  INPT#SHP  for 
duplication.  Maintain  COUNTER. INPT# SHP, 
and  synchronize  TRANSMIT  and  ENDFILE 
conditions  with  INPTISHP  processing. 
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331/413 


332/352 

353/399 


400/413 


415/680 


415/426 


Branch  to  ERR#WPT#SHP  if  too  many  errors, 
or  TRANSMIT  condition  raised,  and  to 
FINISH  on  ENDFILE  condition. 

Translate  PEE  (Port  Entrance/Egress 
Data)  from  INPTtSHP  to  internal 
form  and  store  in  TRNS#SHP, 
simultanteously  setting  validity  flags. 
Translate  origin  port  data. 

Translate  port  of  last  call  data. 

If  port  name  is  missing,  but  other 
data  given,  then  ship  has  arrived  at 
its  final  destination,  so  use  port 
of  destination  name  and  flag  this 
field  in  TRNS#SHP  as  being  the 
final  destination.  Use  Last  Call 
Code  CAR',  ’IN',  'PD',  'SD')  to 
determine  whether  LCALL  date 
should  be  interpreted  as  arrival  or 
departure  date. 

Translate  port  of  destination  name, 
if  it  has  not  been  moved  to  LCALL 
field.  Port  name  abbreviations  must 
be  translated  to  DOT  port  number 
(modified) , and  mon th/day-of -month 
pairs  to  day-of-year. 

Sets  up  the  search  tables  for  all  the 
weeks.  Searches  and  matches  the  ships. 
Calls  MACRO  SHPDYS  which  actually 
calculates  the  days  between  ports  and 
increments  PRTRT . 

Minor  bookkeeping  to  make  the  two  parts 
of  the  program  more  compatible  and 
programming  somewhat  easier  as  values 
are  transferred  from  one  variable  to 
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! 426/431 

432/464 


465/507 


508/532 


533/597 

598/656 


657/662 


a second  variable  that  will  be  used 
in  this  part  of  the  program. 

Test  for  week  13,  that  is  the  13th 
week  of  the  season. 

Two  functions:  sets  up  the  initial 
search  file  for  the  13th  week  of  the 
season, does  the  initial  calculations 
for  any  ship  that  makes  its  first 
appearance  after  the  13th  week  of  the 
season. 

Loads  the  search  file  at  the  beginning 
of  each  week.  Reloads  the  search  file 
as  comparisons  are  made. 

Compares,  takes  a TRNSiSHP  and  searches 
the  file  for  a match.  If  a match  occurs, 
goes  to  calculations. 

Inserts  new  ships  into  search  file. 
Determines  if  the  data  on  a ship  fits 
into  one  of  six  categories  available. 

Then  sets  the  data  up  in  a form  compatible 
with  MACRO  SHPDYS.  SHPDYS  is  called  and 
actual  days  calculated  and 
the  proto- routing  table  incremented. 

Once  the  calculations  are  done,  test 
on  SUM.  If  SUM  is  still  less  than  10 
control  is  sent  back  to  GET#INPT#SHP 
If  SUM  is  greater  than  10,  test  to  see 
if  the  pointer  is  within  5 of  the  last 
ship  on  the  list.  If  pointer  is  not. 

Bend  control  to  GET#INPT#SHP.  If  pointer 
is  within  5 then  send  control  to 
WRrJTE#0UT, 

Writes  out  the  top  of  the  search  file. 
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Resets  the  chain  and  various  pointers 
and  sends  control  back  to  LOAD  where 
the  search  file  is  again  filled  with 
new  ships . 

681/733  ERR# INPT#SHP  program  segment  for 

processing  the  errors  not  correctable 
by  program  logic. 

734/767  ENDFILE.  Closes  ship  input  file,  writes 

out  error  messages,  on  file  ERRMSG 
and  SYSPRINT. 

768/781  Writes  out  the  last  of  the  search  file, 

makes  sure  all  ships  have  been  taken 
off  of  the  input  file  and  written  out 
on  the  new  file  for  next  week. 

782/798  If  INDX  is  equal  to  1,  sends  control 

back  to  INITIALIZE,  closes  the 
appropriate  file.  If  INDX  does  not 
equal  1,  dumps  proto-routing  table  for 
that  week  onto  a temporary  file.  If 
also  MOD (INDX, 4)  = 1,  sends  control  to 
DUMPIPRTRT. 

799/834  Dumps  proto-routing  table  onto  a file 

for  the  appropriate  week.  Sends  control 
back  to  initialize  it,  unless  it  is  the 
first  week  of  the  season,  i.e.,  the  last 
week  processed  and  then  closes  down  all 
the  files,  writes  out  a message  saying 
'normal’  and  ends  the  program. 

INTERNAL  TABLES 

PORTIUASHITBL  Used  for  translating  port  name  abbrevi- 

ation to  DOT  port  number.  Documentation 
elsewhere. 

CAN# PORT  MAPS  DOT  PORT  number  to  canonical  port 

number. 
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DAYS 

IDAYS 


Contains  travel  time  table  for  the 
season.  Documentation  elsewhere. 
Contains  multiple  travel  time  table 
for  the  season.  Documentation  else- 


PRTRT 


CONTROL  VARIABLES 
VALIFLGS . TO# ID 
OCCURENCES: 

NAME 

VALIFLGS. TO# ID 
NEW# VALIFLGS . TO# ID 


where . 

The  proto-routing  table.  A 3x69x69 
table  containing  the  cummulative  days 
spent  by  each  of  the  three  ship  types 
traveling  between  any  two  canonical 
ports . 


HOW  SET  APPLIED  TO 

By  TEST  macro  INPT#SHP.ID 

By  TEST  macro  NEW#INPT#SHP.I 

By  LOADITBLISHP 
macro  using 
VALIFLGS . TO# ID 


FUNCTION: 

BIT  1 

Name  Field 

•l'B  - Valid;  'O'B  - invalid 

BIT  2 

Not  used  - reserved 

for  TYPE  field 

BIT  3 

Flag  Field 

'l'B  - valid;  'O'B  - invalid 

BIT  4 

Year  built  field 

'l'B  - valid;  'O'B  - invalid 

BIT  5 

Type  and  weight 
fields 

'l'B  - both  valid;  'O’B  - at 
invalid 

SPECIAL 

FUNCTIONS  (VALIFLGS.  TO# ID , NEW# VAL# FLGS ) 

* 00000000'B 

ENDFILE  Condition 

'QOOQQQOO'B 

TRANSMIT  Condition 

'OOOQGQOO'B 

First  Record  (void) 

CONTROL 

VARIABLES 

VALIFLGS. TO# PEE 

FORMAT:  BIT (8) 

OCCURENCES: 
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NAME 

TBNStSH? . VALIFLGS . 

TRNSISHP .VALIFLGS . 
TOIPEE . PORTICALL 

TRNSISHP. VALIFLGS. 
TO#PEE.PRT#DST 


HOW  SET 

By  Search  # MONTH# TBL. 

GET# DATE,  and 
SEARCH#  PORT  # HASH#  TBL 

SEARCH# PORT# HASH# TBL 
macros  and  surrounding 
code 


APPLIED  TO 

TRNStSHP.PEE 

PRT#OR 

TRNStSHP.PEE 

PRT#CALL 

TRNStSHP.PEE 

PRT#DST 


FUNCTION: 
BITS  1&2 


Port  number  field 


' ll'B  - fully  valid 
• 10'B  - valid,  but  port  number 


i 

corresponds  to  a port  name 

t 

£ L 

flagged  as  a generic  term  in 

I 

5 

§ 

j l 

the  port  table 

*00'  - invalid  or  blank 

BITS  3&4  Arrival  data  field 

'll'  - fully  valid 

j . ' 

Departure  date 

'10'  - upper  limit,  obtained 

1 i 

field 

by  correlation  on  month  name 

BIT  7 


or  by  use  of  day-of-report, 

- when  month  name  invalid, 
or  actual  data  in  arrival 
field  Of  PRT#LCALL  if  LCALL 
code  was  invalid. 

*01*  - lower  limit,  obtained  by 
using  first  day  of  month  when 
day  of  month  invalid,  or  actual 
date  in  Departure  field  of 
PRT# LCALL  if  LCALL  field  waa 
invalid 

'00'  - filed  missing  or  totally 
invalid,  when  in  Departure 
field  or  Arrival  field  at 
PRT#DST  probably  means  that 
ship  has  not  yet  left  or  arrived, 

’l*  indicates  that  this  port 
was  designated  as  a destination 
port  in  input  data 

’0*  indicates  that  it  was  not  so 
designated 


CONTROL  VARIABLES 

NAME 

GATE 


TESTED  ON 

By  SECOND# PART 
to  control  loading  of 
the  search  table 


VALUES 

*0’  - have  fiist 
TRNSISHP  of  wrek 
Need  to  load  the 
search  table  for 
the  first  time 
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NAME 


KFLG 


TAG 


K 


INDX 


TESTED  ON 


By  SECOND# PART 
WEEK# 13  to  transfer 
control  after  initial 
calculations  on  ship 

By  SECOND# PART 
to  indicate  if  we  are 
searching  forward  or 
backward  through  the 
search  table 

By  SECOND# PART 
Transfer  of  control 
during  loading  of 
search  table 


By  SECOND# PART 
WRAP#UP 


VALUES 

•1'  - have  a TRNS#SHP 
initial  loading.  Has 
started 

'-l*  - need  to  get  a 
TRNS#SHP  after  refill- 
ing of  the  search 
table 

*0'  - we  are  in  week 
13  of  the  season 

•1'-  we  are  in  OTHER 
week  than  13 

'O'  - forward 
' 1 ' - backward 


'O'-  filling  the  tables 
for  the  first  time 
this  week 

'1'  - refilling  the 
table 

Counts  no.  of  weeks 
processed  this  season. 
Sends  control  to  DUMP 
#PRTRT  at  INDX  ? 1 
but  MOD  (INDX,4)=:1 


INPUT  DATA  SETS: 

FILENAME  = SHPINPT 

DDN'AME  = DOTINP52,  DOTINP51,  etc.  Determined  dynamically  during 
execution  by  concatenating  WK#OF#KPRT  to  "DOTINP'. 

CONTENTS:  DOT  Marine  Traffic  Generator  Input  Tape  for  appropriate 

week. 

ENVIRONMENT  Specification:  F RECSIZE  (84)  TOTAL  BUFFERS  (10) 

DCB  Specifications  (on  DD  card)  : TRTCH  = ET,  OPTCD  = C 

FILENAME,  DDNAME  = PRTHSH 
CONTENTS:  Port  Hash  Table 

Referenced  by  DEFINE# PORT# HASH# TBL  macro 
FILENAME,  DDNAME  = NEWPRT 
CONTENTS:  Canonical  Port  Numbers 

ENVIRONMENT  Specification:  FB 
RECSIZE  (200)  BLKSIZE  (5000) 
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FILENAME  - TRAV 


DDNAME  - TRAVTME . FALL 

TRAVTME. WINTER 
TRAVTME. SPRING 
TRAVTME . SUMMER 

CONTENTS:  Number  of  shipdays  between  any  two  canonical  ports.  Also 

multiple  route  table  for  voyages  that  pass  through  the  canals 

ENVIRONMENT  Specification:  FB 

RECSIZE  (80)  BLKSIZE  (3120) 

OUTPUT  DATA  SETS: 

FILENAME,  DDNAME  = ERRMSG 

Print  file  (SYSOUT=A)  containing  all  error  messages  generated. 
FILENAME, DDNAME  = SYSPRINT 

CONTENTS:  System  error  messages,  copies  of  certain  important  error 

messages  from  ERRMSG,  system  debugging  output  . 

Print  File  (SYSOUT=A) 

FILENAME~PROTQ 

DDNAME  = PROTO 4 8 
PROTO 4 4 
PROTO 40 

CONTENTS:  The  Proto- Routing  table  for  each  month. 

ENVIRONMENT  Specification:  FB 

RECSIZE  (72)  BLKSIZE  (3120) 

FILENAME  « PRO#TEM 
DDNAME  = PRO. TEMP 

CONTENTS:  Temporary  proto-routing  table  to  be  used  in  case  a restart 

is  necessary.  At  the  end  of  normal  execution  contains  the  cumulative 
list  for  ports  of  origin. 

ENVIRONMENT  Specification:  FB 

RECSIZE  (80)  BLKSIZE  (3120) 

MACRO  FUNCTIONS: 

TEST  TEST  Checks  ID  field  of  input  record 

and  sets  corresponding  validity  flags. 
COMPARE  COMPARE  Compares  ID  fields  of  two  records  and 

checks  for  identity. 


SIFL 


SET# INVALID#FIELD#LIMIT 
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Used  with  TEST,  DPHT,  SMT,  AND  GD 
macros  to  limit  total  number  of 
invalid  fields  to  be  allowed  before 
branching  to  error  routine. 


DPHT 

DEFINE#PORT#HASH#TBL 

Reads  in  Port  Hash  Table  and  contains 
declarations  for  associated  variables. 

SPHT 

SEARCH#PORT#HASH#TBL 

Searches  Port  Hash  Table  for  a given 
port  name  abbreviation. 

DMT 

DEFINE#MONTH#TBL 

Initializes  Month  Table  to  optimum 
sequence  for  linear  searching. 

SMT 

SEARCH#MONTH#TBL 

Obtains  starting  day  of  given  month 
from  the  Month  Table. 

GD 

GET# DATE 

Adds  day  of  month  to  starting  day  of 
month  to  obtain  day  of  year. 

SHPDYS 

Calculates  ship  days  between  canonical 
ports.  If  canonical  ports  are  identical 
assigns  max.  value  of  2.  Increments 
Proto-Routing  table. 

SYSTEM 

OPERATIONS 

MACRO 

DEFINITIONS: 

REFERENCE  PARAMETERS: 

DSN  = CN5054.PS1.PLI 

.MACLIB 

VOL  = SER  = PUBL03 
UNIT  = 3330-1 

DCB  = (DSORG=PO ,RECEM=FB ,BLKSIZE=3120 ,LRECL=80) 

Sequenced  73/80  WYLBUR  line  numbers 

ACCESS  METHOD: 

1.  Include 

//PL1.SYSLIB  DD  DSN  * CN5054 .PS1.PL.MACLIB 
//  VOL  = SER=PUBL03,UNIT=3330-1,DISP«SHR 

card  among  compiler  DD  cards. 
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2.  Include 

%INCLUDE  (BNTS) , (DST) , . . . (GTS) ; card  specifying 
short  form  of  all  macro  names  to  be  used  in  the 
source  text  for  each  external  procedure. 

I.  If  Gamma  is  to  be  modified. 


f 

e 


< 

i I 

v 1 


1.  Make  changes  in  source  code. 

DSN  = CN5054.PS1.  GAMMA 

VOL  = SER=PUBL03 

2.  RECOMPILE  and  Linkedit.  Gamma/  JCL  to  accomplish  this  is 
located  in  Data  set: 

DSN  = CN5  05  4.  PS  1.  GAMMA.  COMPLIN 
VOL  = SER=PUBL03 

3.  The  load  module  is  placed  in  a library. 

DSN  = CN5Q54.PS1.LINKLIB (GAMMA) 

VOL  = SER“PUBL03 

4.  To  execute  the  load  module  Gamma  the  JCL  is  located  in 
GAMMA. JCL. EXEC. 

DSN  = CN5054.PS1. GAMMA. JCL. EXEC 
VOL  * SER-PUBL03 

The  JCL  for  the  data  set  name  GAMMA. COMPLIN  is  as  follows: 


^ EXEC  PL10PCL  sPftRM .PL I=< 'CHARSET <48> >NOINSQURCE jDPT<TIME> ' > 
" 'NEST ^STORAGE? OFFSET ^AGGREGATE'- > t 

PRRM.LKEB='LIST >MAP'  • ' 

•■'•••'PL I .SYS IN  BB  DSN=C(f5054.P$l  .GAMMA >UNIT»3330~1  f 
'V  VQL=SER=PUBL03 » DISP=<SHR  >KEEP) 

s/PLl  .SYSLIB  BD  BSN=CN50S4.PS1  .PLI  .MFlCLIB >VDL®SER=PUEL03 ? 

UN I T=333 0- 1 I SP=  < SHR  > KEEP  > 

•^LKEB.SYSLMDB  BD  DSN-CN5054.PS!  .LINXLIB  .*UNIT=3330-i » 

" _VnL=SER=PUBLp3 >BISP=<M0B >KEEP) 
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The  JCL  for  GAMMA. JCL. EXEC 


£.  a-JOBLIB  BB  DSN=CN5054.  PS-1 . LINKL  IB*  UN  IT-3330- 1? 

3.  A'  VDL=SEP=sPUBL03* BISP^SHR 

4.  //GO  EXEC  P6M*6fiMMA* PhRM='  ISP  <1  02K)  /D" 

5.  -aSYSPRINT  DD  SYSDUT=*  > 

6.  //EVEN  BB  BSN=CN5 054 . PS  1 . TEMP . PORT . SHP .FILE1*  UN I T=333  0- 1 * 

7.  -a-  VDL=SEP=PUBL03* BISP® <QLB*  KEEP>  , BCB=  <RECFM=FB»  LRECL=38> 

8.  ABODES  DB  BSN«CN5054.PS1. TEMP. FORT. SHP. FILE2»UNIT=3330-1 * 

9.  -A’  VOL=S'ER=PUBL 03 * DISP®  <DLD*KEEP>  * BCB=  <RECFM=FB»  LRECL=38> 

10.  //NEWPPT  BB  BSN=CN5054. PS1.NEWPRT * 

1 1 . //  VQL=SEP*PUBL03* BISP=  CSHR* KEEP)  *UN IT =333 0-1  * 

12.  A BC  B=  <REC  FM=F  B * LREOL=£  0 0 * ELKS I ZE=5  0 0 O') 

1?.  //EPRMSG  BB  SYSaUT=fi*BCB=<BLKSIZE=£64G*LRECL=132* 

14.  //  RECFM=FB> 


15. 

16. 
17. 
IS. 
19. 
£0. 
£1. 
L-  L—  . 
i"  :• . 

c',4. 
£5. 
£6. 
£7 . 
£8. 

29. 

30. 

31. 


34. 

35. 

36. 


40. 

41. 
4c'  • 

43. 

44. 
15. 
4b. 

47. 

48. 

49. 

50. 


■j  ;• . 

54. 

55. 
5b. 


//PRTHASH  BB  BSN=CN5054.  PS  1 . PPTHSH*  VDL®SER=PUBL03* 

//  BISP®SHR»  UNIT =333  0- 1 

aTRRV  BB  D5N=CN505*. PSl . TRftVTME . WINTER * V0L=SER=PUBLQ3* 

//  UNI  T=3330  -1 > D ISP=  <DLD«  KEEP.; 

//BQTINP01  BB  VOL=S ER=PS 1 0 0 1 * UNIT=3330* DSN*CN5054. PSl . COPYOl > 

.A-  BCB=  <RECFM=F  B > LRECL-84 , ELKS  I ZE=3 1 08>.  > B I SP=SHR 

■"•-'  DDT  I NP  0£  BB  VGL=SEP-PS  1 0 0 1 * UN  I T=?33  0 > DSN=CN5  054.  PS  1 . CDPY  0£  * 

//  BCB= <RECFM«F B * LRECL-84 * BLK S I ZE=3 1 08> , D I SP=SHR 

//DOTINP03  BB  VOL=SER=PSI 001  *UNIT=?330* DSN=CN5054. PS1 . COPY 03 * 

A"  DOBS’  <PECFM=FB* LPECL=84 * BLKSIZE-31 08)  > DISP=SHR 

// DOT I NP 04  B D VOL=$  ER=PS I 00 1 * UN  1 T-33  3 0 - DSN=CN5 054 . PS 1 . COPY  04 * 

//  DCB= <RECFM®F B * LREC L=84 . BLKS I ZE*3 1 08> » D I SP-SHR 

//DOT I NP  05  BB  VOL=SER-  PS  1 0 0 1 * UN  1 T=333  0 * DSN=CN5  054 . PS  1 . COPY  05  * 

//  DCB® <RECFM=FB» LRECL«84* BLKSIZE-31 08)  * DISP=SHR 

//DOT  IMP  06  DD  VQL=SER=PSI  001 * UN  I T =333  0 * BSN=CN5  054 . F S 1 . COPY 06* 

//  DCB®  <RECFM*FB»  LRECL=84 * ELK'S  I ZE=3 1 08>  > B I SP-SHR 

//DOT I NP 07  DB  VOL=SEP=PS 1 001 * UNI T-3330* DSN-CN5054. PS1 . COPY 07* 

//  DCB=<RECFM=FB*LRECL=84*BLKSIZE=3ie8>*niSP®SHR 

//DOT  I NP 08  DB  VOL=SER=PSI 001  * UNIT-3330*  B$N=CN5054. PSl . CDPY 08  * 

//  DCB=cRECFM=FB*LRECL®84»£LKSIZE=3108>*BISP=SHP 

//DOT I NP 09  BB  VOL=SER®PSI 001 * UNIT®3330* DSN-CN5054. PSl . CDPY 09 * 

//  DCB=<RECFM=FB*LRECL=84* BLKSIZE=3108>  5 DISP-SHR 

//DQTINP1 0 BB  VOL=SER=sPSI  001 * UNIT-3330*  DSN=CN5054.  PSl . COPY1 0* 

//  DCB®  <RECFM=FB»LRECL=34*BLKSIZE=31 08>  * DISP=SHR 

//DDT I NP 1 1 DB  VDL=SER=PS 1 0 0 1 * UN I T®333  0 » DSN-CN5 054 . PS 1 . CDPY 1 1 , 

//  DCB=  <.RECFM«FB»  LRECL=84*  BLKS  IZE=31  08>  * DISP-SHR 

//DDT  I NP  1 £ DD  VDL-SEP-PS  1 0 0 1 * UN  I T=333  0 * DSN*CN5  054  .PSl.  CDPY  1 £ * 

//  DCB*  <RECFM-FB*  LRECL-84*  BLKSIZE=31  OS)  * DISP-SHR 

//DOT I NP 1 3 DB  VOL=SER®PS 1 0 0 1 * UN I T=333 0 * DSN-CN5 054 .PSl. CDPY 1 3 * 

//  DCB®  <RECFM®FBs  LRECL-84*  BLKSIZE®31 08)  * DISP-SHR 

■'/PROTO  0 1 DB  DSN=CN5  054  .PSl.  PPOTD  0 i > UN  IT-333  0- 1 > 

//  DISP= (NEW*  KEEP>  » DCB= tRECFN-FB ? LRECL=80> 

BLKSIZE=35£0> » VQL»SER»PUBL03« SPACE® <TRK» <10* 5> > 

••  /PRDTO  05  BB  DSN==CN5  054 . PS  1 . PROT  D 05*  UN  1 1=3330-1  * 

D ISP= (NEW*  KEEP* * DCB® <RECFN=FB* LPECL-80* 

//  BLKS  I ZE=35£0>  * VOL=S£P==PUIl  03*  SPACE-  <TRK*  <10*5)  ) 

//PRDTO 09  DD  DSN=CN5054. DS1 . PROTO 09*  UN  I T»3330-l * 

//  DISP= (NEW* KEEP) * DCB* <RECFM*FB»LRECL*80» 

BLKS  I ZE=3530'1  » VOL=SER=PUBL  03?  SPACE®  <TRK»  <10*5>> 
//PRD“T EN  DD  VC1L=SER=PUBL03> UNI T=3330-  1 * 

BSN=CN5054. PSl . PRO. TEMP*  BISP- «DLB*  KEEP> * 

BCE’s  <RECFM=FB » LF:ECL®8 0 * BLKS  I ZE  *3 1 2 0> 
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An  execution  time  parameter  is  passed  to  Gamma  on  the  PARM  card. 
Code  is  as  follows: 

A — * Fall?  sets  WK#OF#R?RT«53 

Fudge=365 

B — ► Summer?  sets  WK#OF#RPRT=40 

Fudge=286 

C — ► Spring?  sets  WK#OF#RPRT=27 

Fudge=195 

D — «►  Winter?  sets  WK#0F#RPRT=14 


Fudge=104 

I E — b»  Test  case;WK#OF#RPRT=53 

L„. 


OT'PfiTal 


K.  PROGRAM:  EPSILON 

REFERENCE  PARAMETERS  (SOURCE  PROGRAM) 

DSNAME  = CN5054.PS1. EPSILON 
VOL  = SER  * PUBL03 
UNIT  = 3330-1 

DCB  = (RECFM  * FB,  LRECL  =80,  BLKSIZE  = 3120) 

FUNCTION:  EPSILON  reads  a monthly  Proto-Routing  Table  and  from 
it  generates  a SHP.DAYS  Table.  The  Proto-Routing 
Table  is  indexed  by  ship  type,  egress  port,  and  entrance 
port.  EPSILON  examines  each  element  of  the  Proto- 
Routing  Table.  If  for  each  ship  type,  t,  the  (t,A,B)-th 
element  of  the  table  is  zero,  then  no  action  is  taken. 

If  for  some  t the  (t,A,B)-th  element  is  non-zero, 
then  the  (A,B)-th  element  of  the  Routing  Table  is 
examined. 

i.  If  the  (A,B)-th  element  is  a blank  followed  by  an  integer, 
n then  the  route  from  A to  B is  a canonical  route  and 
the  contents  of  the  (t,A,B)~th  position  in  the  Proto- 
Routing  Table  are  added  to  the  (t,n)-th  element  of  the 
SHP.DAYS  Table  where  n is  the  (A,B)-th  element  of  the 
Siglum  Table.  This  is  done  for  each  t where  the  (t,A,B)-th 
elements  of  the  Proto-Routing  Table  is  non-zero, 

ii.  If  the  (A,B) -th  element  of  the  Routing  Table  is  N 

followed  by  a zero  then  A=B  and  the  contents  of  the 
(t,A,B)-th  element  of  the  Proto-Routing  Table  are  moved 
into  the  (t,A)-th  element  of  the  RESIDUAL  TABLE  for  random 
distribution  in  the  vicinity  of  port  A. 

iii.  If  the  (A,B)-th  element  of  the  Routing  Table  is  a P 

followed  by  an  integer,  M,  then  the  RECURSIVE  INDEX  LIST 
is  started.  P indicates  that  the  voyage  from  A to  B is 
a compound  route.  That  is,  the  route  is  not  direct.  M is 
the  canonical  port  number  of  the  end  of  the  first  segment 
of  the  route.  The  (A,M)-th  element  of  the  Siglum  Table 
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stored  in  the  RECURSIVE  INDEX  LIST.  Then  the  (M,B)-th 
element  of  the  Routing  Table  is  examined.  If  it  is  a 
canonical  route  number,  that  is  blank  followed  by  an  integer, 
then  the  (M,B)-th  element  of  the  Siglum  Table  is  stored 
in  the  RECURSIVE  INDEX  LIST,  and  the  list  is  terminated. 

If  the  (M,B)-th  element  of  the  Routing  Table  is  a B 
followed  by  an  integer  K,  then  the  MULTIPLE  ROUTING  TABLE 
is  used  to  determine  which  of  the  two  pointers  is  correct. 

If  the  (M,B) -th  element  of  the  ROUTING  TABLE  is  a P 
followed  by  an  integer  J,  then  this  step  is  repeated 
with  the  (M,J)-th  element  of  the  Siglum  Table  being 
added  to  the  RECURSIVE  INDEX  LIST,  and  the  recursion 
continues  until  the  list  is  terminated.  When  termination 
occurs,  then  the  RECURSIVE  INDEX  LIST  is  used  to  extract 
distances  from  the  PORT  DISTANCE  TABLE.  The  extracted 
distances  are  used  to  proportion  the  ship  days  among  the 
various  component  routes.  That  is  one  assumes  the  number 
of  days  spent  going  from  point  A to  point  M is  that  same 
fraction  of  the  total  days  spent  going  from  A to  B as 
the  distance  from  A to  M is  to  the  total  distance  from 
A to  B. 

If  the  (A,B)-th  element  of  the  ROUTING  TABLE  is  a B 
followed  by  an  integer  L,  then  the  various  ship  types 
influence  the  choice  of  routes.  One  examines  the  L-th 
row  of  the  MULTIPLE  ROUTING  TABLE  if  t=3,  i.e.,  ship  type 
is  large  tanker,  the  first  entry  of  row  L indicates  the 
next  canonical  port  for  the  route.  If  t*3  then  the 
second  entry  of  the  L-th  row  indicates  the  next  canonical 
port. 

After  processing  the  entire  PROTO- ROUTING  TABLE,  the  program 
writes  out  the  SHP.DAYS  TABLE,  zeroes  it,  and  proceeds  to 
the  table  for  the  next  month. 
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INPUT  DATA  SETS 


Filename  = FT33F001 

DSNAME  = CN5  0 5 4 , PS 1 . RTNGTBL . NEW 

CONTENTS:  Routing  Table 

Filename  = FT34F001 

DSNAME  = CN5054 .PS1. SGLMTBL 

CONTENTS:  Siglum  Table 

Filename  = FT36F001 

DSNAME  = CN5054.PS1.PRTDIST.FALL  (SUMMER, SPRING, WINTER) 
CONTENTS:  Port  Distance  Table 


Filename 
DSNAME 
CONTENTS : 


FT08F001 

CN5054 .PS1 .PR0T048 ,35,22,09 

First  Proto-Routing  Table  of  the  Season 


Filename 
DSNAME 
CONTENTS : 


FT09F001 

CN5054 .PSl .PR0T044 ,31,18,05 

Second  Proto-Routing  Table  of  the  Season 


Filename  = 
DSNAME 
CONTENTS  - 

OUTPUT  DATA 

Filename  = 
DSNAME 
CONTENTS : 

Filename  = 
DSNAME  = 
CONTENTS : 

Filename  « 
DSNAME 
CONTENTS : 


FT10F001 

CN5054 . PSI . PR0TO40 , 27,14,01 

Third  Proto-Routing  Table  of  the  Season 

SETS 

FT20F001 

SHP.  DAYS 48, 35, 22, 09 

Ship  Days  Table  for  the  first  month  of  the  season 
FT21FC01 

SHP .DAYS  44 , 31,18,05 

Ship  Days  Table  for  the  second  month  of  the  season 
FT22F001 

SHP. DAYS  40, 27 ,14,01 

Ship  Days  Table  for  the  third  month  of  the  season 


MAJOR  VARIABLES 


jo 

0 


RTA(69,69)  Alphabetic  information  from  the  Routing  Table 

RTN (69,69)  Numeric  information  from  the  Routing  Table 

PRT (3/69,69)  Proto- Routing  Table 

MRT(69 ,69)  Multiple  Routing  Table 

PDT (500)  Port  Distance  Table 

SIG (69 ,69)  Siglum  Table 

RESIdO , 69)  Residual  Table  for  random  distribution  in  the 
vicinity  of  ports 

RESINL (20)  The  RECURSIVE  INDEX  LIST,  to  temporarily  hold  Siglum 
values  for  compound  routes. 

SDT (3,500)  SHP . DAYS  Table  indexed  by  ship  type  and  Siglum 


RESINL (20) 


SDT (3,500)  SHP. 
PROGRAM  STRUCTURE 
29/42 


i 

43/44 

IS: 

!l 

46/92 

1 [! 

fc  w... , 

1 

49/51 

t 1 ; 

53/91 

L; 

57/58 

1 Li 

60/61 

FUNCTION 

Reads  in  the  Routing  Table,  the  Multiple 
Routing  Tables,  the  Siglum  Table,  the 
Port  Distance  Table  for  the  appropriate 
season  • 

Reads  in  the  Proto-Routing  Table  for  each 
month  of  the  season  • 

Performs  the  function  of  EPSILON  • 

Checks  if  Proto-Routing  Table  is  zero  for 
these  two  ports. 

Performs  the  correct  calculation  based  on 
the  alphabetic  value  of  the  Routing  Table 

Routing  Table  has  a blank  in  the  alpha- 
betic position  • 

Routing  Table  has  an  N in  the  alphabetic 
position  . 
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63/65 

66 

67/83 

87/90 

93/95 

96/102 

103/105 

106 


Routing  Table  has  B in  the  alphabetic 
position. 

Routing  Table  has  a P in  the  alphabetic 
position 

Accumulates  Siglum  in  the  RECURSIVE  INDEX 
LIST  for  cases  P and  B. 

Calculates  Ship  Days  Table  for  cases  P and  B. 
Writes  the  Ship  Days  Table  and  the  Residual 
Table  onto  the  disk  pack. 

Zeros  the  SHIP. DAYS  and  Residual  Tables 
for  the  next  month. 

Increments  counters  for  controlling  the 
number  of  months  processed  and  the  logical 
devices  for  I/O. 

Tests  on  the  number  of  months  processed. 


SYSTEM  OPERATIONS 


The  JCL  is  located  in  the  data  sets  JCL. EPSILON.  Below  is  a copy 
of  the  JCL. 


2.  //  EXEC  FTG1CG 

3.  //FGRT.SYSIN  PP  PSHsCNSO^.PS  1 . EPSILON ? I'M T*333 0-1  ? 

4.  //  VGL=SEP=PL'PL03?PISP=<SHP*KEEP> 

5.  //CG.FT33F001  PP  PSN=CH5054 .PS1 .PTN6TBL .NEW • 

6.  //  UMIT*3330-1  •VCL=SEP*PL‘BL03?PISP«<SHPfKEEP> 

7.  //(5G.FT34F001  PP  PSN=CN5054 .PS  1 .S6LMTBL *UNIT*3330-1 r 

S.  //  VGL=?EP=Pl'EL03.PISP=<SHP>KEEP> 

3.  /''6Q.FT36F001  PP  PSN=CN5054 .PS1  .PRTPIST .FPLL  rUNIT»333P-l  ? 

10.  //  VCL=SEP=P«.‘BL03*BISP=<SHP?KEEP> 

13.  //GO.FT08FOO1  PP  PSN«CN505a  .PS1  .PPGTDS  ?UNIT=3330-1  t 

13.  //  VOL=SEP.sPUBL03fPISP=<OLPfKEEP>  fPCP=<PECFM»FB» 

14.  //  LP£CL=9  0 • ELKS 1 ZE =352  0 > 

14.1  //GQ.FT09F001  PP  PSN«CNS054 .PS1 .PPGTOf fUNIT=3330-l * 

14.2  /V  VGL=SEP=PUBL03*PISP=<GLP *KEEP>  »PCB=(P£CFf*=FB ? 

14.3  //  LPECL=8P  ?BLKS I2E=3520> 

14. a //SQ.FT 1 OF 001  PP  PSN«CN5054.PS1 .PRQTOf ♦UNIT*3330-l r 


I 14.5  //  VOL=SEP=PUBL 03 • P I SP®  < DLP ? KEEP>  * PCB=  < PECFM=FP • 

| [ 14.8  //  LPEC  L=8  0 • BLK  S 1 2E =353  0 > 

f ] 15.  //GG.FT20F001  PP  PSN=CN5054.PS1 . SHP .BPYS48 rUNIT*3330-l t 

l * 16.  //  VOL=SEP»PUBL03*PISPs<;DLP.KEEP>?PCB=<PECFN=FBfLPECL=7£* 

f , 17.  //  BLKSI2E=3120>>SPPCE«<TPK.<10?5>) 

f 18.  //6G.FT21F0C1  PP  PSN=CN5054.PSl.SHP.BPYS44.UrUT*3330~l ? 

f 1 • 19.  " VOL = SEP  =PUPL  03  * P I S P=  < DLP  ? KEEP > ? PC B53 < P£CFN®FB » LPECLa?2 r 

| 20.  //  BLKSI2E=31c,0>‘SPPCE=iTPt'X10*5>> 

||  21.  //<?O.FT22FOOl  PP  PSN*CN505d  .PS1  .SHP  .PPYS40  ?L'NITB53330-1  • 

1 . 22.  //  VOL  =SEP=PL'BL 03 » P I SP»  < OLP *KEEP > •PCBSS < PECFH*FB »LPECL*72 1 


L.  PROGRAM;  THSTA1 
REFERENCE  PARAMETERS  - SOURCE  PROGRAM 
DSNAME  = CN50  54 . PS  2 . SPLINE2 
VOL  = SER  * PUBL03 
Unit  » 3330 

DCS  » (RECFM  = FB,  BLKSIZE  = 3120,LRECL  « 80) 

REFERENCE  PARAMETERS  - JCL 

DSNAME  = CN5054.PS2. SPLINE. JCL 
VOL  = SER  * PUBL03 
Unit  = 3330 
WYLBUR  Edit  Format 

FUNCTION;  Generates  a set  of  X values  to  be  used  in  Theta 

for  the  distribution.  It  takes  a set  of  data  points 
and  determines  the  Type  I*  Cubic  Spline  in  which  »t 
interpolates.  The;;  it  uses  two  iterations  of  New  ton ‘s 
method  to  determine  intermediate  values  between  any 
two  successive  sextile  points.  Thus  producing  a 
set  of  13  X values.  This  is  done  for  each  season. 
INPUT  FILES; 

DSNAME  CONTENTS 

CN5054.PS2.SEXTILB.SP  A table  of  sextile  points 
CN5054.PS2.SEXTILE.su 
CNS054 . PS2 . SEXTILE . FALL 
CN5054 . PS2 . SEXTILE .WINTER 
OUTPUT  FILES 

DSNAME  CONTENTS 

CN5054.PS2.X13SP  Table  containing  13  X values  correspond- 

CN5054.PS2.X13SU  ing  to  each  route. 

CN5054 . PS2.  X13FALL 
CN5054 . PS2 . X13WXOTSR 
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MAJOR  VARIABLES 


M{J) 

XTEMP 

X (13) 


PROGRAM  STRUCTURE 
LINES 
1/5 

5/15 

16/34 


35/last 


Coefficients  c*  ^,c  spline 

X value  genera  iter  one  iteration 

of  the  Newton's  method 

The  13  X values,  six  of  which  are  the 

original  sextile  points  normalized 

between  [-1,1] 


FUNCTION 

Dimensions  various  vectors  and  reads  in 
a set  of  sextile  points 
Normalizes  the  sextile  points  and  checks 
if  they  are  monotonic. 

Determines  the  coefficient  of  the  cubic 
spline  from  the  sextile  points,  i.e.,  the 
mj ' a. 

Performs  two  iterations  of  the  Newton 
method  and  routes  the  generated  values 
plus  the  original  ones  on  a file. 


M.  PROGRAM;  X13  TABLE 
REFERENCE  PARAMETERS: 

DSNAME  = CN50M.PS2.X135U,CN5054.PS2.X13SP 

CN5054 . PS2 .X13FALL,CN5054 .PS2 .X13WINTER 
VOL  = S2R  •■=  PUBL03 
UNIT  = 3330 

DCB  - ( RECFM  = FB,  LRECL  » 80,BLKSIZE  - 3120) 
tOURCE:  Generated  by  THETAl 

USAGE:  Used  by  THETA  for  distribution  of  shipping  data  by  one 

degree  square, 

FORMAT:  Matrix  30x13  or  29x]3  c wording  to  the  season. 

Each  group  spans  two  lines  7 entries  are  on  the  first 
line  in  (7F10.5)  and  six  entries  on  the  subsequent  line. 
For  each  route  there  is  a corresponding  set  of  ,!3x 
values. 

Used  uy  p mm  THETA  in  the  distribution  of  the 
density  across  the  route. 


FUNCTION: 


N.  PROGRAM:  SEXTILE  TABLE 

REFERENCE  PARAMETERS: 

DSNAME  = CN5054.PS2, SEXTILE. SU, 

CN5054. PS2. SEXTILE. SP,CN5054.PS2. SEXTILE. FALL, 
CN5054.SP2. SEXTILE. WINTER 
VOL  = SER  = PUBL03 
UNIT  ® 3330 

DCB  - (RECFM  = FB,  LRECL  = 80,BLKSIZE  » 3120) 

SOURCE:  Ocean  Routes  Envelope  Report 

USAGE:  Used  by  program  THETA1  to  build  13  X values  for  each 

set  of  sextile  points 

FORMAT:  SEXTILE  table  7x29  matrix  or  7x30.  Varies  according 

to  seasons.  Each  group  spans  I line  with  7 entries  in 
format  (7F10.5). 

FUNCTION:  Used  by  THETA1  to  generate  the  13  X values. 


o 


PROGRAM:  THETA 


REFERENCE  PARAMETERS  - SOURCE  PROGRAM; 

DSNAME  » CN5054 .PS2 .THETAW 
VOL  = SER  ■ PUBL03 
UNIT  = 33330 

DCB  = (RECFM  * FB,  BL.KSIZE  = 3120,  LRECL  * 80) 
REFERENCE  PARAMETERS  - JCL 

DSNAME  = CN5Q54.PS2.THETA.JCLW 


VOL  = SER  ■ PUBL03 
UNIT  = 33330 

DCB  = (RECFM  = FB,  BLKSISE  = 3120,  LRECL  = 80) 
FUNCTION:  Takes  the  Ship-Days  - Table,  created  by  Epsilon, 

applies  the  growth  factor  to  the  ship  populations 
and  converts  the  data  from  ship-days  per  route  to 
ships  per  one  degree  square.  The  distributions 
are  done  between  canonical  ports,  for  the  various 


types  of  vessels  for  each  month. 

INPUT  DATA  SETS: 

DSNAME  CONTENTS 

CN5054 . PS1 . SGLMNDX  Inverse  Siglum  Table  & Siglum  Map 
CN5054 ,PC1 . ENVINDX  Index  Envelope  by  Route  and  Index  Envelope 


by  Siglum 

CN5054.PS2.X13SP  The  13  X values  generated  by  sextilepoints  per  season 
CN5Q54  .PS). SHP.DAYS  48  Ship  days  table,  per  month 

CN5054.PS2.POINTR.SP  Pointer  to  the  13  X values  by  route,  per  season 
CN5054.PS1.PRTDIST. SPRING  Port  Distance  Table,  one  for  each  season 


OUTPUT  DATA  SETS 

DSNAME  CONTENTS 

CN5054. KERCH. COUNTS. ARRAY. SPRING  Distribution  by  1°  square 

FALL 

WINTER 

SUMMER 
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PROGRAM  STRUCTURE 


LINES 

1/11 


11/17 

18/49 

50/55 

56/62 

63/67 


68/87 


88/91 

92/99 


FUNCTION 

Dimensions  various  arrays.  That  ia  SSIPDS 
for  the  ship  days  table,  PRTDST  the  port  dis- 
tance table,  SPLINE  the  Spline  Table,  SIGMP2 
the  Siglum-map,  INSIG  the  Index  Envelope  by 
Siglum,  IEBR  the  .index  Envelope  by  Route. 

The  Route  Trees  are  read  route  at  a time  and 
each  component  is  stored  in  a vector.  They 
consist  of  LENGTH,  the  edge  length,  LAT,  the 
latitude  of  the  vertics,  LONG,  longitude  of 
the  vertice,  BPOINT,  the  backpointer,  IDENTS, 
the  identifier  of  the  vertices,  WIDTH,  the 
vertex  width,  EDWITH,  the  edge  width  of  the 
segment,  DEN,  an  array  that  has  the  1°  square 
distribution. 

Initializes  the  distribution  matrix  DEN  to  zero. 
Reads  in  all  the  input  data  sets. 

Applies  the  growth  factor  to  the  ship  days 
table 

Begin  to  process  the  routes  one  at  a time. 

Checks  the  line  number  of  the  present  route 
by  use  of  a counter  called  ISTART  to  that  of 
the  Index  Envelope  by  Route. 

Checks  to  make  sure  that  a route  does  not 
exceed  41  lines.  It  then  proceeds  to  read  in 
the  designated  number  of  line  of  the  route. 

Sets  Flag2  to  False  if  the  route  crosses  the 
date  line. 

Accesses  the  appropriate  set  of  13  X values 
for  the  given  route. 


100/108 


109/110 


111/119 


120/145 


146/155 


156/164 

165 

166/171 

172/183 


184/408 


Looks  at  the  first  number  of  the  Sigluro  Map 
for  the  given  route.  This  indicates  the  number 
of  Siglia  assigned  to  that  route.  It  checks 
to  make  sure  that  it  is  not  zero  and  if  it  is 
it  processes  the  next  route. 

It  refers  to  the  Inder  Envelope  by  Siglum  for 
the  given  Siglia.  This  tells  us  the  line  numbers 
in  the  Route  Trees  for  the  two  canonical  ports. 
Next  it  computes  the  density,  i.e.,  the  number 
of  ship  days  per  month.  Also  it  determines  the 
distance  for  the  given  Siglia, 

Next  it  determines  the  peak  or  the  common  branch 
in  the  trees  by  doing  a comparison  on  the  identi- 
fiers. 

Calls  the  subroutine  to  actually  distribute 
the  number  of  ships  across  the  segment  of  the 
route. 

Writes  out  the  180x360  matrix  on  a data  set. 

The  subroutine  that  distributes  the  number  of 
ship  days  across  the  route. 

Dimensions  and  declares  various  variables. 

Checks  if  this  route  crosses  the  date  lines 
and  sets  a flag.  If  it  does  another  flag  is 
set  up  to  determine  which  segment  actually 
crosses  the  date  line. 

Processes  the  first  segment  according  to  the 
various  cases.  If  the  vertex  widths  of  the 
2 vertices  are  zero,  one  set  of  calculations 
are  performed.  If  one  vertex  is  zero  and  one 
is  non-zero,  one  will  be  sent  to  another  part 
of  the  program  for  various  calculations.  Still 
another  set  of  equations  are  used  when  both 
vertices  have  a non-zero  vertex  width.  All 
of  these  equations  can  be  obtained  from 

Chapter  4. 
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DATA  FILES 


1.  TABLE:  Routing  Table  (Alpha)  (Final  Version) 

REFERENCE  PARAMETERS : 

DSNAME  = CN5054 .PSl.RTNGTBL.NEW 
VOL  * SER  - PUBL03 
UNIT  * 3330 

DCB  « (RECFM=FB,LRECL=72,BLKSIZE=3096  ) 

SOURCE:  Generated  by  Program  Alpha  from  the  Original  Routing 

Table  and  the  Mediterranean  additions,  with  later 
corrections  made  simultaneously  to  the  original. 

USAGE:  Used  by  program  Betal  to  generate  the  siglum  Table. 

FORMAT:  1.  69x69  matrix  entered  in  69  groups  of  69  elements. 

Each  group  spans  3 lines  with  24  elements  each  on 
lines  1 and  2,  and  21  elements  on  line  3.  Each 
entry  contains: 

(a)  A 1-letter  code;  specifying  meanings  of  the 
following  number: 

P Canonical  Port 

Blank  Canonical  Route 
B Multiple  Routing  Table  Index 

N Null  Entry 

(b)  A 2-digit  number,  with  meaning  as  specified 
by  (a) . 

2.  Multiple  Routing  Table  15x2  matrix,  entered  in 
15  lines  of  2 entries  each.  Each  entry  contains  a 
2-digit  code,  denoting  a canonical  port. 

FUNCTION:  Let  a,  8 denote  two  canonical  ports  with  l<a, 

B569.  Then  the  (a, 8)  element  of  the  Routing  Table 
determines  the  route  from  port  a to  port  8 by 
specifying  either  an  intermediate  port  (denoted 
by  P) , a canonical  route  linking  a to  8 (denoted 
by  a blank) , a multiple  routing  table  to  be  checked 
(denoted  by  a B) , or  identity  of  ports  a and  8 


(denoted  by  an  N) . If  a multiple  routing  table  is 
specified,  let  y be  the  numerical  code  given.  Then 
the  (y,l)  and  (y,2)  entries  in  the  multiple  routing 
table  gives  the  intermediate  port  for  each  of  the  two 
options,  supertanker  and  non- supertanker. 

CORRECTION  PROCEDURE : 

Either  enter  changes  simultaneously  in  the  final  version, 
original  version,  and  Mediterranean  supplement  of 
Routing  Table,  or  enter  changes  in  the  original  version 
and  Mediterranean  supplement,  then  rerun  program  ALPHA. 
After  any  changes,  the  following  programs  must  be 
rerun:  GAMMA1.  After  changes  affecting  ports 

linked  by  direct  routes,  the  following  programs 
must  be  rerun:  BETA1,  BETA,  BETA2,  GAMMA2.  If 

these  modifications  cause  the  number  of  Sigla  to 
exceed  506,  programs  BETA  and  BETA 2 will  need  to 
be  modified  internally  as  well. 


2. 


TABLE: 


Routing  Table  (Alpha)  (Original  Version) 


REFERENCE 

SOURCE: 

USAGE: 

FORMAT: 

FUNCTION: 

CORRECTION 


PARAMETERS: 

DSNAME  = CN5054.PS1.RTNGTBL.ORG 
VOL  * SER  - PUBL03 
UNIT  » 3330 

DCB  = (RECFM=FB,LRECL=60,BLKSIZE=3120) 

Entered  from  the  Ocean  Routes  Envelopes  Report,  with 
corrections  (attached) . 

Used  by  program  Alpha  to  build  the  final  version  of  the 
Routing  Table. 

1.  Routing  Table  50x50  matrix,  entered  in  50  groups 
of  50  entries  each.  Each  group  spans  3 lines  with 
20  entries  each  on  lines  1 and  2 and  10  entries  on 
line  3.  Otherwise,  like  Routing  Table  (final 
version) . 

2.  Multiple  Routing  Table  4x2  matrix,  arranged  as  . 
in  the  Routing  Table  (final  version) . 

As  in  the  Routing  Table  (final  version) , except  that  the 
Mediterranean  canonical  ports  are  excluded. 

STATUS  PROCEDURE: 

See  Routing  Table  (final  version) . 
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3. 


TABLE: 


Routing  Table  (Alpha)  (Mediterranean  Additions) 


REFERENCE  PARAMETERS: 

DSNAME  = CN5054 .PSl.RTNGTBL.ADD 
VOL  = SER  * PUBL03 
UNIT  = 3330 

DCB  = (RECFM*FB,LRECL=72,BLKSIZE=3096) 

SOURCE:  Entered  from  the  Ocean  Routes  Envelopes  report,  with 

corrections  (attached) . 


USAGE:  Used  by  program  Alpha  to  build  the  Routing  Table 

(final  version). 

FORMAT:  1.  Routing  table  24x24  matrix,  entered  in  24  lines 

of  24  entries  each. 


FUNCTION:  It  is  like  the  Routing  Table  (final  version),  but  it  in- 
cludes only  Mediterranean  canonical  ports.  Indexing  of  the 

port  numbers  differs  from  that  of  final  version, 
and  is  related  to  it  by  the  Mediterranean  port/ 
route  renumbering  table. 


CORRECTION  STATUS  PROCEDURE: 

See  Routing  Table  (final  version) . 
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4.  TABLE:  Mediterranean  Port/Route  Renumbering  Table 

REFERENCE  PARAMETERS: 

DSNAME  = CN5054.PS1.WYLIB  (PORTNO) 

VOL  = SER  * WYLIB  1 
Wylbur  Edit  Format 

SOURCE:  Generated  by  program  Alpha  while  building  the  final 

Routing  Table. 

USAGE: 

port  numbers  assigned  to  the  Mediterranean*  s 
additional  ports  and  routes. 

FORMAT:  Lists  of  new  canonical  port,  route  number  and 

the  corresponding  Mediterranean  additional  port 
numbers  and  routes,  as  used  in  the  ORE  report. 


5. 


TABLE: 


Route  i:  ;ees  (Original  Form) 


REFERENCE 

SOURCE: 

USAGE: 

FORMAT : 


PARAMETERS: 

DSNAME  = CN5054.PS1. ENVELOPE 
VOL  = SER  = PUBL03 
UNIT  = 33*  ' 

DCB  - (REC^-«FB,LRECL=80,BLKSIZE=3120) 

Sequenced  in  columns  73  through  80. 

Entered  from  the  ORE  report,  with  subsequent  corrections. 

Used  by  program  Beta  to  generate  the  final  form  of 
Route  Trees  and  various  indices. 

84  groups  containing  a variable  number  of  lines, 
corresponding  to  each  of  the  84  canonical  routes. 

Each  group  is  headed  by  a marker  card  with  an'R' 
in  column  1 and  the  2-digit  route  number  in  columns 
3/4.  Each  line  within  the  group  represents  a separate 
vertex,  except  when  alternate  data  must  be  entered 
for  the  various  seasons.  In  this  case,  a group 
of  four  lines  represents  a vertex,  with  the  season 
code  ( * P * for  spring,  'S'  for  summer,  'F'  for  fall, 

'W'  for  winter)  punched  in  column  1.  Data  fields 
on  each  line  and  formats  are  as  follows: 

Columns 

1 Route  Header  or  season  marker 

2/5  Edge  width  (F4.2) 

6/10  Latitude  (2-digit  degree,  2-digit  minute, 

1-letter  hemisphere  code  (N  or  S) 

11/16  Longitude  (3-digit  degree,  2-digit  minute, 

1-letter  hemisphere  code  (W  or  E) 

17/20  Vertex  width  (F4.2) 

21  Always  blank 

22/23  If  the  vertex  is  identical  with  a canonical  port, 
the  canonical  port  number  was  inserted,  other- 
wise it  remained  blank 
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Columns 


FUNCTION: 


CORRECTION 


CORRECTION 


24  Always  Blank 

25/48  Binary  identifier 

49  If  the  vertex  is  identical  with  a canon- 
ical port/  a direction  code,  as  follows: 
Blank  Entrance  and  Egress  point 

1 Entrance  point  only 

2 Egress  point  only 

It  provides  the  basic  tree  structure,  which  was 
obtained  from  the  ORE  Report.  It  is  used  to 
create  the  final  fora  of  the  Route  Trees  by 
program  Beta. 

STATUS 

Errors  detected  by  BETA  program  corrected. 
Errors  detected  by  examining  output  from  BETA2 
and  DELTA1  have  been  corrected. 

PROCEDURE: 


Enter  corrections,  then  rerun  programs  BETA, 
BETA2  and  GAMMA2, 


6. 


TABLE: 


Route  Trees  (Final  Form) 


VERSIONS: 

REFERENCE 


SOURCE: 

USACS: 

FORMAT: 


FUNCTION : 


4 (one  for  each  season) 

PARAMETERS: 

DSNAME  » CN5054.PS1.OUTENV. SPRING,  CN5054.PS1 

OUTENV. SUMMER,  CN5054 .PSl.OUTENV. FALL, 

CN5054. PSl.OUTENV.  WINTER 

VOL  = SER  = PUBL03 
UNIT  = 3330 

DCB  * (RECFM=FB,  LR£CL»55,  BLKSIZE«3135) 

Generated  by  program  Beta  from  Route  Trees 
(original  version) 

Used  by  program  BETA2  to  build  the  PORT  DISTANCE  TABLE,  and 
also  in  Theta  to  distribute  the  ships  across  the  routes. 

Approximately  800  lines,  each  containing  data  for 
a single  vertex  as  follows: 

Columns 

1/24  Vertex  Identifier  (Binary) 

25/30  Vertex  Latitude  (F6.2)  (Decimal  Degrees) 

(Negative  - Southern  Hemisphere, 

Positive  - Northern  Hemisphere) 

31/37  Vertex  longitude  (F3.2)  (Decimal  Degrees) 

(Negative  - Eastern  Hemisphere  , 

Positive  - Western  Hemisphere) 

38/41  Vertex  Width  (F4.2) 

42/47  Edge  Length  (F6.2)  (la  Latitudinal  Degrees) 

42/51  Edge  Width  (F4.2) 

52/55  Back  Pointer  (14) 

Determine  the  Siglum  assigned  fco  a given  pair  of  canonical 
ports  from  the  Siglum  Table.  Look  up  this  Siglum  on 
the  INDEX  TO  ENVELOPE  BY  SIGLUM.  The  two  line 
numbers  contained  in  this  table  point  to  the  ver- 
ticiea  in  the  ROUTE  TREE  ENVELOPES  corresponding  to 
these  ports.  From  each  of  these  vertices,  trace 
the  Route  Trees  backwards  using  the  back  pointer 
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until  the  first  common  parent  node  is  encountered* 
This  may  be  most  easily  determined  by  comparing  the 
identifiers.  Alternatively,  use  the  INDEX  TO 
ENVELOPE  BY  ROUTE  to  locate  the  root  node  of  the 
tree  structure  for  any  given  route.  Look  up  the 
Siglum  Map  as  well,  using  the  route  number.  The 
list  of  sigla  given  here  for  that  route  may  now 
be  used  to  determine  the  terminal  leaves  of  the 
Route  Tree  by  looking  them  up  in  the  INDEX  TO 
ENVELOPE  BY  SIGLUM. 

CORRECTION  STATUS: 

As  for  Route  Trees  (original  version) . 


CORRECTION  PROCEDURE: 

Make  all  changes  to  the  original  version  of  the  Route 
Trees  (or  to  Routinq  Table) , then  rerun  BETA  to  update 
indices  as  well.  The  following  programs  must  then 
be  rerun:  BETA2,  GAMMA2. 


f 


7.  TABLE:  Siglum  Table 

REFERENCE  PARAMETERS: 

DSNAME  « CN5054 . PS1 . SGLMTBL 
VOL  = SER  = PUBL03 
UNIT  = 3330 

DCB  • (RECFM=FB,  BLKSIZE=3096 , LRECL=72) 

SOURCE:  1.  Siglum  Table:  constructed  by  .program  BETA1 

from  final  version  of  RTNGTBL. 

2.  Multiple  Routing  Table:  carried  from  Routing 

Table,  final  version. 

USAGE:  Will  be  used  in  EPSILON 

FORMAT:  1.  Siglum  Table:  69x69  matrix,  corresponding  to 

each  possible  pairing  of  world  canonical  ports, 
entered  in  69  groups  of  6S  entries  each.  Each  group 
spans  4 lines,  with  18  entries  in  each  of  the  first 
3 lines  and  15  in  the  last.  Each  entry  consists  of: 

a.  A 1- letter-code,  determining  meaning  of  the 

accompanying  number: 

P Canonical  Port 

Blank  Siglum 

B Multiple  Routing  Table  Index 

N Null  entry 

b.  A 3-digit  number,  with  meaning  as  specified 
by  (a) . 

2.  Multiple  Routing  Table:  See  Routing  Table, 

final  version. 

FUNCTION:  Like  the  Routing  Table,  except  that  route  numbers  have 

been  replaced  by  sigia.  To  every  ordered  pair  of 
canonical  ports  linked  by  a canonical  route  a 
different  siglum  is  assigned.  This  siglum  may 
be  used  later  to  retrieve  the  route  trees,  and  ship 
table. 

CORRECTION  STATUS: 

As  for  Routing  Table. 
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CORRECTION  PROCEDURE 


Make  all  corrections  to  routing  table,  then  rerun 
BETA1  to  update  indices  as  well.  The  following 
procedures  must  also  be  rerun:  BETA,  BETA2 , GAMMA2. 


8. 


TABLES:  Siglum  Indices 


REFERENCE 

SOURCE: 

USAGE: 

FORMAT: 


PARAMETERS: 

DSNAME  = CN5054.PS1.SGLMNDX 
VOL  = SER  = PUBL03 
UNIT  = 3330, 

DCB  a (RECFM=FB,  LRECL=72,  BLKSIZE=3096) 

Built  by  program  BETA1  while  constructing  the 
Siglum  Table  • 

Used  by  program  BETA  while  constructing  route 
envelope  tree  indices,  and  by  GAMMA  while  check- 
ing them. 

1.  Inverse  Siglum  Table:  Headed  by  word  'INV# 

SGLM#TBL' , followed  by  a list  of  number  pairs,  with 
12  pairs  per  line  and  2 spaces  between  pairs.  Each 
pair  corresponds  to  a successive  siglum,  and  contains 
two  2-digit  numbers  specifying  the  pair  of  cano- 
nical ports  (origin  and  destination,  in  order) 
associated  to  that  siglum.  (WARNING:  The  current 

version  of  BETA  expects  up  to  506  entries  in  the 
INV# SGLMtTBL , corresponding  to  a maximum  of  506 
sigla.  If  this  limit  is  exceeded,  BETA  must  be 
modified.  If  less  than  506  sigla  exist,  INV# 

SGLMtTBL  may  need  to  be  padded  with  blank  lines 

to  simulate  this  number.) 

2.  Siglum  Map:  Headed  by  word  'SGLMIMAP* . Con- 

tains 85  canonical  routes,  each  requiring  three 
lines.  Each  group  contains  a header  entry  plus 

45  entries  (15  per  line)  in  14  format  beginning  in 
Column  3.  The  header  entry  is  in  12  format  and 
begins  in  Column  1 of  the  first  line  of  each  group. 
The  header  entry  specifies  the  number  of  port  pairs 
linked  by  that  canonical  route  (or  equivalently, 
the  number  of  sigla  assigned  to  each  route). 
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These  sigla  themselves  follow  as  a list  in  the 
regular  entries,  with  excess  entry  positions 
set  to  0. 

FUNCTION:  Let  x be  a given  siglum  number.  Then  the  x-th  entry 

in  the  inverse  siglum  Table  gives  the  port  pair 
associated  with  that  siglum.  (If  this  port  pair 
were  looked  up  in  the  siglum  Table,  the  corres- 
ponding entry  there  would  be  x).  Let  y be  a 
canonical  route  number.  Then  the  y-th  group  of 
entries  in  the  siglum  Map  gives  all  sigla  corres- 
ponding to  that  route. 

CORRECTION  STATUS: 

As  for  Routing  Table. 

CORRECTION  PROCEDURE: 

Make  all  changes  to  Routing  Table,  then  rerun 
BETA1,  to  update  Siglum  Table  as  well.  Following 
programs  must  then  also  be  rerun:  BETA,  BETA2, 

GAMMA1, 
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TABLE:  Envelope 


REFERENCE 

SOURCE: 

USAGE: 

FORMAT: 


FUNCTION: 


CORRECTION 


PARAMETERS: 

DSNAME  = CN5054 . PS1 . ENVINDX 
VOL  = SER  * PUBL03 
UNIT  * 3330 

DCB  = (RECFM=FB,  LRECL=72,  BLKSIZE«3096) 

Built  by  program  BETA  from  Route  Envelope  Drees 
and  Siglum  Indices. 

Used  by  program  BETA2  when  building  Port  Distance 
Table,  and  by  GAMMA  and  THETA  whenever  referencing 
Route  Trees. Table.  Also  for  debugging  Route  Trees  and 
Routing  Table. 

1.  Index  to  Envelope  by  Route  headed  by  word 

' INDX#ENVLP#BY#ROUTE* , followed  by  84  numbers  in 
(14, 2X)  format  with  12  numbers  per  line,  forming 
7 lines. 

2.  Index  to  Envelope  by  Siglum  immediately  follows 
index  to  envelope  by  route.  Headed  by  word  'INDX# 
ENVLP # BY # SGLM ' followed  by  one  pair  of  numbers  for 
each  siglum  in  (14,  14,  4X)  format  with  6 pairs 
per  line.  (There  are  currently  501 

sigla,  hence  501  such  pairs). 

See  Route  Trees  (final  version).  A *0'  anywhere 
except  at  the  end  of  the  Index  to  Envelope  by 
Siglum  indicates  that  a leaf  was  missing  or  not 
correctly  identified  as  a canonical  port  in  the 
tree  structures  (original  version) , or  that  an 
incorrect  route  was  specified  in  the  routing 
table.  (A  0 at  the  end  occurs  if  the  number  of 
sigla  < 506) . 

STATUS: 


As  for  Route  Trees  and  Routing  Table 


I-'  - 


I 


CORRECTION  PROCEDURE: 

Make  all  changes  in  original  version  of  tree  structure 
or  in  routing  table.  Then  rerun  BETA1  (for  changes 
in  routing  table)  and  BETA  (always) . Before  re-run- 
ning BETA,  pad  the  slglum  index  with  blank  lines,  if 
necessary  to  simulate  506  sigla.  Following  programs 
must  then  also  be  rerun:  BETA2 , GAMMA2. 

L 


: 
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TABLE: 


Port  Distance  Table 


VERSIONS: 

REFERENCE 

SOURCE: 

USAGE: 

FORMAT: 

CORRECTION 

CORRECTION 

FUNCTION: 


4 (one  for  each  season) 

PARAMETERS: 

DSNAME  - CN5054.PS1.PRTDIST. SPRING 
VOL  » SER  » PUBL03 
UNIT  * 3330 

DCB  « (LRECL=72 , BLKSIZE=»3036,  RECFM«FB) 

Generated  by  BETA  2 program  from  route  trees 
(final  version,  copy  appropriate  to  given  season) 
and  Index  to  Envelope  by  Siglum. 

Used  by  GAMMA2  to  construct  estimated  travel  time 
table,  and  also  in  Epsilon  to  build  the  Ship  Days  Table. 

One  pair  of  numbers  in  (IX, 14, IX, F6. 2)  format  for  each 
siglum  (currently  approximately  500  in  all) , with 
6 pairs  in  a line.  The  first  number  in  each  pair 
is  the  number  of  the  siglum  to  which  the  pair 
corresponds,  the  second, the  distance  (in  units  of 
latitudinal  degrees)  between  the  pair  of  canonical 
ports  corresponding  to  that  siglum. 

STATUS: 

As  for  Routing  Table  . 

PROCEDURE: 

Make  all  changes  in  Routing  Table  or  Route  Tress 
then  re-run  BETM , BETA,  and  BETA2.  Program  GAMKA2 
must  then  also  be  rerun  to  construct  new  Estimated 
Travel  Time  Table. 

Determine  siglum  corresponding  to  port  pair  linked 
by  direct  route,  look  up  in  Port  Distance  Table  to 
determine  distance  between  them.  Unreasonable 
distances  indicate  poor  routings  by  Routing  Table 
or  route  trees. 
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TABLE: 


Estimated  Travel  Time  Table 


VERSIONS : 
REFERENCE 

SOURCE: 

USAGE: 

FORMAT: 


FUNCTION: 


4 (one  for  each  season) 

PARAMETERS: 

DSNAME  = CN5054.PS1. TRAVTME . SPRING 
VOL  = SER  = PUBL03 
UNIT  « 3330 

DCB  = (R£CFM=FB,  BLKSIZE=3096 , LRECL=72) 

Built  by  Program  GAMMA2  from  Siglum  Table  and  Port 
Distance  Table. 

To  be  used  by  GAMMA  when  needed  to  reconstruct 
arrival/departure  dates. 

69x69  matrix  each  entry  corresponding  to  every 
possible  corresponding  port  pair,  whether  linked 
by  direct  route  or  not.  These  entries  are  arranged 
in  69  groups  of  69  entries  each,  with  each  group 
spanning  4 lines  with  18  entries  each  in  lines 
1-3  and  15  entries  in  line  4.  (The  Siglum  Table 
is  laid  out  similarly) . Each  entry  is  in  14  for- 
mat and  gives  either  (if  positive)  the  expected 
travel  time  in  integral  days  between  the  two  ports 
or  else  (if  negative)  an  index  to  the  Multiple 
travel  Time  Table.  This  is  followed  immediately 
by  the  Multiple  Travel  Time  Table,  consisting  of 
pairs  of  numbers  in  (14,  14)  format,  one  pair  per 
line,  giving  the  estimated  travel  times  for  each 
of  the  2 type  classes  (super-tanker,  non-super  - 
tanker) . 

Let  (x,y)  denote  a pair  of  canonical  ports.  Then 
the  (x,y)  entry  gives  the  estimated  travel  time 
between  these  to  ports,  unless  it  is  negative.  In 
that  case,  let  its  absolute  value  be  Z,  then  the 
Zth  pair  in  the  Multiple  Travel  Time  Table  gives 
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the  estimated  travel  time  between  those  two  ports 

for  each  of  the  type  classes.  Unreasonable  travel  j 

times  indicate  poor  Routing  Table  or  Route  Trees. 

Travel  time  estimates  assume  velocity  of  12  knots.  J 

CORRECTION  STATUS s j 

As  for  Routing  Table  and  Route  Trees.  L 

L 

CORRECTION  PROCEDURE: 

Make  changes  in  Routing  Table  and  Route  Trees.  ( 

Then  rerun  programs  BETA,  BETA1,  BETA2,  and  GAMMA2. 


[ 


r 


i: 


12.  TABLE:  Port  Table 

REFERENCE  PARAMETERS: 

DSNAME  = CN5054.PS1. NEWPORT. 

VOL  = SER  = PUBL03 
UNIT  = 3330 

DCB  = (RECFM=FB,  BLKSIZE=3120,  LRECL=80) 

Sequenced  73/80 

{Old  version:  DSNAME  « CN5054 . PS1 .PORTTBL. ) 

SOURCE:  Hand  corrections  to  output  from  the  ORE  report  program, 

which  used  DOT  port  table. 

USAGE:  Used  by  GAMMA1  to  build  port  Hash  Table,  and  by 

GAMMA  to  assign  ports  to  proper  canonical  ports. 

FORMAT:  The  line  number  of  each  entry  in  the  port  Table  is 

the  port  number  used  by  DOT.  The  first  field 
contains  a new  port  number,  which  is  usually  the 
same  but  differs  in  the  following  cases: 

1.  Port  abbreviation  occurs  with  several  spellings  - 
port  numbers  reset  to  assign  the  same  port  number  to 
all  variants. 

2.  Pert  name  ia  a generic  term  (such  as  East  Coast 
u.s.)  - port  number  is  flagged  with  a negative  sign. 

The  other  fields  give  latitude,  longitude  of  port  and 
assigned  canonical  port. 

CORRECTION  STATUS: 

Mediterranean  canonical  port  assignments,  enter 
variant  spellings  and  generic  names  flagged.  Errors 
discovered  by  DELTA2  are  corrected  in  NEWPORT,  with  generic 
names  not  flagged  in  PORTTBL. 

CORRECTION  PROCEDURE: 

Enter  corrections  directly,  then  rerun  GAMM&l. 
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13. 


TABLE: 


Port  Hash  Table  (Machine  Readable  Form) 


iV 


I 

r 


REFERENCE  PARAMETERS: 


SOURCE: 

USAGE: 


FORMAT: 


DSNAME  = CN5054.PS1.PRTHASH, 

VOL  = SER  = PUBL03 
UNIT  = 3330 

DCB  = (RECFM=F  , LRECL=10000) 

Built  by  program  GAMMA 1 from  port  table. 

Used  by  program  GAMMA  to  determine  the  port  number 
assigned  to  each  port  abbreviation. 

Machine  readable  only.  Each  bin  consists  of  10 
bytes,  as  follows: 

Byte  Contents 

1/6  Port  Abbreviation  (character) 

7/8  Port  number  (fixed  binary  (15)  ) 

9/10  Overflow  pointer  (fixed  binary  (15)  ) 

The  prime  area  currently  contains  2999  bins,  and 
the  overflow  area  contains  2001  bins. 


CORRECTION  STATUS: 


Does  reflect  flagging  of  variant  spellings  and 
generic  names  in  the  port  table. 

CORRECTION  PROCEDURE: 

Make  changes  (if  any)  in  port  ?nble,  then  rerun 
program  GAKHA1. 

FUNCTION:  Take  character  string  value  of  the  port  abbreviation, 

treat  it  as  a number,  then  divide  it  by  2999  (or  other 
prime  number,  it  this  is  changed),  use  this  as  bin  index . 
Compare  port  abbreviations,  if  they  are  the  same# 
then  the  port  number  is  as  given  ii  that  bln. 

Otherwise,  use  overflew  pointer  as  the  new  bin  index 
(to  overflew  area)  and  repeat  this  until  either  port 
abbreviation  is  found,  or  overflow  printer  * 0# 
in  which  case  port  abbreviation  is  not  in  port 
table. 


it 


in 

wit 


:'1  ?■ 


<*''  ' 

wt 


■yt 

5*  * 

W* 


m 

-V  V 


If 


I 

[ 


\ 

! 

) 
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(CONT) 13.  TABLE:  Port  Hash  Table  (People  - readable  form). 

REFERENCE  PARAMETER: 

DSNAME  * CN5Q54.PS1.HASHOUT 
VOL  * 3ER=  POBL03 
UNIT  * 3330 

DCB  * (RECFM  = FB,  LRECL  ■ 18,  BLKSIZE  * 1800) 
SOURCE:  Built  by  program  GAMMA  1 from  PORT  TABLE. 

USAGE  « People  - readable  form  of  PORT-HASH  TABLE  - 
no  other  function. 

FORMAT: 

Columns  Field 

1/6  Port  Abbreviation  (Character) 

7/8  Blank 


9/12  Port  Number  (14) 

13/14  Blank 

15/16  Overflow  Index  (14) 

First  2999  records  belong  to  prime  region,  remainder  to 
| [ overflow. 


j FUNCTION:  Same  as  PORT  HASH  TABLE  (machine- readable  form). 

| I CORRECTION  STATUS,  PROCEDURE:  Same  as  PORT  HASH  TABLE  (machine- 


14 


TABLE:  Sextile  Array 

REFERENCE  PARAMETERS: 


DSNAME  = CN5054.PS2.X13FALL»CN5054..PS2.X13W/ 
CN5054.PS2.X13SP,CN5054 .PS2.X13SU 
VOL  = SER  = PUBL  03 

DCN  » (RECFM  = FB,  LRECL  = 80,  BLKSIZE  = 3120) 

SOURCE:  A supplementary  program  entitled  Newton  which 

processed  the  sextile  points  created  by  the  ORE. 

USAGE:  Used  by  program  Theta  to  help  produce  the 

Merchant  Counts  Array 

FORMAT:  An  array  which  varies  in  sizes  from  29x13  to  31x31 

according  to  seasons.  The  first  group  of  thirteen 
elements  spans  two  lines.  The  first  line  has  7 
entries  in  F10.5  format  and  the  second  line  has 
6 elements  in  the  same  format. 

FUNCTION:  It  takes  the  Type  I'  Cubic  Splines  consisting  of 
seven  sextile  points  from  the  ORE  work.  These 
points  are  then  normalized  to  the  interval  §-1,1  ] , 
and  the-  five  intermediate  sextile  points  which  were 
determined  by  Newtons  method  will  describe  the  dis- 
tribution. This  data  is  stored  in  the  Sextile  Array 
and  is  accessed  via  the  Cubic  Spline  Pointer  Table. 
The  last  set  of  points  is  the  general  default  spline. 


li 
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TABLE:  Cubic  Spline  Pointer  Table 

REFERENCE  PARAMETERS: 

DSNAME  ■ CN5054.PS2.PGINTR.SP,CN5054 .PS2.P0INTR.  SUMMER 

CN5054. PS2.P0INTR. FALL, CN5Q54.PS2.POINTR.  WINTER 

VOL  = SER  = PUBLQ3 

DCB  = (RECFM  = FB,  LRECL  = 80,  BLKSIZE  = 3120) 

SOURCE:  Hand  generated  according  to  given  route 

USAGE;  Used  by  program  THETA  as  one  of  the  inputs  to 
help  create  the  Merchants  Count  Array 


FORMAT:  There  are  7 lines  with  12  numbers  on  each  line 

in  13, 3X  format 

FUNCTION:  This  table  contains  the  pointers  to  the  Sextile 


Array.  This  table  contains  85  entries,  one  for 
each  route,  one  per  season. 


TABLE:  Proto  Routing  Table 

REFERENCE  PARAMETERS: 

DSNAME  = CN5054 . PS1 . PR0T048 , CN5054 . PS1 . PR0T044, 

CN5054 . PS1.PROTO40 ,CN5054.PSl,PROTO35 » 
CN5054.PS1.PROTO31 ,CN5G54.PS1.PRQT027 , 

CN5054 .PS1. PR0T022 ,CN5Q54 .PS1 .PR0T018 , 

CN5054 . PS1. PR0TQ14 ,CN5054.PS1 .PRQTO09, 

CN5054 . PS1 . PROTO05 ,CN5054 . PS1 . PROTOOl 

VOL  = SER  = PUBL03 

DCB  = (REFM  = FB,  LRECL  = 80,  BLKSIZE  « 3520) 

SOURCE:  Builc  by  program  GAMMA  while  processing  the 

Lloyds  data  base. 

USAGE:  Used  by  program  EPSILON  to  produce  the  Ships 

Days  Table  by  month. 

FORMAT:  There  are  three  69x69  arrays  one  for  merchants, 

tankers  less  than  70,000gt,  and  tankers  greater 
than  70,000gt.  The  first  line  of  each  array 
indicates  which  month  is  being  processed  for 
example  WEEK  = 35.  The  second  line  indicates  the 
shiptype  for  example,  TAU  = 1.  The  69x69  array 
signifies  the  number  of  canonical  ports,  of  a 
given  shiptype  in  a given  month.  The  format  is  as 

follows:  line  1 contains  message  WEEK  = , 

line  2 has  TAU  = . The  69  groups  of  69  entries 

each  spans  five  lines  with  16  entries  each  in 
lines  1-4  and  5 entries  in  line  5.  Each  entry  is  in 
15  format.  Therefore  the  format  is  as  follows: 

4 (1615/)  , 515 . The  code  for  the  ship  types  and  months 
are  as  follows: 

TAU  = 1 signifies  that  the  following  table  consists 

of  the  number  of  ship  days  for  merchant  vessels 
between  pairs  of  canonical  ports 


TAU  = 2 signifies  that  the  following  table  conaiata 
of  the  number  of  ship  days  for  tankers  under 
70,000  gt«  between  pairs  of  canonical  ports 
TAU  = 3 {signifies  that  the  following  table  consists 
of  the  number  of  ship  days  for  tankers  over 
70,000  gt.  between  pairs  of  canonical  ports 


WEEK  » 

1 

represents 

the 

month 

of 

January 

WEEK  = 

5 

represents 

the 

month 

of 

February 

WEEK  = 

9 

represents 

the 

month 

of 

March 

WEEK  = 

14 

represents 

the 

month 

of 

April 

WEEK  * 

18 

represents 

the 

month 

of 

May 

WEEK  * 

22 

represents 

the 

month 

of 

June 

WEEK  = 

27 

represents 

the 

month 

of 

July 

WEEK  « 

31 

represents 

the 

month 

of 

August 

WEEK  » 

35 

represents 

the 

month 

of 

September 

WEEK  « 

40 

represents 

the 

month 

of 

October 

WEEK  = 

44 

represents 

the 

month 

of 

November 

WEEK  = 

48 

represents 

the 

month 

of 

December 

Let 

: (x,y)  denote 

a pair  of 

canonical  ports. 

Then  the  (x,y)  entry 

gives 

the  total  number 

of  days  that  ships  have  been  enrcute  between 

ports  x and  y. 


17.  TABLE:  Ship  Days  Table 

REFERENCE  PARAMETERS: 

PSNAME  ■ CN5054.PS1 . SHP . DAYS48 , CN5054. PS1. SHP.DAYS44 , 
CN5054.PS1.SHP.DAYS40,CN5054.PS1.SHP.DAYS35, 
CN5054.PS1 . SHP. DAYS 31 ,CN5 054 . PS1 . SHP.DAYS27, 
CN5054 .PS1 . SHP.DAYS22 ,CN5054. PS1. SHP. DAYS18, 
CN5054 . PS1 . SHP .DAYS14 ,CN5054 . PSl . SHP . DAYS09, 
CN5054 . PS1 .SHP . DAYS05 , CN5054 .PS1 . SHP . DAYS01 

VOL  » SER  = PUBL03 

DCB  = (RECFM  = FB,  LRECL  =72,  BLKSIZE  = 3096) 

SOURCE:  Built  by  program  EPSILON  from  the  Proto  Routing 

Table,  Siglum  Table  and  Routing  Table. 

USAGE:  Used  by  program  Theta  to  produce  the  Merchants 

Count  Array 

FORMAT:  The  format  of  the  Ship  Days  Table  is  similar  to 

that  of  the  Proto-Routing  Table.  There  are  again 
three  different,  ship  types;  TAU  = 1,  TAU  = 2,  TAU  =» 
3.  Each  class  is  composed  of  501  elements  and  is 
listed  according  to  siglia  i.e.,  canonical  ports 
that  are  linked  by  a direct  route.  Following  this 
lies  the  Residual  Table,  again  for  each  ship  type. 
Heading  the  table  is  the  specified  TAU,  Then  the 
500  element  vector  spans  63  lines  where  the  first 
62  lines  contain  eight  elements  in  F9.1  format 
and  line  63  has  five  entires  in  the  same  format. 
Then  there  is  a header  entitled  RESIDUAL  TABLE. 

The  following  69  entries,  corresponding  to  the  can- 
onical ports,  which  span  six  lines.  The  first  five 
lines  contain  12  entires  in  16  format  and  the  sixth 
line  has  nine  entries  in  16  format. 
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17. 


TABLE:  Ship  Days  Table 

REFERENCE  PARAMETERS: 

DSNAME  = CN5054. PS1. SHP.DAYS48 ,CN5054 . PS1.SHP.DAYS44 , 
CN 5054. PS1. SHP .DAYS 40, CN5 054 . PS1. SHP.DAYS35, 
CN5054.PSl.SHP.DAYS:tl,CN5054.PSl.SHP.DAYS27, 
CN5  054 . PS1 . oHP . DAYS  22 , CN5054 . PS1 . SHP . DAYS10 , 
CN5054 .PS1. 5HP .DAYS14 , CN5054 .PSl .SHP .DAYS09, 
CN5054.PS1.SHP„DAYS05,CN5054 .PS1.SHP.DAYS01 
VOL  = SER  « PUBL03 

DCB  = (RECFM  = FB,  LRECL  = 72,  BLKSI2E  = 3096) 

SOURCE:  Built  by  program  EPSILON  from  the  Proto  Routing 

Table,  Siglum  Table  and  Routing  Table. 

USAGE:  Used  by  program  Theta  to  produce  the  Merchants 

Count  Array 


L 

ii 

L 

Ij 

0 

0 


FORMAT:  The  format  of  the  Ship  Days  Table  is  similar  to 

that  of  the  Proto-Routing  Table.  There  are  again 
three  different  ship  types;  TAU  *»  1,  TAU  = 2,  TAU  = 
3.  Each  class  is  composed  of  501  elements  and  is 
listed  according  to  siglia  i.e.,  canonical  ports 
that  are  linked  by  a direct  route.  Following  this 
lies  the  Residual  Table,  again  for  each  ship  type. 
Heading  the  table  is  the  specified  TAU.  Then  the 
500  element  vector  spans  63  lines  where  the  first 
62  lines  contain  eight  elements  in  F9.1  format 
and  line  63  has  five  entires  in  the  same  format. 
Then  there  is  a header  entitled  RESIDUAL  TABLE. 

The  following  69  entries,  corresponding  to  the  can- 
onical ports,  which  span  six  lines.  The  first  five 
lines  contain  12  entires  in  16  format  and  the  sixth 
line  has  nine  entries  in  16  format. 
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FUNCTION: 


The  500  element  vector  represents  the  total 
number  of  ship  days  for  the  given  siglia.  The 
Residual  Table  stores  the  number  of  ships  to  be 
distributed  for  local  traffic  and  'non-route' 
traffic  about  the  canonical  ports. 
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MEMORANDUM  FOR  DISTRIBUTION  LIST 

Subj : DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

Ref:  (a)  SECNAVINST  5510.36 

Enel:  (1)  List  of  DECLASSIFIED  LRAPP  Documents 

1 . In  accordance  with  reference  (a),  a declassification  review  has  been  conducted  on  a 
number  of  classified  LRAPP  documents. 

2.  The  LRAPP  documents  listed  in  enclosure  (1)  have  been  downgraded  to 
UNCLASSIFIED  and  have  been  approved  for  public  release.  These  documents  should 
be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N772)  letter  N772A/6U875630,  20  January  2006. 

DISTRIBUTION  STATEMENT  A:  Approved  for  Public  Release;  Distribution  is 
unlimited. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


^ 

BRIAN  LINK 
By  direction 
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